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22nd  International  Symposium  on  Ballistics 
Vancouver  BC,  Canada 
14-18  November  2005 


Agenda 


Tuesday.  15  November  2005 

Invited  Presentation:  by  From  Columbia  To  Discovery:  Understanding  the  Impact  Threat  to  the  Space  Shuttle,  by  James  D.  Walker,  Southwest  Research  Institute 

General  Oral  Session  #1 

•  Plasma  Ignition  of  a  30mm  Cannon,  Richard  A.  Beyer,  Andrew  L.  Brant,  Joseph  J.  Colburn,  US  Army  Research  Laboratories 

•  Numerical  Computations  of  Subsonic  and  Supersonic  Flow  Choking  Phenomeana  in  Grid  Finned  Projectiles,  Nicolas  Parise,  SNC  Technologies,  Inc;  Alain 
Dupuis,  Precision  Weapons  Section,  Defense  Research  and  Development  Canada 

•  The  Use  of  Electric  Power  in  Active  Armour  Applications,  Martin  van  de  Voorde,  R.  Boeschoten,  TNO  Defence,  Security  and  Safety 

•  Prevention  of  Sympathetic  Detonation  Between  Reactive  Armor  Sandwiches,  Andreas  Holzwarth,  Fraunhofer-Institut  Fur  Kurzzeitdynamik 

Terminal  Ballistics  Oral  Session  #7 

•  Bullet  Impact  on  Steel  and  Kevlar®/Steel  Armor  -  Experimental  Data  and  Hydrocode  Modeling  with  Eulerian  and  Lagrangian  Methods,  Dale  S.  Preece, 
Vanessa  S.  Berg,  and  Loyd  R.  Payne,  Sandia  National  Laboratories 

•  Progress  on  the  NDE  Characterization  of  Impact  Damage  in  Armor  Materials,  Joseph  M.  Wells,  JMW  Associates 

•  The  influence  of  sabot  threads  on  the  performance  of  KE  penetrators,  Nick  Lynch  and  John  Stubberfield,  QinetiQ 

•  The  Use  of  Electric  Power  in  Active  Armour  Applications,  Martin  van  de  Voorde,  R.  Boeschoten,  TNO  Defence,  Security  and  Safety 

•  Prevention  of  Sympathetic  Detonation  Between  Reactive  Armor  Sandwiches,  Andreas  Holzwarth,  Fraunhofer-Institut  Fur  Kurzzeitdynamik 

Exterior  Ballistics  Oral  Session  #1 

•  Micro-Adaptive  Flow  Control  Applied  to  a  Spinning  Projectile,  Dr.  Jubaraj  Sahu,  U.S.  Army  Research  Laboratory 

•  Micro-Adaptive  Flow  Control  and  Nonlinear  Aerodynamics  (Combustion  Gas  Generators),  Dr.  Jubaraj  Sahu  and  Ms.  Karen  Heavey,  U.S.  Army  Research 
Laboratory 

•  Aerodynamic  Characteristics  of  a  Grid  Finned  Projectile  from  Free-Flight  Tests  at  Supersonic  Velocities,  Alain  Dupuis,  DRDC  -  Valcartier;  Claude  Berner, 
French-German  Research  Institute 

•  Recent  Computations  and  Validations  of  Projectile  Unsteady  Aerodynamics,  Roxan  Cayaz,  Eric  Carette,  Giat  Industires;  Remy  Thepot,  Patrick  Champigny, 
Office  National  D'Etudes  et  de  Recherches  Aerospatiales 

•  Research  of  Flight  Characteristics  of  Rod-Type  Projectile  with  Triangular  Cross-section,  Dr.  Wenjun  Yi,  Prof.  Xiaobing  Zhang,  and  Prof.  Jianping  Qian, 
Ballistic  Research  Laboratory  of  China 


Wednesday.  16  November  2005 


Exterior  Ballistics  Oral  Session  #2 

•  Impact  of  Nose-Mounted  Micro-Structures  on  the  Aerodynamics  of  a  Generic  Missile,  Dr.  Daniel  Corriveau,  Defence  R&D  Canada  (DRDC  -  Valcatier) 

•  Analyses  of  Gliding  Control  for  an  Extended-Range  Projectile,  Prof.  Zhongyuan  Wang,  Prof.  Houqian  Xu,  Dr.  Jinguang  Shi,  Dr.  Wenjun  Yi,  Prof.  Shaosong 
Chen,  Ballistic  Research  Laboratory  of  China 

•  Theoretical  Design  for  a  Guided  Supersonic  Projectile,  Pierre  Wey,  Claude  Berner,  Eckhart  Sommer,  Vollcer  Fleck,  Henry  Moulard,  French-German  Research 
Institute  of  Saint-Louis  (ISL) 

Interior  Ballistics/Launch  Dynamics  Oral  Session  #1 

•  Ceramic  Gun  Barrel  Technologies,  Larry  Burton,  Jeff  Swab,  Rob  Carter,  Ryan  Emerson,  U.S.  Army  Research  Laboratory,  Weapons  &  Materials  Research 
Directorate 

•  M865  TID  Improvement  Study,  Kerry  Henry,  Army  Research  Development  &  Engineering  Center;  Jason  W.  Gaines,  General  Dynamics-OTS 

•  Two-Dimensional  Modeling  of  Mortar  Interal  Ballistics,  Clive  R.  Woodley,  David  Finbow,  QinetiQ;  Vladimir  Titarev,  Eleuterio  Toro,  Umeritek  Limited 

•  The  Mechanism  Analysis  of  Interior  Ballistics  of  Serial  Chamber  Gun,  Dr.  Sanjiu  Ying,  Charge  Design  Laboratory  of  China;  Prof.  Xiaobing  Zhang,  Prof. 
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Yaxiong  Yuan,  and  Dr.Yan  Wang,  Ballistic  Research  Laboratory  of  China 
Terminal  Ballistics  Oral  Session  #2 

•  Behind  Armor  Debris  Computations  with  Finite  Elements  and  Meshless  Particles,  Gordon  R.  Johnson  and  Robert  A.  Stryk,  Network  Computing  Services,  Inc. 

•  Experimental  and  Numerical  Study  of  the  Penetration  of  Tungsten  Carbide  Into  Steel  Targets  During  High  Rates  of  Strain,  Eva  K.  Friis,  Nammo  Raufoss  AS; 
Oyvind  Froyland  and  John  F.  Moxnes,  FFI  (Norwegian  Defence  Research  Establishment) 

•  Mine  Neutralisation  with  Small  Calibre  Projetile  Impact,  Mark  Dijkstra,  J.H.  Meulman,  TNO  Defence  Safety  and  Security 

Vulnerability ,  Lethality  and  Wound  Ballistics  Oral  Session 

•  The  Application  of  Critical  Perforation  Analysis  (CPA)  to  Military  Personal  Armour  Research  and  Evaluation,  Catherine  H.  Crawford  and  Philip  Gotts, 
Defence  Clothing  Research  and  Project  Support 

•  Fragment  Patterns  Behind  Concrete  Structures  Caused  by  KE  Projectiles,  S.  Lampert,  Fene  Jeanquartier,  D.  Hoffmann,  and  B.  Lehmann,  Armasuisse 

•  Mine  Neutralisation  with  Small  Calibre  Projetile  Impact,  Mark  Dijkstra,  J.H.  Meulman,  TNO  Defence  Safety  and  Security 


Thursday.  17  November  2005 


Warhead  Mechanisms  Oral  Session  #/: 

•  Soft-Recovery  of  Explosively  Formed  Penetrators,  David  Lambert  and  Matthew  Pope,  Air  Force  Research  Laboratory,  Munitions  Directorate;  Stanley  E. 

Jones  and  Jonathan  Muse,  Aerospace  Engineering  and  Mechanics,  University  of  Alabama 

•  The  Gurney  Velocity:  A  “Constant”  Affected  by  Previously  Unrecognized  Factors,  Joseph  E.  Backofen,  BRIGS  Co. 

•  The  Influence  of  Post  Detonation  Burning  Process  on  Blast  Wave  Parameters  in  Air,  Meir  Mayseless  and  I.  Belsky,  IDF,  Armor  Branch;  E.  Muzychuk,  IMI, 
Central  Laboratory 

Interior  Ballistics/Launch  Dynamics  Oral  Session  #2 

•  3-D  Finite-Element  Gun  Launch  Simulation  of  a  Surrogate  Excalibur  155-mm  Guided  Artillery  Projectile  -  Modeling  Capabilities  and  its  Implications,  M.R. 
Chowdhury  and  A.  Frydman,  US  Army  Research  Laboratory;  J.  Cordes,  L.  Reinhardt  and  D.  Carlucci,  US  Army  ARDEC,  Analysis  and  Evaluation  Division 

•  Caseless  Ammunition  and  Advances  in  the  Characterization  of  High  Ignition  Temperature  Propellant  (HITP),  Paul  Shipley,  AAI  Corporation;  Erin  K. 
Hardmeyer,  US  Army  ARDEC;  and  Ben  Ashcroft,  Alliant  Technical  Systems 

•  Ballistic  Launch  to  Space,  Ed  Schmidt  and  Mark  Bundy,  Army  Research  Laboratory 

•  A  Novel  Launcher  for  Cavitating  Weapons,  Chris  Weiland  and  Pavlos  Vlachos,  Mechanical  Engineering  Department,  Virginia  Polytechnic  Institute  and  State 
University;  and  Jon  Yagla,  Engagement  Systems  Department,  Naval  Surface  Warfare  Center 

Warhead  Mechanisms  Oral  Session  #2 

•  The  Role  of  Rayleigh  Taylor  Instability  in  Shaped  Charge  Jets  Formation  and  Stability,  Dr.  Simcha  Miller,  Mr.  Gershon  Kliminz,  Rafael  Ballistic  Center, 

•  Simulation  of  Cylinder  Expansion  Tests  Using  an  Eulerian  Multiple-Material  Approach,  Laura  Donahue,  R.C.  Ripley,  Martec,  Ltd 

•  Application  of  Powder  Tantalum  Material  for  Explosively  Formed  Penetrator  (EFP)  Warhead,  Richard  Fong,  Mike  Hespos,  William  Ng,  and  Steven  Tang,  US 
Army  ARDEC 

•  Oilwell  Perforators:  Theoretical  Considerations,  Brenden  Grove,  Schlumberger  Reservoir  Completions  Center 

•  The  Study  on  Lethality  Simulation  Method  for  Fragmentation  Warhead,  Yang  Yunbin,  Qu  Ming,  and  Qian  Lixin,  Institute  of  Structual  Mechanics,  China 
Academy  of  Engineering  Physics 


Friday.  18  November  2005 


Terminal  Ballistics  Oral  Session  # 3 : 

•  Performance  Evaluation  of  Multi-Threat  Body  Armour  Systems,  B.  Anctil  and  M.  Keown,  Biokinetics  and  Associates  Ltd.;  G.  Pageau,  M.  Bolduc,  and  D. 
Bourget,  Defence  R&D  Candada  -  Valcartier 

•  The  Residual  Damage  in  CFRP  Composite  After  Ballistic  Impacts  (Experiments  &  Simulations),  Koen  Herlaar,  TNO  Defence,  Security  and  Safety 

•  The  Effect  of  Boundary  Conditions  on  the  Ballistic  Performance  of  Textile  Fabrics,  Colin  R.  Cork,  University  of  Manchester,  School  of  Materials 

General  Oral  Session  #2 

•  Wind  Tunnel  Verification  of  the  Performance  of  a  Smart  Material  Canard  Actuator,  Paul  Weinacht,  William  F.  Drysdale,  Travis  Bogetti,  and  Rod  Don,  US 
Army  Research  Laboratory;  James  T.  Arters,  Jack  R.  Vinson,  Aaron  R.  Hickman,  University  of  Delaware;  Lamar  Auman,  US  Army  Aviation  and  Missile 
RD&E  Center;  Oded  Rabinovitch,  Technion  Israel  Institute  of  Technology 

•  The  Fragmentation  of  Metal  Cylinders  Using  Thermobaric  Explosives,  William  Andrews,  Royal  Military  College  of  Canada;  Michael  Dunning,  Defence  R&D 
Canada  -  Suffield;  and  Kevin  Jaansalu,  Montana  Tech  (University  of  Montana) 

•  A  Novel  Test  Methodology  to  Assess  the  Performance  of  Ballistic  Helmets,  B.  Anctil  and  M.  Keown,  Biokinetics  and  Associates  Ltd.;  G.  Pageau  and  D. 
Bourget,  Defence  R&D  Candada  -  Valcartier 

•  Ballistic  Analysis  of  Bulgarian  Electroslag  Remelted  Dual  Hard  Steel  Armor  Plate,  William  Gooch,  Matthew  Burkins,  and  David  Mackenzie,  US  Army 
Research  Laboratory  Weapons  and  Materials  Research  Directorate;  Stefan  Vodenicharov,  Institute  of  Metal  Science,  Bulgarian  Academy  of  Sciences 
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SYMPOSIUM  SCOPE  AND  OBJECTIVES 

The  objective  of  the  22nd  International  Symposium  on  Ballistics  is  to  focus  on  potential  technical 
advances  and  break-throughs  in  the  21st  century  in  the  general  areas  of: 

Interior  Ballistics 

Launch  Dynamics 

Exterior  Ballistics 

Projectile  and  Warhead  Design 

Terminal  Ballistics 

Vulnerability 

Modeling  and  Simulation 
Wound  Ballistics 


Over  200  papers  will  be  presented  by  authors  from  26  countries. 
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SYMPOSIUM  PROGRAM 

Monday.  November  14.  2005 


2:00  pm  -  5:00  pm 

Registration 

5:00  pm  -  6:30  pm 

Reception  in  Exhibit  Area 

Tuesday.  November  15.  2005 


7:00  am 

Continental  Breakfast  and  Registration 

7:00  am  -  6:00  pm 

Exhibits  Open 

8:00  am 

Opening  Remarks 

C.  Samuel  Campagna ,  National  Defense  Industrial  Association 

8:10  am 

Welcome  and  Opening  Remarks 

William  Flis,  DE  Technologies,  Inc.  &  Brian  Scott,  US  Army  Research  Laboratory 

8:20  am 

Keynote  Address 

Dr.  Robert  Walker,  Director-General,  Research  and  Development  Programs 
(DGRDP),  Defense  Research  and  Development  Canada 

9:05  am 

Invited  Presentation 

From  Columbia  to  Discovery:  Understanding  the  Impact  Threat  to  the  Space  Shuttle 
James  D.  Walker,  Southwest  Research  Institute 

9:50  am 

Morning  Break 

General  Oral  Session  #1 

Chairpersons:  B.  Janzon  and  J.  Carleone 

10:20  am 

Plasma  Ignition  of  a  30mm  Cannon 

Richard  A.  Beyer ;  Andrew  L.  Brant ,  Joseph  J.  Colburn,  US  Army  Research 
Laboratories 

10:40  am 

Numerical  Computations  of  Subsonic  and  Supersonic  Flow  Choking  Phenomena  in 
Grid  Finned  Projectiles 

Nicolas  Parise,  SNC  Technologies,  Inc.;  Alain  Dupuis,  Precision  Weapons  Section, 
Defense  Research  and  Development  Canada 

11:00  am 

Multiple  Explosively  Formed  Penetrator  (MEFP)  Warhead  Technologies  for  Mine  and 
Improvised  Explosive  Device  (IED)  Neutralization 

Richard  Fong,  William  Ng,  Steve  Tang,  LaMar  Thompson,  U.S.  Army  Armament 
Research,  Development  and  Engineering  Center 

11:20  am 

The  Use  of  Electric  Power  in  Active  Armour  Applications 

Martin  van  de  Voorde,  R.  Boeschoten,  TNO  Defence,  Security  and  Safety 

11:40  am 

Prevention  of  Sympathetic  Detonation  between  Reactive  Armor  Sandwiches 

Andreas  Holzwarth,  Fraunhofer-lnstitut  fur  Kurzzeitdynamik 

12:00  pm 

Lunch 
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1:30  pm 

1:30  pm 

1:50  pm 
2:10  pm 

2:30  pm 

2:50  pm 

3:10  pm 
3:40  pm 

3:40  pm 
4:00  pm 

4:20  pm 

4:40  pm 

5:00  pm 

5:20  pm 


3:10  pm  Exterior  Ballistics  Poster  Session 

Chairpersons:  Z.  Wang  and  P.A.  Karsten 

Terminal  Ballistics  Oral  Session  #1 
Chairpersons:  E.  Linden  and  C.  Anderson 

Bullet  Impact  on  Steel  and  Kevlar®/Steel  Armor  -  Experimental  Data  and  Hydrocode 
Modeling  with  Eulerian  and  Lagrangian  Methods* 

Dale  S.  Preece,  Vanessa  S.  Berg,  Mathew  A.  Risenmay,  Sandia  National 
Laboratories 

Progress  on  the  NDE  Characterization  of  Impact  Damage  in  Armor  Materials 
Joseph  M.  Wells,  JMW  Associates 

Design,  Analysis,  and  Testing  of  an  Unconfined  Ceramic  Target  to  Induce  Dwell 
Timothy  J.  Holmquist,  Network  Computing  Services,  Inc.;  C.  Anderson,  Jr., 
Southwest  Research  Institute;  Thilo  Behner,  Ernst-Mach-lnstitut 

The  Influence  of  Sabot  Threads  on  the  Performance  of  KE  Penetrators  against 

multiple  plate  targets 

Nick  J.  Lynch,  J.  Stubberfield,  QinetiQ 

Visualization  of  Wave  Propagation  and  Impact  Damage  in  a  Polycrystalline 
Transparent  Ceramic  -  AION 

Elmar  Strassburger,  Fraunhofer  Institut  fur  Kurzzeitdynamik;  Parimal  Patel,  James 
W.  McCauley,  US  Army  Research  Laboratory;  Douglas  W.  Templeton,  US  Army 
TARDEC 

Afternoon  Break 

5:20  pm  Terminal  Ballistics  Poster  Session  #1 
Chairpersons:  A.  Diederen 

Exterior  Ballistics  Oral  Session 
Chairpersons:  W.  Reinecke  and  A.  Dupuis 

Advanced  Time-Accurate  CFD/RBD  Simulations  of  Projectiles  in  Free  Flight 
Jubaraj  Sahu,  US  Army  Research  Laboratory 

Aerodynamic  Characterstics  of  a  Grid  Finned  Projectile  from  Free-Flight  Tests  at 
Supersonic  Velocities 

Alain  Dupuis,  DRDC  -  Valcartier;  Claude  Berner,  French-German  Research 
Institute 

Recent  Computations  and  Validations  of  Projectile  Unsteady  Aerodynamics 
Roxan  Cayaz,  Eric  Carette,  Giat  Industries;  Remy  Thepot,  Patrick  Champigny, 
Office  National  d’Etudes  et  de  Recherches  Aerospatiales 

The  Derivation  of  Spin  Stabilised  Projectile  Yaw  Rates  and  Ballistic  Model 
Coefficients  Using  Conventional  CW  Doppler  Radar  Systems 

John  Tate,  FLEET 

Research  of  Flight  Characteristics  of  Rod-Type  Projectile  with  Triangular 
Cross-Section 

Wenjun  Yi,  Xiaobing  Zhang,  Jianping  Qian,  Ballistic  Research  Laboratory  of 
China,  Nanjing  University  of  Science  &  Technology 

Adjourn  for  the  Day 
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Wednesday.  November  16.  2005 


7:00  am 

Continental  Breakfast  and  Registration 

7:00  am  -  5:00  pm 

Exhibits  Open 

8:00  am 

Administrative  Remarks 

8:10  am  -  9:50  am 

Terminal  Ballistics  Poster  Session  #2 

Chairpersons:  J.  Riegel  and  E.  Hirsch 

Exterior  Ballistics  Oral  Session  #2 

Chairpersons:  P.  Nel  and  E.  Schmidt 

8:10  am 

Impact  of  Nose-Mounted  Micro-Structures  on  the  Aerodynamics  of  a  Generic  Missile 
Daniel  Corriveau,  Defence  R&D  Canada  (DRDC  -  Valcatier) 

8:30  am 

Ballistic  Simulations  and  Wind  Tunnel  Testing  of  120  mm  Mortar  Bomb  Tail  Fin 
Geometries  -  In  Search  for  Extra  Range 

Jukka  Tiainen,  Ari  Makkonen,  Patria  Weapon  Systems  Oy;  Mikko  Korhonen, 
Timo  Sailaranta,  TKK/Laboratory  of  Aerodynamics 

8:50  am 

Bringing  Solid  Fuel  Ramjet  Projectiles  Closer  to  Application  -  An  Overview  of  the 
TNO/RWMS  Technology  Demonstration  Programme 

Ronald  G.  Veraar ,  TNO  Defence,  Security  and  Safety  Research  Group  Rocket 
Technology;  Guido  Giusti,  Rheinmetall  Waffe  Munition  Schweiz  AG 

9:10  am 

Analysis  of  Gliding  Control  for  an  Extended-Range  Projectile 

Zhongyuan  Wang,  Houqian  Xu,  Jinguang  Shi,  Wenjun  Yi,  Shaosong  Chen, 
Ballistic  Research  Laboratory  of  China,  Nanjing  University  of  Science  &  Technology 

9:30  am 

Theoretical  Design  for  a  Guided  Supersonic  Projectile 

Pierre  Wey,  Claude  Berner,  Eckhart  Sommer,  Volker  Fleck,  Henry  Moulard, 
French-German  Research  Institute  of  Saint-Louis  (ISL) 

9:50  am 

Morning  Break 

10:20  am  - 12:00  pm 

Warhead  Mechanisms  Poster  Session 

Chairpersons:  R.  Fong  and  F.  Mostert 

Interior  Ballistics/Launch  Dynamics  Oral  Session  #1 

Chairpersons:  C.  Candland  and  C.  Woodley 

10:20  am 

Ceramic  Gun  Barrel  Technology 

Lawrence  W.  Burton,  Jeffrey  J.  Swab,  Ryan  Emerson,  Robert  Carter,  US  Army 
Research  Laboratory,  Weapons  &  Materials  Research  Directorate 

10:40  am 

M865E3  Cold  Target  Impact  Dispersion  Study 

Kerry  Henry,  Army  Research  Development  &  Engineering  Center;  Jason  W. 
Gaines,  General  Dynamics-OTS 

11:00  am 

An  Alternative  Technique  to  Evaluate  and  Characterize  Pressure  Waves  in  Large 
Calibre  Guns 

Victor  Schabort,  Denel  Land  Systems  Western  Cape 

11:20  am 

Two-Dimensional  Modelling  of  Mortar  Internal  Ballistics 

Clive  R.  Woodley,  David  Finbow,  QinetiQ;  Vladimir  Titarev,  Eleuterio  Toro, 
Numeritek  Limited 
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11:40  am 


12:00  pm 


1:30  pm 


1:50  pm 


2:10  pm 


2:30  pm 


2:50  pm 


3:10  pm 


3:40  pm 


4:00  pm 


4:20  pm 


4:40  pm 


The  Mechanism  Analysis  of  Interior  Ballistics  of  Serial  Chamber  Gun 
Sanjiu  Ying,  Charge  Design  Laboratory  of  China,  Nanjing  University  of  Science  & 
Technology;  Xiaobing  Zhang ,  Qaxiong  Yuan ,  Yan  Wang ,  Ballistic  Research 
Laboratory  of  China,  Nanjing  University  of  Science  &  Technology 

Lunch 


Terminal  Ballistics  Oral  Session  #2 
Chairpersons:  M.  Mayseless  and  T.  Holmquist 

Behind  Armor  Debris  Computations  with  Finite  Elements  and  Meshless  Particles 
Gordon  R.  Johnson ,  Robert  A.  Stryk,  Network  Computing  Services,  Inc. 

Expirmental  and  Numerical  Study  of  the  Penetration  of  Tungsten  Carbide  into  Steel 
Targets  During  High  Rates  of  Strain 

Eva  K.  Friis,  Nammo  Raufoss  AS,  Oyvind  Froyland,  John  F.  Moxnes,  FFI 
(Norwegian  Defence  Researdh  Establishment 

Fragmentation  Behavior  of  Tungsten  Alloy  Cubes  on  Normal  Aluminum  Plate  Targets 
Karl  Weber,  Fraunhofer-lnstitut  fur  Kurzzeitdynamik,  Ernst-Mach  Institut 

The  Failure  Kinetics  of  High  Density  DEDF  Glass  Against  Rod  Impact  at  Velocities 
From  0.4  to  2.5  km/s 

Thilo  Behner,  V.  Hohler ;  M.  Moll,  Fraunhofer  Institut  fur  Kurzzeitdynamik  (Ernst- 
Mach  Institut);  Ch.  E.  Anderson  Jr.,  Southwest  Research  Institute;  D.  L.  Orphal, 
International  Research  Associates,  Inc.;  D.  W.  Templeton,  US  Army 
RDECOM-TACOM 

Mine  Neutralisation  with  Small  Calibre  Projectile  Impact 

Mark  Dijkstra,  J.H.  Meulman,  TNO  Defence,  Safety  and  Security 

Afternoon  Break 

Vulnerability,  Lethality  and  Wound  Ballistics  Oral  Session 
Chairpersons:  R.  Vaziri,  A.  Persson 

The  Application  of  Critical  Perforation  Analysis  (CPA)  to  Military  Personal  Armour 
Research  and  Evaluation 

Catherine  H.  Crawford ,  Philip  Gotts,  Defence  Clothing  Research  and  Project 
Support 

An  Efficient  Mechanistic  Approach  to  Modelling  the  Ballistic  Response  of  Multi-Layer 
Fabrics 

Ali  Shahkarami,  Reza  Vaziri ,  Anounsh  Poursartip,  Composites  Group, 
Departments  of  Civil  Engineering  and  Materials  Engineering  The  University  of  British 
Columbia;  Navin  Tajani,  DuPont  Advanced  Fibers  Systems 

Pencilling -A  Novel  Behind  Armour  Blunt  Trauma  Injury 

Eluned  A.  Lewis ,  Defence  Clothing  Research  and  Project  Support;  Ian  Horsfall , 
Celia  Watson,  Engineering  Systems  Department,  Royal  Military  College  of  Science, 
Cranfield  University 

Scaling  the  Dynamic  Response  of  Armored  Vehicle’s  Floor  Subjected  to  a  Large 
Buried  Charge 

Avidov  Neuberger,  MOD,  Tank  Program  Management;  S.  Peles,  IMI,  Central 
Laboratory  Division;  D.  Rittel,  Technion,  Israel  Institute  of  Technology,  Faculty  of 
Mechanical  Engineering 


22nd  International  Symposium  on  Ballistics 


5:00  pm 

Fragment  Patterns  Behind  Concrete  Structures  Caused  by  KE  Projectiles 

Rene  Jeanquartier ;  D.  Hoffmann ,  S.  Lam  pert ,  B.  Lehmann,  Armasuisse 

5:20  pm 

Adjourn  for  the  Day 

Thursday.  November  17.  2005 


7:00  am 

Continental  Breakfast  and  Registration 

8:00  am  - 11:00  am 

Exhibits  Open 

8:00  am 

Adminstrative  Remarks 

8:10  am  -  9:50  am 

Interior  Ballistics/Launch  Dynamics  Poster  Session 

Chairpersons:  C.  Woodley  and  C.  Candlant 

Warhead  Mechanisms  Oral  Session  #1 

Chairpersons:  M.  Murphy  and  P.Y.  Chanteret 

8:10  am 

Soft-Recovery  of  Explosively  Formed  Penetrators 

David  E.  Lambert,  Matthew  Pope,  Air  Force  Research  Laboratory,  Munitions 
Directorate;  Stanley  Jones,  Jonathan  Muse,  University  of  Alabama,  Aerospace 
Engineering  and  Mechanics 

8:30  am 

The  Gurney  Velocity:  A  “Constant”  Affected  by  Previously  Unrecognized  Factors 
Joseph  E.  Backofen,  BRIGS  Co. 

8:50  am 

Influence  of  Post  Detonation  Burning  Process  on  Blast  Wave  Parameters  in  Air 

Meir  Mayseless,  E.  Muzychuk,  IDF,  Mil.;  M.  Mayseless,  1.  Belsky,  IMI,  Central 
Laboratory 

9:10  am 

Steerable  Fragment  Masses 

Manfred  Held,  TDW/EADS 

9:30  am 

Penetration  Performances  of  Tungsten-Copper  Shaped  Charge  Liner 

Seong  Lee,  Eun  Pyo  Kim,  Youngmoo  Kim,  Sung  Ho  Lee,  Moon-Hee  Hong, 
Joon-Woong  Noh,  Agency  for  Defense  Development 

9:50  am 

Morning  Break 

10:20  am  - 12:00  pm 

Vulnerability/Lethality/Wound  Ballistics  Poster  Session 

Chairpersons:  W.  Gooch 

Interior  Ballistics/Launch  Dynamics  Oral  Session  #2 

Chairpersons:  B.  Burns  and  R.  Cayzac 

10:20  am 

3-D  Finite-Element  Gun  Launch  Simulation  of  a  Surrogate  Excalibur  155-mm  Guided 
Artillery  Projectile  -  Modeling  Capabilities  and  its  Implications 

M.R.  Chowdhury,  A.  Frydman,  US  Army  Research  Laboratory;  J.  Cordes,  L. 
Reinhardt,  D.  Carlucci,  US  Army  ARDEC,  Analysis  and  Evaluation  Division 

10:40  am 

Method  of  Calculating  Initial  Firing  Data  of  Artillery  Laser  Terminal-Guidance  Weapon 
System 

Feipeng  Zeng,  Liren  Liu,  Faculty  of  Artillery  Command,  Nanjing  Artillery  Academy 

22nd  International  Symposium  on  Ballistics 


11:00  am 

Caseless  Ammunition  &  Advances  in  the  Characterization  of  High  Ignition 
Temperature  Propellant 

Patricia  M.  O’Reilly,  Erin  Hardmeyer,  Chad  Sensenig,  US  ArmyARDEC;  Ben 
Ashcroft,  Alliant  Techsystems;  Dave  Cleveland,  The  Johns  Hopkins  University, 
OApplied  Physics  Laboratory;  Bo  Engel,  Paul  Shipley,  AAI  Corporation 

11:20  am 

Ballistic  Launch  to  Space 

Edward  Schmidt,  M.  Bundy,  US  Army  Research  Laboratory 

11:40  am 

A  Novel  Launcher  for  Cavitating  Weapons 

Chris  J.  Weiland,  Pavlos  P  Vlachos,  Dept  of  Mechanical  Engineering 

Virginia  Tech;  Jon  J  Yagla,  Mechanical  Engineer  Engagement  Systems  Department 
Naval  Surface  Warfare  Center 

12:00  pm 

Lunch 

Warhead  Mechanisms  Oral  Session  #2 

Chairpersons:  R.  Brown  and  M.  Held 

1 :30  pm 

The  Role  of  Rayleigh  Taylor  Instability  in  Shaped  Charge  Jets  Formation  and  Stability 
Simcha  Miller,  Gershon  Kliminz,  Rafael  Ballistic  Center 

1 :50  pm 

Simulation  of  Cylinder  Expansion  Tests  Using  an  Eulerian  Multiple-Material  Approach 
Laura  K.  Donahue,  R.C.  Ripley,  Martec,  Ltd. 

2:10  pm 

Application  of  Powder  Tantalum  Material  for  Explosively  Formed  Penetrator  Warhead 
Richard  Fong,  William  Ng,  Steven  Tang,  Michael  Hespos,  US  Army  Armament 
Research,  Development  and  Engineering  Center 

2:30  pm 

Oilwell  Perforators:  Theoretical  Considerations 

Brenden  M.  Grove,  Schlumberger  Reservoir  Completions  Center 

2:50  pm 

Planar  Cutting  Jets  from  Shaped  Charges 

Geoffery  EB  Tan,  T.K.  Lam,  Y.K.  Tham,  DSO  National  Laboratories 

3:10  pm 

The  Study  on  Lethality  Simulation  Method  for  Fragment  Warhead 

Yang  Yunbin,  Qu  Ming,  Qian  Lixin,  Institute  of  Structural  Mechanics,  China 
Academy  of  Engineering  Physics 

3:30  pm 

Adjourn  for  the  Day 

4:00  pm  -  5:30  pm 

Reception 

Friday.  November  18.  2005 


7:00  am 

Continental  Breakfast  and  Registration 

8:00  am 

Adminstrative  Remarks 

Terminal  Ballistics  Oral  Session  #3 

Chairpersons:  1.  Cullis  and  D.  Nandlall 

8:10  am 

Performance  Evaluation  of  Multi  threat  Body  Armour  Systems 

B.  Anctil,  M.  Keown,  Biokinetics  and  Associates  Ltd.;  G.  Pageau,  M.  Bolduc  D. 
Bourget,  Defence  R&D  Canada  -  Valcartier 

22nd  International  Symposium  on  Ballistics 


8:30  am 

Finite  Element  Simulations  and  Experiments  to  Determine  the  Residual  Damage  of  a 
CFRP  Composite  Material  After  Ballistic  Impacts 

Koen  Herlaar,  M.  Van  der  Jagt-Deutekom,  TNO  Defence,  Security  and  Safety 

8:50  am 

The  Effect  of  Boundary  Conditions  on  the  Ballistic  Performance  of  Textile  Fabrics 
Colin  R.  Cork,  University  of  Manchester,  School  of  Materials 

9:10  am 

Terminal  Ballistic  Effects  of  Low  Density  Materials  Used  as  Confinement  Plates  for 
Explosive  Reactive  Armours 

Hanspeter  Kaufmann,  RUAG  Land  Systems;  Andre  Koch,  Armasuisse 

9:30  am 

Quantification  of  the  Effect  of  Using  the  Johnson-Cook  Damage  Model  in  Numerical 
Simulations  of  Penetration  and  Perforation 

Charles  E.  Anderson  Jr.,  T.  R.  Sharron,  Southwest  Research  Institute;  Timothy  J. 
Holmquist,  Network  Computing  Services,  Inc. 

9:50  am 

Morning  Break 

General  Oral  Session  #2 

Chairpersons:  V.  Sanchez-Galvez  and  P.  Cuniff 

10:20  am 

Comparisons  of  Internal  Ballistics  Simulations  of  the  AGARD  Gun 

Clive  R.  Woodley,  QinetiQ;  Alain  Carriere,  Patrice  Franco,  Dieter  Hensel,  Julien 
Nussbaum,  Institut  Franco-Allemand  de  Recherches  de  Saint-Louis  (ISL);  Tatjana 
Groge,  Ernst-Mach-lnstitut  (EMI);  Stefan  Kelzenberg,  Fraunhofer-lnstitut  fur 
Chemische  Technologie  (ICT),  Baptiste  Longuet,  DGA/DCE/ETBSr3 

10:40  am 

Wind  Tunnel  Verification  of  the  Performance  of  a  Smart  Material  Canard  Actuator 
Paul  Weinacht,  William  F.  Drysdale,  Travis  Bogetti,  Rod  Don,  US  Army  Research 
Laboratory;  James  T.  Arters,  Jack  R.  Vinson,  Aaron  R.  Hickman,  University  of 
Delaware;  Lamar  Auman,  US  Army  Aviation  and  Missile  RD&E  Center;  Oded 
Rabinovitch,  Technion  Israel  Institute  of  Technology 

11:00  am 

Fragmentation  of  Metal  Cylinders  Using  Thermobaric  Explosives 

M.R.  Dunning,  Defence  Research  and  Development  -  Suffield,  W.S.  Andrews, 
Department  of  Chemistry  and  Chemical  Engineering,  Royal  Military  College  of 
Canada;  K.M.  Jaansalu,  Department  of  Metallurgical  and  Materials  Engineering,  The 
University  of  Montana 

11:20  am 

A  Novel  Test  Methodology  to  Assess  the  Performance  Ballistic  Helmets 

B.  Anctil,  M.  Keown,  Biokinetics  and  Associates  Ltd.;  D.  Bourget,  G.  Pageau, 
Defence  R&D 

11:40  am 

Ballistic  Analysis  of  Bulgarian  Electroslag  Remelted  Dual  Hard  Steel  Armor  Plate 
William  Gooch,  Matthew  Burkins  and  David  Mackenzie,  US  Army  Research 
Laboratory,  Weapons  and  Materials  Research  Directorate;  Stefan  Vodenicharov, 
Institute  of  Metal  Science,  Bulgarian  Academy  of  Sciences 

12:00  pm 

Presentation  of  Awards 

The  Rosalind  and  Pei  Chi  Chou  Award  for  Young  Authors 

The  Neil  Griffiths  Memorial  Award 

The  Louis  and  Edith  Zernow  Award 

12:15  pm 

Invitation  to  the  23rd  International  Symposium  on  Ballistics,  Tarragona,  Spain,  2007 

12:25  pm 

Closing 

William  Flis,  DE  Technologies,  Inc.  &  Brian  Scott,  US  Army  Research  Laboratory 

POSTER  SESSIONS  START  HERE 

Exterior  Ballistics  Poster  Session 
1:30  pm  -  3:10  pm  Tuesday,  November  15 

1913  Fractional  Calculus  for  Design  of  Aerodynamic  Missile’s  Autopilot  and  Digital  Realization 

Bangchu  Zhang,  Chenming  Li,  Zipeng  Han ,  Zou  Yun,  Fuming  Xu,  Ballistic  Research  Laboratory  of  China,  Nanjing 

University  of  Science  &  Technology 

1915  The  Simulation  of  Rocket  Trajectory  in  Simulink 
Xin  Changfan,  Nanjing  University  of  Science  &  Technology 

1924  Establishing  a  Pitch  Damping  Test  Capability  at  CSIR  Defencetek 
Fabrizio  Dionisio,  CSIR,  Defencetek 

1944  The  Investigation  About  Using  Different  Guidance  Laws  on  Improving  Impact  Point  Deviation  of  a  Rocket 
Handong  Zhao,  Fang  Wang,  Qingshang  Liu,  Key  Laboratory  of  Instrumentation  and  Dynamic  Measurement,  North 
University  of  China 

1951  Numerical  Integration  Method  Based  on  4th  Largrange  Polynomial  of  Strap-Down  INS  System 
Guoguang  Chen,  Xiaoli  Tian,  Changfan  Xin,  Yaqi  Bao,  North  University  of  China 

1952  Research  on  Real  Time  Trajectory  Measure  Device  of  Range 
Changfan  Xin,  Guoguang  Chen,  Xiaoli  Tian,  North  University  of  China 

1953  Research  on  Attitude  Control  Strategy  of  Glide  Range  Extend  Rocket 
Xiaoli  Tian,  Guoguang  Chen,  Changfan  Xin,  North  University  of  China 

1954  Optimal  Algorithm  of  Glide  Range  Extend  Rocket’s  Trajectory 
Guoguang  Chen,  Xiaoli  Tian,  Changfan  Xin,  North  University  of  China 

1978  Practical  Propulsion  by  Directed  Energetic  Processes 

Joseph  P.  Backofen,  BRIGS  Co. 

2014  Investigating  the  Method  of  Obtaining  Ammunition  Roll  Attitude  by  Detecting  the  Geomagnetic  Vector 
Hongsong  Cao,  Guoguang  Chen,  Department  of  Mechatronics  Engineering,  North  University  of  China 

2065  External  Ballistic  Trajectory  Computations  for  Direct/Indirect  Fire  Weapon  Systems 
David  J.  Norton,  General  Dynamics  Canada 

2087  The  Influence  of  Laser  Rangefinder  Parameters  on  the  Hit  Probability  in  Direct  Tank  Fire 
Vladimir  Cech,  OPROX,  Inc.,  Jiri  Jevicky,  Department  of  Mathematics,  University  of  Defense 

2123  Flight  Dynamics  Modeling  and  Experiment  for  Composite  Concepts.  Application  to  Ribbon  Aerodynamic 
Stabilization 

Christopher  Grignon,  S.Heddadj,  Giat  Industries 

2128  Onboard  Measurements  with  Magnetic  Sensors:  Determination  of  the  Attitude  and  the  trajectory  Position 

V.  Fleck,  E.  Sommer ;  S.  Changey,  French-German  Research  Institute  (ISL);  D.  Beauvois,  Ecole  Superirure  d’Electrecite 

(Supelec) 

21 33  Aerodynamic  Characteristics  of  a  Long  Range  Spinning  Artillery  Shell  Obtained  from  3D  Magnetic  Sensors 
V.  Fleck,  E.  Sommer ;  C.  Berner,  French-German  Research  Institute  (ISL);  A.  Dupuis,  DRDC 

2143  Experimental  Testing  and  Numerical  Simulation  of  Separation  Disturbances  for  Two-Stage  Kinetic  Energy  Missiles 
Nicolas  Parise,  SNC  Technologies,  Inc.;  Richard  Lestage,  Francoise  Lesage,  Precision  Weapons  Section,  Defense 
Research  and  Development  Canada  Valcartier 

2147  Numerical  Study  on  the  Base  Drag  Characteristics  of  a  Base  Bleed  Projectile  with  a  Central  Propulsive  Jet 
Chang-Kee  Kim,  Agency  for  Defense  Development;  J.Y.  Choi,  Pusan  National  University,  Department  Aerospace 
Engineering 


2169  Solid  Fuel  Ramjet  (SFRJ)  Propulsion  for  Artillery  Projectile  Applications  -  Dynamic  Testing  Progress 
Anton  Stockenstrom ,  Dynax 

3010  Pitch  and  Bending  During  In-Flight  Extension 

W.  G.  Reinecke ,  Institute  for  Advanced  Technology;  M.  G.  Miller,  Physical  Sciences,  Inc. 

4004  Improvements  in  Aerodynamic  Design  for  KE  Less-Lethal  Projectiles 
Jamie  H.  Cuadros,  Arts  &  Engineering 


Terminal  Ballistics  Poster  Session  #1 
3:40  pm  -  5:20  pm  Tuesday,  November  15 

1001  Numerical  Simulations  of  Silicon  Carbide  Tiles  Impacted  by  Tungsten  Carbide  Spheres 

Constantine  G.  Fountzoulas,  Jerry  C.  LaSalvia,  Bryan  A.  Cheeseman,  Weapons  and  Materials  Research  Directorate; 
Michael  J.  Normandia ,  Ceradyne,  Inc. 

1007  Shock  Mitigation  for  Blast  Protection  Using  Hertzian  Tapered  Chains 

Robert  Doney,  US  Army  Research  Laboratory;  Surajit  Sen,  Department  of  Physics,  State  University  of  New  York  at 
Buffalo 

1011  A  Predictive  Model  for  the  Dwell/Penetration  Transition  Phenomenon 
Jerry  C.  LaSalvia,  US  Army  Research  Laboratory 

1012  Effect  of  Ceramic  Thickness  on  the  Dwell/Penetration  Transition  Phenomenon 
Jerry  C.  LaSalvia,  US  Army  Research  Laboratory 

1014  The  Development  of  Hybridized  Thermoplastic-Based  Structural  Materials  with  Applications  to  Ballistic  Helmets 
Shawn  Walsh,  Brian  R.  Scott,  David  M.  Spagnuoio,  AMSRD-ARL-WM-MB 

1015  Time  Resolved  Observation  of  the  Deformation  and  Surface  Strain  of  a  Textile  Fabric  Subject  to  Ballistic  Impact 
Brian  Scott,  Peter  Dehmer,  US  Army  Research  Laboratory;  Timothy  Schmidt,  Trilion  Quality  Systems 

1016  Analytic  Design  Trends  of  Fabric  Armor 

Brian  Scott,  Chian-Fong  Yen,  US  Army  Research  Laboratory 

1018  High-Speed  Photographic  Study  of  Wave  and  Fracture  Propagation  in  Fused  Silica 

Elmar  Strassburger,  Fraunhofer-lnstitut  fur  Kurzzeitdynamik,  Ernst-Mach-lnstitut  (EMI);  Parimal  Patel,  James  W. 

McCauley,  US  Army  Research  Laboratory;  Douglas  W.  Templeton,  US  Army  TARDEC 

1022  Low  Velocity  Ballistic  Properties  of  Shear  Thickening  Fluid  (STF)-Fabric  Composites 

M.  J.  Decker,  R.  G.  Egres,  N.  J.  Wagner,  University  of  Delaware,  Dept,  of  Chemical  Engineering  and  Center  for 
Composite  Materials;  E.  D.  Wetzel,  U.S.  Army  Research  Laboratory 

1023  An  Approximate  Solution  of  the  Long-Rod  Penetration  Equations 
William  Walters,  Cyril  Williams,  ARL,  Terminal  Effects  Division 

1901  Tubular  Projectile  Interaction  with  Stationary  and  Moving  Oblique  Plates 

Olof  Andersson,  Swedish  Defence  Research  Agency  (FOI),  Weapons  and  Protection  Division 

1911  A  Study  on  the  Moving  Features  of  Double-Layer  Explosive  Reactive  Armor  with  Definite  Angle  by  Numerical 
Simulation  and  Experiments 

Zhengxiang  Huang,  Xianfeng  Zhang,  Gang  Li,  School  of  Mechanical  Engineering,  Nanjing  University  of  Science  & 
Technology 

1927  Mechanics  of  Structural  Design  of  EPW  Warhead 

X.W.  Chen,  Institute  of  Structural  Mechanics,  China  Academy  of  Engineering  Physics 

1928  Armour  Qualification  Utilizing  Maximum  Likelihood  Ballistic  Limit  Calculation 
Moshe  Ravid,  Shlomo  Galperin;  Rimat  Advanced  Technologies,  Ltd. 


1934  Perforation  of  Concrete  Targets  by  an  Eroding  Tungsten-Alloy  Rod 
Stephan  Lam  pert,  Rene  Jeanquartier,  Armasuisse 

1 955  Ballistic  Properties  of  Single-Melt  Titanium-6Aluminum-4Vanadium  Alloy  Plate 

Brij  J.  Roopchand,  US  Army  Tank-Automotive  and  Armament  Command,  Armament  Research,  Development,  & 
Engineering  Center 

1987  Preliminary  Investigations  of  Potential  Light  Weight  Metallic  Armour  Applications 
Martin  van  de  Voorde,  A.M.  Dierderen,  K.  Herlaar,  TNO  Defence,  Safety  and  Security 

2001  Oblique  Warhead  Penetration  and  Perforation  of  Multi-Layered  Metallic  Targets 

Yongxiang  Dong ,  Feng  Shunshan,  Wang  Fang,  State  Key  Laboratory  of  Explosion  Science  &  Technology,  Beijing 
Institute  of  Technology 

2019  Influence  of  Projectile  Material  on  Yawed  Long  Rod  Projectiles  Penetrating  Oblique  Plates 
Ewa  Liden,  Swedish  Defence  Research  Agency  (FOI),  Weapons  and  Protection  Division 

2035  Advanced  Aliphatic  Polyurthane  Resins  for  High  Durability  and  Superior  Ballistic  Performance  Laminated  Glass 
Francisco  Folgar,  INTER  Materials,  LLC 

2037  Impact  and  Penetration  of  B4C  Ceramic,  Aluminum,  and  Berlyllium  by  Depleted  Uranium  Rods  at  2.0  KM/S 
Scott  A.  Muiiin,  James  D.  Walker,  Carl  E.  Weiss,  Southwest  Research  Institute;  Paui  O.  Lesiie,  Los  Alamos  National 
Laboratories 

2122  A  Comparison  of  Some  Analytical  and  Empirical  Models  for  Kinetic  Energy  Penetration  of  Semi-Infinite  and  Finite 
Thickness  Steel  Targets 

Nick  J.  Lynch,  J  T  Miiis,  QinetiQ 

2181  Computed  Tomography  of  High-Speed  Events 

Karsten  Michaei,  Phiiip  Heiberg,  Fraunfoher  Institute  for  High  Speed  Dynamics,  Ernst-Mach-lnstitut 

2186  Characterization  of  Behind-Armor  Debris  Particles  from  Tungsten  Penetrators 
Brad  A.  Pedersen,  S.  Bless,  Institute  for  Advanced  Technology 

Terminal  Ballistics  Poster  Session  #2 
8:10  am  -  9:50  am  Wednesday,  November  16 

2050  On  the  Critical  Thickness  of  Ceramic  to  Shatter  WC-Co  Bullet  Cores 

Paul  J.  Hazell,  Engineering  Systems  Department,  Cranfield  University,  Royal  Military  College  of  Science;  C.  J.  Roberson, 
Advanced  Defence  Materials  Limited 

2060  The  Effect  of  Spaced  Armour  on  the  Penetration  of  Shaped  Charge  Warheads 
James  D.  Shattock,  Cranfield  University 

2072  Modeling  Impact  and  Penetration  Using  a  Deterministic  and  Probabilistic  Design  Tool 

David  S.  Riha,  Jason  B.  Pleming,  Ben  H.  Tucker ;  Scott  A.  Muiiin,  James  D.  Waiker,  Carl  E.  Weiss,  Southwest 

Research  Institute;  Edward  A.  Rodriguez,  Paui  O.  Lesiie,  Los  Alamos  National  Laboratories 

2106  On  the  Ballistic  Efficiency  of  the  Three  Layered  Metallic  Targets 

Stanislav  Role,  Military  Technical  Institute  of  Protection;  Jaroslav  Buchar,  Mendel  University;  Giovanni  Cozzani,  OTO 
MELARA  S.p.A  Vojtech  Hruby,  University  of  Defence 

2107  Effect  of  the  Temperature  on  the  Ballistic  Efficiency  of  Plates  Made  From  Cast  Iron 
Stanislav  Role,  Military  Technical  Institute  of  Protection;  Jaroslav  Buchar,  Mendel  University 

2113  Displacement  Device  to  Measure  the  Acceleration  of  the  Bulge  of  RHA  Plates  Under  Anti-Tank  Mine  Blast 
Manfred  Held,  TDW/EADS;  Peter  Heeger,  WTD;  Josef  Kiermeir,  CONDAT 


2116  Defeating  Mechanisms  of  Explosive  Reactive  Armour  Sandwiches 
Manfred  Held,  TDW/LFK/EADS 

2121  Comparisons  of  Unitary  and  Jacketed  Rod  Penetration  into  Semi-Infinite  and  Oblique  Plate  Targets  at  System 
Equivalent  Velocities 

John  Stubberfield,  N  J  Lynch,  QinetiQ;  /  Wallis,  QinetiQ  Farnborough 

2126  Finite  Element  Simulations  and  Experiments  of  Ballistic  Impacts  on  High  Performance  PE  Composite  Material 
Koen  Herlaar ;  M.  Van  der  Jagt-Deutekom,  TNO  Defense,  Security  and  Safety 

2129  The  Use  of  Foam  Structures  in  Armoured  Vehicle  Protection  Against  Landmines 

David  A.  Cendon,  Vincente  Sanchez-Galvez,  Francisci  Galvez,  Alejandro  Enfedaque,  Departamento  de  Ciencia  de 
Materiales,  E.T.S.I.  Caminos,  Canales  y  Puertos,  Universidad  Politecnica  de  Madrid,  Spain 

2136  The  Numerical  Simulation  of  the  Impact  of  an  Aluminum  Cylinder  into  a  Steel  Cone 
Izak  M.  Snyman,  Defencetek  Landwards  Programme,  CSIR 

2139  Characterization  of  Al  6061-T6  using  Split  Hopkinson  Bar  Tests  and  Numerical  Simulations 

Amal  Bouamoul,  Manon  Bolduc,  DRDC  -  Valcartier 

2145  Finite  Element  Modeling  of  Light  Armoured  Vehicles  (LAV)  Welds  Heat  Affected  Zones  Sublected  to  an 
Anti-Vehicular  (AV)  Blast  Landmine  Loading:  A  Summary  of  the  Numerical  Model  and  Experiment 
Patrice  Gaudreault,  Defence  Research  &  Development  Canada;  Amal  Bouamoul,  Robert  Durocher ;  Benoit  St-Jean, 
DRDC  Valcartier 

2151  Designer  Projectiles  by  Density  Variation:  Towards  the  Nano-Projectile 
John  P.  Curtis,  QinetiQ 

2160  Numerical  Simulation  for  the  Front  Section  Effect  of  Missile  Warhead  on  the  Target  Perforation 

Ho  Soo  Kim,  Ki-Sun  Yeom ,  Seong  Shik  Kim,  Agency  for  Defense  Development  (ADD);  Larry  Sotsky,  US  Army  ARDEC 

2180  The  Penetration  Process  of  Projectiles  into  Long  Bars  in  the  Axial  Direction 
Dan  Yaziv,  G.  Gans,  Y.  Reifen,  RAFAEL 

2199  The  Electromagnetic  Launch  Trends  Utilization  for  Shaped  Charge  Jets  Penetration  Depth  Decrease 
S.V.  Demidkov,  Effective  Soft  Ltd. 

3003  Simulation  of  the  Perforation  of  Low  Mass  Long  L/D  Rods  Against  Finite  RHA  Plates 
P.  Church,  I.  Cullis,  A  Bowden,  D  Gibson,  QinetiQ,  Ltd. 

3011  Deflection  and  Fracture  of  Tungston  Rods  by  Yawed  Impact 

S.  Bless,  R.  Russell,  Instutite  for  Advanced  Technology;  K.  Tarcza,  US  Army  ARDEC;  E.  Taleff,  Department  of 
Mechanical  Engineering,  The  University  of  Texas  at  Austin;  M.  Huerta,  The  University  of  Texas  at  El  Paso 

3012  Anomalies  in  the  Strength  of  Alumina  under  Dynamic  Compression 

T.  Beno,  S.  Bless,  S. Nichols,  Instituite  for  Advanced  Technology 

3014  On  the  3D  Visualization  of  Ballistic  Damage  in  TI-6AL-4V  Applique  Armour  with  X-Ray  Computed  Tomography 
J.M.  Wells,  JMW  Associates;  W.H.  Green,  N.L.  Rupert,  US  Army  Research  Laboratory,  Weapons  &  Materials  Research 
Division;  John  M.  Winter ;  Jr.,  ORISE  Contractor  at  WMRD;  S.J.  Cimpoeru,  DSTO  Melbourne 

4012  Analytical  Models  for  Foam,  Ice  and  Ablator  Impacts  into  Space  Shuttle  Thermal  Tiles 
James  D.  Walker ;  Sidney  Chocron,  Walt  Gray,  Southwest  Research  Institute 

4013  CTH  Simulations  of  Foam  and  Ice  Impacts  into  the  Space  Shuttle  Thermal  Protection  System  Tiles 
Sidney  Chocron,  Walt  Gray,  James  D.  Walker,  Southwest  Research  Institute 

4016  Damage  Created  in  Composite  Sheet  by  Explosives  -  Effects  of  Fibre  Type,  Explosive  Mass  and  Attenuating 
Material 

M.  R.  Edwards,  R.  Unwin,  Centre  for  Materials  Science  and  Engineering,  Cranfield  University 


Warhead  Mechanisms  Poster  Session 
10:20  am  - 12:00  pm,  November  16 

1914  Experimental  Investigation  of  Equivalent  Blast  Characteristics  for  Aluminiferous  Explosive  in  Shallow  Underwater 
Wenbin  Gu,  Jianqing  Liu,  Qingli  Su,  Weiping  Zhou,  Ballistic  Research  Laboratory  of  China,  Nanjing  University  of 
Science  &  Technology 

1 935  Break-up  of  Copper  Shaped  -  Charge  Jets:  A  Combined  Experimental/Numerical/Analytical  Approach 
Jacques  Petit,  Centre  d’Etudes  de  Gramat;  V.  Jeandaude,  C.  Fressemgeas,  Laboratoire  de  Physique  et  Mecanique  des 
Materiaux,  Universite  de  Metz/CNRS 

1 949  A  Theoretical  Analysis  for  Initial  Fragment  Velocity  and  Peak  Overpressure  of  a  Blast  Fragmentation  Device 
Jin  Jianming,  Institute  of  Structural  Mechanics,  China  Academy  of  Engineering  Physics 

1963  Scaling  the  Dynamic  Response  of  Armored  Vehicle’s  Floor  Subjected  to  a  Large  Buried  Charges 

Avidov Neuberger,  IMOD,  MANTAK,  Tank  Program  Management;  S.  Petes,  IMI  Central  Laboratory  Division;  D.  Rittel, 

Technion,  Israel  Institute  of  Technology,  Faculty  of  Mechanical  Engineering 

1 984  High-Speed  Flash  X-Ray  Computed  Tomo-Cinematography 

Philip  Heiberg,  Karsten  Michael,  Fraunhofer  Institute  for  High-Speed  Dynamics,  Ernst-Mach-lnstitut 
2002  The  Influence  of  Parameters  Other  Than  Liner  Velocity  on  Shaped  Charge  Jet  Coherence 

Frederik  Mostert,  CSIR  Defencetek;  C.  J.  Terblanche,  Denel  Land  Systems  -  Western  Cape;  M.  F.  Maritz,  Department  of 
Applied  Mathematics,  University  of  Stellenbosch 

2053  Comparison  of  Vulnerability  and  Performances  of  Insensitive  Munitions  (IM)  and  Non  IM  Directed  Energy 
Warheads 

Frederic  Peugeot,  MSIAC,  NATO  HQ 

2059  Trumpet  Shaped  Liners’  Influence  on  Slug  Properties 
Eitan  Hirsch,  Consultor;  Meir  Mayseless,  IMI  Central  Laboratory 

2077  A  Study  on  the  Structure  of  Small  Caliber  EFP 

Chen  Zhigang,  Zhao  Taiyong,  Hou  Xiucheng,  Dong  Surong,  The  Research  Institute  of  Explosive  Demolition  & 

Defence  Technology,  North  University  of  China;  You  Zheng,  Department  of  Precision  Instruments  and  Mechanics, 
Tsinghua  University 

2078  A  Framework  for  the  Analyses  and  Visualization  of  X-Ray  Computed  Tomography  Image  Data  using  a  Compute 
Cluster 

Jeffrey  R.  Wheeler,  US  Army  Research  Laboratory;  William  H.  Green,  US  Army  Research  Laboratory,  Weapons  & 
Materials  Research  Division;  Michael  Schuresko,  Baskin  School  of  Engineering,  University  of  California;  Michael  Patrick 
Lowery,  Computational  Mathematics  Department,  University  of  California 

2137  An  Artillery  Shell  for  Anti-Bunker  Applications  (155  ABS) 

Remi  Boulanger,  Giat  Industries;  Anders  Vangen  Jordet,  Dagfinn  Hoff,  NAMMO  Raufoss 

2141  Development  of  a  TMRP-6  Surrogate  Mine 

Wes  Baillargeon,  A.  Sirois  and  G.  McIntosh;  DRDC  Valcartier 

2148  Use  of  Foamed  TNT  Mixtures  as  a  Dispersion  Charge  of  Submunitions 

Jun  Sik  Hwang,  S.-K.  Kwon,  C.-K.  Kim,  S.-W.  Kwon,  S.-S.  Kim,  S.-H.  Moon,  Explosive  Trains  and  Gun  Propellant 
Team,  Agency  for  Defense  Development 

2152  Wall  Breaching  Tandem  Warhead 

Andreas  Helte,  Torgny  Carlsson,  Hakan  Hansson,  Svante  Karlsson,  Jonas  Lundgren,  Lars  Westerling,  Hakan 
Ornhed,  Swedish  Defence  Research  Agency,  FOI,  Weapons  and  Protection  Division 

2156  An  Analytical  Penetration  Model  for  Jets  with  Varying  Mass  Density  Profiles 

Milton  F.  Maritz,  Stellenbosch  University;  Klaus  D.  Werneyer,  TTP  Products;  Frederik  J.  Mostert,  Defencetek 


2168  Peripheral  Initiation  Technology  Development 

Arthurs.  Daniels,  Ernest  L.  Baker,  William  J.  Poulos,  Vladimir  M.  Gold,  B.  Fuchs,  US  Army  ARDEC 

2175  Shaped-Charge  Jet  Stability  Calculations:  The  Role  of  Initial  and  Boundary  Conditions 
James  S.  Stolken,  S.  Christian  Simonson,  Mukul  Kumar,  Lawrence  Livermore  National  Laboratory 

3007  Enhanced  Focused  Fragmentation  Warhead  Study 

Richard  Fong,  William  Ng,  Peter  Rottinger,  Steve  Tang,  US  Army  Armaments  Research,  Development  and  Engineering 
Center 

4015  Microstructure  and  Properties  of  the  Explosively  Formed  Petals  in  Aluminium  Alloys 
M.  R.  Edwards,  J.  M.  Cassar,  Centre  for  Materials  Science  and  Engineering,  Cranfield  University 

Interior  Ballistics/Launch  Dynamics  Poster  Session 
8:10  am  -  9:50  am  November  17 

1017  In-Bore  Mechanics  Analysis  of  the  M855  Projectile 
Joseph  T.  South,  James  F.  Newill,  US  Army  Research  Laboratory 

1020  Dynamic  Strain  Measured  in  a  105-mm  Composite  Gun  Barrel  -  A  Fiction  or  Reality 
Jerome  T.  Tzeng,  US  Army  Research  Laboratory 

1916  A  Vector  Way  for  Calculating  Propellant’s  Combustion  Performance 

Wei  Zhifang,  Department  of  Mechanical  and  Electronic  Engineering,  North  University  of  China 

1931  Interior  Ballistics  Code  Applied  to  ETC  Concept:  Computations  and  Validations 
Gilles  Legeret,  Dominique  Boisson,  Giat  Industries 

2004  The  FHIBS  Internal  Ballistics  Code 

Clive  R.  Woodley,  Steve  Billett,  QinetiQ;  Caroline  Lowe,  Department  of  Applied  Mathematics  and  Theoretical  Physics, 
Centre  for  Mathematical  Studies,  University  of  Cambridge;  William  Speares,  The  Cylinders;  Eleuterio  Toro  Laboratory  of 
Applied  Mathematics,  Faculty  of  Engineering,  University  of  Trento 

2005  Modelling  the  Ignition  of  40mm  Gun  Charges 
Clive  R.  Woodley,  QinetiQ 

2020  Thermo-Mechanical  Erosion  Study  of  the  120mm  Chromium  Coated  Gun  Barrel:  Computation  and  Validation  of  the 
Heat  Exchange  Boundary  Condition 

Dominique  Boisson,  Gilles  Legeret,  Roxan  Cayzac,  Giat  Industries 

2031  MOBIDIC-NG:  A  1 D/2D  CFD  Code  Suitable  for  Interior  Ballistics  and  Vulnerability  Modelling 
Baptiste  Longuet,  Pascal  Millet,  Eric  Taiana,  ETBS;  Patrick  Della,  Pieta  Christiane  Reynaud,  SNPE  Materiaux 
Energetiques  CRB;  Patrice  Franco,  Alain  Carrere,  Institut  Franco-Allemand  de  Recherches  de  Saint-Louis  (ISL);  Gilles 
Legeret,  Dominique  Boisson,  Giat  Industries;  Alexandre  Papy,  ERM  ABAL  30 

2070  Barrell  Life  Results  of  the  5.56  mm  XC77A1  Cartridge 
Etienne  Munger,  SNC  Technologies,  Inc. 

2117  Further  Investigation  of  the  Effect  Known  as  Electrothermal  Pyrolysis 
Steve  R.  Fuller,  M.J.  Taylor,  QinetiQ 

2150  Determination  of  Force  and  Temperature  Impact  on  Missile’s  Fuel  Charge  in  Process  of  Ignition 
Dmitriy  Orlov,  GDT  Software  Group 

2161  Unsteady  Intermediate  Ballistics:  2D  and  3D  CFD  Modelling,  Application  to  Sabot  Separation 
Roxan  Cayzac,  Eric  Carette,  Giat  Industries,  Division  Munitions;  Thierry  Alziary  de  Roquefort,  Universite  de  Poitiers, 
Laboratoire  d’Etudes  Aerodynamiques;  Philippe  Bidorini,  Emmanuel  Bret,  Pascal  Delusier,  Serge  Secco,  DGA/ETBS, 
Direction  de  I’Expertise  Technique 


2166  Large  Caliber  Firing  with  Electro  Thermal-Chemical  Ignition  (ETI) 

Jonathan  D.  Shin,  John  J.  O’Reilly,  David  T.  Keyser,  US  Army  Research,  Development  and  Engineering  Center  - 
TACOM;  Jahn  Dyvik,  United  Defense  L.R 

3006  Rail  Gun  Test  Projectile  for  Improved  Developmental  Testing  of  Precision  Munition  Electronics 

T.  Myers,  D.  Cariucci,  J.A.  Cordes,  US  Army  ARDEC,  Analysis  and  Evaluation  Division,  Fuze  and  Precision  Munitions 

Technology  Directorate 

4010  Improved  Mortar  Barrel  Thermal  Model 

M.  Pocock,  C.  Guyott,  Frazer-Nash  Consultancy  Ltd;  P.  Locking,  BAE  Systems,  Land  Systems 

Vulnerability/Lethality/Wound  Ballistics  Poster  Session 
10:20  am  - 12:00  pm  November  17 

1855  On  Incorporating  XCT  into  Predictive  Ballistic  Impact  Damage  Modeling 
Joseph  M.  Wells,  JMW  Associates 

1878  New  Soft-Target  Failure  Criteria  for  System-Analytical  Considerations 

Markus  J.  Estermann,  RUAG  Defence,  Warhead  Division;  Beat  P.  Kneubuhl,  Aramasuisse 

1941  Protecting  Vehicles  from  Landmine  Blasts 

Sheri  L.  Hlady,  Denis  Bergeron,  Defence  R&D  Canada  -  Suffield;  Rene  Gonzalez,  US  Army,  PM  Light  Tactical  Vehicles 

1957  Office  of  Naval  Research  Limb  Protection  Program 
Graham  K.  Hubler,  NRL 

1980  Survivability  and  Lethality  Assessment  Software  Based  on  Virtual  Mode  Technology 
Lu  Yonggang,  Qian  Lixin,  Yang  Yubin,  Liu  Tong,  Institute  of  Structural  Mechanics,  CAEP 

1981  “TBM-Xpert”  -  A  New  Endgame  Code:  Features  and  Validation 

Werner  Arnold,  EADS-TDW  Gesellschaft  fur  verteidigungstechnische;  E.  Rottenkolber,  NUMERICS  GmbH 

1989  The  Use  of  Ballistic  Knowledge  in  Ammunition  Safety  Cases 
Martin  van  de  Voorde,  TNO  Defence,  Safety  and  Security 

2011  A  Note  on  the  Roecker-Ricchiazzi  Model  of  Penetrator  Trajectory  Instability 
William  J.  Flis,  DE  Technologies,  Inc. 

2022  Numerical  Calculation  and  Simulation  of  Missile  Jet-Airplane  Interaction 
Feipeng  Zeng,  Faculty  of  Artillery  Command,  Nanjing  Artillery  Academy 

2111  Need  for  Enhanced  Protection  Against  Blast  Threats  for  Soldiers  Exposed  to  Roadside  Improvised  Explosive 
Devices  (lEDs) 

Frangois-Xavier  Jette,  Jean-Philippe  Dionne,  Aris  Makris,  Med-Eng  Systems,  Inc.;  Karl  Masters,  PEO  Soldier; 
Christine  Perritt,  PM  Soldier  Equipment 

2119  WitnessMan:  The  Software  Tool  to  Design,  Analyse  and  Assess  a  Witness  Pack  with  Respect  to  Military  and 
Medical  Effects  on  an  (Un)protected  (Dis)mounted  Soldier. 

Theo  L.A.  Verhagen,  R.  Kemper,  H.  Huisjes,  S.G.  Knijnenburg,  A.  Pronk,  M.H.  van  Klink,  TNO  Defence,  Security  and 
Safety 

2154  Injury  Risks  Resulting  from  Deminer  Position 

Frangois-Xavier  Jette,  Jean-Philippe  Dionne,  Ismail  El  Maach,  Aris  Makris,  Matt  Ceh,  Med-Eng  Systems,  Inc.;  Denis 
Bergeron,  Defence  R&D  Canada  Suffield 

2163  RPG  Mitigation  for  Military  Vehicles 

Karl  Pfister,  Dipl.  Ing  (FH)  Armatec  Survivability  Corporation 


2164  Protection  Against  Closely-Spaced  Impacts  by  Small  Arms  Bullets 

Michael  J.  Iremonger,  Cranfield  University,  Royal  Military  College  of  Science;  Abdullah  Alsalmi 

301 3  Vulnerability  Evaluations  of  30mm  Airburst  Ammunition 

Quoc  Bao  Diep,  Eimund  Smedstad,  Nammo  RaufossAS;  Nick  Rogers,  System  Design  Evaluation  (SDE) 

3018  Challenges  and  a  Solution  in  Determining  Land  Mine  or  IED  Neutralization  Effectiveness 
Robert  Colbert ,  Mark  Majerus,  William  dark,  DE  Technologies,  Inc. 

4011  Numerical  and  Experimental  Analysis  of  the  Detonation  of  Sand-Buried  Mines 
N.  Heider,A.  Kiomfass,  Fraunhofer-lnstitut  fur  Kurzzeitdynamik,  Ernst-Mach-lnstitut 


The  symposium  registration  fees  are: 


Regular  Late/Onsite 
after  10/28/05 

$950  (US)  $1045  (US) 


IBC  Committee  Appreciation  Dinner 
Guest  at  Both  Receptions 
Guest  at  One  Reception 


$75  (US) 
$75  (US) 
$50  (US) 


The  symposium  registration  fee  includes  attendance  at  all  sessions,  bound  symposium  proceedings  with  CD, 
continental  breakfasts,  coffee  breaks,  lunches,  receptions,  and  administrative  costs.  The  registration  fee  will 
also  include  a  compact  disc  (CD)  which  contains  a  cumulative  database  of  titles,  authors  and  abstracts  of  all  of 
the  22  Ballistics  Symposia. 

To  register  online  for  this  conference  visit:  http://reaister.ndia.org/interview/reaister.ndia7~Brochure-6210.  You 
can  also  visit  the  NDIA  website  at  www.ndia.org  and  select  “Schedule  of  Events”.  Then  select  2005  November 
and  scroll  down  to  the  22nd  International  Symposium  on  Ballistics.  Once  there,  select  the  blue  “Register”  link  in 
the  lower  left  hand  corner  of  your  screen.  Review  your  information  and  then  select  “submit”  one  time  only  and 
then  select  “confirm”.  On-line  registration  will  close  at  5:00  pm  EST  on  October  28,  2005.  You  must  register 
on-site  after  this  date. 

-or- 

You  may  fax  the  completed  registration  form  contained  in  this  brochure  to  (703)  522-1885. 

-or- 

You  may  mail  the  completed  registration  form  contained  in  this  brochure  to:  Event  #  6210,  National  Defense 
Industrial  Association,  2111  Wilson  Boulevard,  Suite  400,  Arlington,  VA  22201 -3061. 

Payment  must  be  made  at  the  time  of  registration.  Registrations  will  not  be  taken  over  the  phone. 

Cancellations  and  Refund  Policy 

Registrants  who  cannot  attend  the  22nd  International  Symposium  on  Ballistics  must  provide  written  notification 
via  email  to  bbommelje@ndia.org  or  fax  to  (703)  522-1885  on  or  before  September  16,  2005  to  avoid  a 
cancellation  fee. 

Cancellations  received  between  September  16,  2005  and  October  28,  2005  will  receive  a  refund  minus  a  $75 
cancellation  fee.  No  refunds  will  be  given  to  cancellations  recieved  after  October  28,  2005  however, 

SUBSTITUTIONS  ARE  WELCOME  IN  LIEU  OF  CANCELLATIONS. 

You  must  have  a  government  picture  identification  (drivers  license,  passport,  military  ID,  etc.)  to  receive  a 
symposium  badge.  Badges  must  be  worn  at  all  times  during  the  symposium. 

Special  Needs 

NDIA  supports  the  Americans  with  Disabilities  Act  of  1990.  Attendees  with  special  needs  should  call 
(703)  522-1820  prior  to  October  3,  2005. 

Hotel  Accommodations 

A  limited  block  of  rooms  have  been  reserved  at  the  Fairmont  Waterfront  Hotel.  The  industry  room  rate  is  $219 
Canadian  (approximately  $180  US).  The  government  symposium  room  rate  is  approximately  $114  Canadian 
($94  US).  Please  call  (604)  691-1991  to  make  reservations. 

In  order  to  ensure  the  discounted  NDIA  rate,  please  make  reservations  early  and  ask  for  the  NDIA  room  block. 
Rooms  will  not  be  held  after  Tuesday,  October  11,  2005  and  may  sell  out  before  than.  Rates  are  also  subject 
to  increase  after  this  date. 


*The  government  room  rate  applies  only  to  active  duty  military  and  civilian  government  employees.  It  is  not 
available  to  government  contractors ,  retired  military  or  retired  civilian  government  employees.  ID  cards  and/or 
travel  orders  will  be  required  at  check-in  to  verify  rate  eligibility. 


GENERAL  INFORMATION 


Symposium  Attire 

Appropriate  dress  for  this  symposium  is  business  for  civilians  (coat  and  tie)  and  class  A  uniform  for  military. 
Inquiries 

For  more  information  regarding  the  symposium  contact  Britt  Bommelje,  Meeting  Planner  at  (703)  247-2587  or 
bbommelje@ndia.org. 

Promotional  Partnerships 

Increase  your  company  or  organization  exposure  at  this  premier  event  by  becoming  a  Promotional  Partner.  A 
Promotional  Partnership  ($5,000)  will  add  your  company  name  to  the  back  cover  of  the  on-site  brochure  as  well 
as  main  platform  recognition  throughout  the  conference,  signage  at  all  events  including  the  opening  reception 
and  a  350-word  organization  description  in  the  conference  agenda.  For  more  information,  please  contact  Sam 
Campagna  at  703-247-2544  or  scampagna@ndia.org. 

www.defensejobs.com 

The  Defense  Industry’s  leading  employment  website;  find  a  job,  post  a  job  listing,  post  a  resume,  and  search 
resumes.  For  more  information  please  contact  info@defensejobs.com  or  (703)  247-9461 .  Please  visit 
www.defensejobs.com 

IBC  Committee  Appreciation  Dinner 

The  IBC  Committee  Appreciation  Dinner  will  be  held  on  Friday,  November  18,  2005  in  Vancouver.  This  dinner 
is  open  to  IBC  Committee  members  and  their  guests  only.  If  you  and  your  guest  would  like  to  attend,  please 
make  a  note  of  it  on  the  registration  form.  There  is  a  $75  charge  per  person  to  attend. 


“The  Department  of  Defense  finds  this  event  meets  the  minimum  regulatory  standards  for  attendance 
by  DoD  employees.  This  finding  does  not  constitute  a  blanket  approval  or  endorsement  for  attendance. 
Individual  DoD  component  commands  or  organizations  are  responsible  for  approving  attendance  of  its 
DoD  employees  based  on  mission  requirements  and  DoD  regulations.  ” 
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22nd  International  Symposium  on  Ballistics 

Vancouver  Convention  Center  ~  Vancouver,  BC,  CANADA 

November  14-18,  2005 
Event  #6210 

Organization  Information 


Organization  Name  (as  it  should  appear  on  booth  sign  --  limited  to  40  characters  and  spaces) 

Point  of  Contact  (for  fees  and  exhibitor  service  kit) 

Title 

Street  Address 

City 

State 

Zip  Code 

Telephone 

Fax 

E-mail 

Reserve  your  booth  and  register  your  exhibit  staff  at  http://exhibits.ndia.org 


Exhibit  Space  Information: 
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Is  it  possible  to  find  light-weight  and  flexible 
body  armours  capable  of  defeating  multiple 


Why  should  wo  care? 


Common  threats  to 
military  and  law 
enforcement  personnel: 

•  Small  arm  projectiles 

•  Fragmentation  from 
explosive  devices  (e.g. 
lEDs) 

•  Knife,  spike,  flechette 

Can  they  be  protected 
against  all  these  threats 
with  a  single  body 
armour  system? 


Our  Strategy 


Established  multiple  performance  requirements 

Select  appropriate  test  methods 

Contact  mfgs  to  propose  solutions  and  provide 
test  samples 

Evaluate  armour  materials 


Performance  Requirements  & 
T  est  Methods 
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Drop  mass  method 
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1  Armour  Sample 

Description 
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PSDB  Level  KR1 
NIJ  Stab  Level  2 
NIJ  Stab  Level  2 
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Woven  aramid  and  laminated  polyethylene,  26  layers  NIJ  Ballistic  Class  II 

total  NIJ  Stab  Level  2 

Woven  PBO  and  aramid,  laminated  polyethylene,  50  NIJ  Ballistic  Class  IIIA 

layers  total  NIJ  Spike  Level  2 


10** 
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Custom 

2.7 

Coated  woven  aramid,  woven  aramid,  laminated 
polyethylene,  41  layers  total 

Custom 

11.9 

2** 

Coated  woven  aramid,  woven  PBO,  laminated 
polyethylene,  51  layers  total 

Custom 

12.3 

13** 

Woven  PBO,  unidirectional  aramid,  38  layers  total 

Custom 

6.3 
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Steel  sheets  and  woven  fabric,  15  layers  total 


Coated  woven  aramid,  30  layers 


Coated  woven  aramid,  30  layers 


Coated  woven  polyethylene,  27  layers 


2  types  of  woven  aramid,  32  layers  total 


2  types  of  woven  aramid,  31  layers  total 


Multi-layers  of  dense  woven  aramid,  18  layers 
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Areal  Density  (kg/m 


Stab  Resistance  (Level  2) 


20.0 
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cgk  Armour  Sample 


Sample  1  tested  at  Level  1 
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Flechette  Resistance 


Armour  Sample 
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Fragmentation  /  Bullet  Resistance 


Ballistic  Limit 


1000 


Sphere 


Armour  Sample 


Fragmentation  /  Bullet  Resistance 


Ballistic  Limit  Back  Face  Signature 
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Armour  Sample  Armour  Sample 
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Performance  Comparison 


140 

120 

100 

80 

60 

40 

20 

0 


9 


4 


5 


3  2  10 

Armour  Sample 
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■  Sphere 
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Vso/AD  -  FSP  &  Sphere  (m7kg.s) 


Vso/AD  -  Flechette  (m  /kg.s) 


Perforamnce  Rating 


Perforamnce  R 


Conclusions 


Current  technologies  can  provide  minimum  protection 
against  stab,  flechette,  and  ballistic  threats. 

a.d.  >  9.9  kg/m2  for  3  stab  threats  (level  2) 

Preliminary  assessment  of  flechette  resistance  possible 
with  drop  mass  method 

Performance  optimization  possible  through  #  of  layers, 
sequence,  material  combination 

Desired  requirements  can  not  be  achieved  for  a.d.  <10 
kg/m2  (9mm  HP  Bofors  and  flechette) 


Way  ahead 


Test  all  samples  for  stab  Level  1 

Refine  drop  mass  method:  reduce  weight  to  increase 
velocity,  modify  flechette  simulator 

Complete  ballistic  test  evaluation  for  samples  No.  1 1 
and  12 

Evaluate  semi-rigid  solutions  (metallic  or  ceramic  tiles) 
for  high  performance  bullets 

Optimize  performance  with  minimum  aerial  density 

Consider  zones  with  different  protection  levels 
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The  Problem 


New  lightweight 
composite  helmets  have 
increased  protection 
against  penetration 

Result  in  large  back- 
face  deformation 

Increased  risk  of 
serious  skull/brain 
injuries 

No  widely  accepted 
evaluation  procedure 


Our  Strategy 


rrentrco 


Implement  an  impact  force  me 
headform 


Based  on  injury  model  develo 
(2003)  using  localized  skull  pr 

Develop  test  procedure 


Conduct  experimental  trials  wit 
helmet  models 


J 

Define  injury  function 


Measurement  System  Requirements 


Dynamic  Loads 

•  Force  <  5,000  N 

•  Duration  <2  ms 

PVDF  gauge 


Load  cell 


■  Evaluation  under 
ballistic  loading 
conditions 


Instrumentation  Selection 


Peak  Pressure  -  output  (KPa) 


Force  applied  vs.  measured 
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Peak  Force  -  input  (N) 


■  PVDF-Ktech  ♦  PVDF-Dynasen  ▲  Loadcell-Kistler 


Pressure  sensitive  film  to 
measure  loading  area 


Impact  Force  Measurement  Headform 


Aluminum  Cover 


Shaped  Impact 
P4XJ&  .. 


34  mm 


48  mm 


Aluminum  Cover 


Shaped  Impact 
P4XJ&  .. 


LwA  C  c  lb 


Helmet  Performance  Evaluation 


3  Combat  Helmet  Models 
Similar  Ballistic  Limit  (V50) 

350  m/s  <  v  <  450  m/s 


Helmet  A  Helmet  B  Helmet  C 


Headform  Response 


Peak  Force  of  Individual  Load  Cells 


Peak  Impact  Force  (N) 


Individual  Peak  Forces 


Peak  Impact  Force  (N) 


Individual  Peak  Forces 


orce  (N)  Force  (N) 


Impact  Force  Measurements 
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Centre  Kl9nt 


Time  (ms) 


Helmet  A 
405  m/s 


Helmet  C 
392  m/s 


ulse  (N.ms)  Impulse  (N.ms) 


Skull  Fracture  Injury  Function 


Impact  velocity  (m/s) 


Helmet  B  /  PMHS 
(Bass  et  al.  2003) 
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Force  vs.  velocity 


Injury  risk 
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2000  4000  6000 

Peak  Force  (N) 


valid  only  for  concentrated  load 


Load  Distribution 


Average  Peak  Pressure 


Conclusions 


Miniature  load  cell  suitable  to  measure  helmet 
backface  loading 

Instrumented  headform  was  able  to  quantify  the 
performance  of  ballistic  helmets 

Can  be  used  to  predict  the  risk  of  skull  fracture 


Limitations 


Injury  function  valid  only  for  concentrated  load 

Contact  area  >  sensing  area 

Peak  force  must  be  within  sensing  area 

Does  not  address  distributed  forces  (rigid 
helmets) 


Way  ahead 


Additional  impact  locations  (e.g.  front,  rear) 

Consider  other  measurements  (head 
acceleration) 

Review  data  analysis  procedure  (peak,  sum 
impulse) 

Calibration  procedure 

Laboratory  re-enactments  of  injurious  cases 
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JXgj/7  Research  Goal 


To  assess  the  ability  of  TBXs  to  generate  fragments. 
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fhV  Basic  TBX  Concept 


Metal 

particles 


Explosive 

dispersal 


Fuel-rich 

detonation 


Shock 


Reflected  shock 
enhances  mixing 


Confined 

mixing 


*  * 


Experimental  Charge 


a 


vN 


Detonator 


Booster 


TBX 


6  mm  casing 


140  mm 


PMR-2A 


Trial  Variables 


Cs 


*4 


102  mm 


a 


,\ 


6  mm  casing 


Explosives 

-  TBX  1,  2,  3 

-  C4  (baseline) 

Wall  thickness 

-  3.8  mm,  6  mm, 

8  mm,  and  9.5  mm 

Casing  material 

-  1026  steel 

-  Ductile  cast  iron  (DCI) 

-  Grey  cast  iron  (GCI) 


rs 
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TBX  formulations 


TBX  1 


Monopropellant  and 
magnesium  particles 


TBX  2 


TBX  3 


Nitromethane  and 

aluminium  (60/40  %wt) 


Monopropellant,  aluminium 
and  RDX 


Witness  Packs 
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Secondary  Combustion 


First  fragment  impacts 


f\J  Witness  Pack  Analysis 

L^ 


Software  converts  hole  size  and 
depth  of  penetration  to  mass  and 
velocity 


Mass  distribution  results  were 
compared  to  literature  values 

The  average  velocities  were 
compared  to  other  methods. 
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C4-filled,  6  mm,  1026  steel  casing 
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Average  Fragment  Mass 


• 

C4  (1026  steel) 

o 

TBX  1  (1026  steel) 

▼ 

TBX  2  (1026  steel) 

V 

TBX  3  (1026  steel) 

■ 

TNT  (PMR-2A) 

□ 

TBX  1  (Ductile  cast  iron) 

♦ 

TBX  1  (Grey  cast  iron) 

O 


O 

▼ 


o 


6  7 

Wall  thickness  (mm) 


T~ 

8 


9 


10 


Flash  X-ray  Site  (1/2) 


Rim  cassette 


Metal  basin 


X-ray  source  1 


X-ray  source  2 


Charge 


Plywood  charge  box 


Flash  X-ray  Site  (2/2) 


TBX  1,  6  mm,  1026  steel,  100|Js 


C4,  6  mm,  1026  steel,  90jis 


X-ray  Images 


TBX  1,  3.8  mm,  1026  steel,  IOOjis 


TBX  1,  6  mm,  ductile  cast  iron,  IOOjis 


Predicting  Fragment  Velocity 


Gumey  equation 

-  Simple  and 
long-standing 

-  Geometry-specific 


M  =  mass  of  casing 
C  =  mass  of  explosive 

4lE  =  Gurney  constant 


SplitX 

-  Gurney-based  computer  code  that  considers: 
End  confinement 

Shock  wave  interaction  and  propagation 
User-controlled  geometry 


Fragment  Velocity  (m/s) 
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Explosive 

Wall 

thickness 

Casing 

Gurney’ 
(+/-  12%) 

SplitX’ 

Witness'3 

packs 

X-ray  a 
images 
(+/-  50) 

C4 

6  mm 

1026 

1460 

1360 

610 

1400 

C4 

8  mm 

1026 

1260 

1150 

590 

— 

TBX  1 

3.8  mm 

1026 

1010 

990 

450 

1000 

TBX  1 

6  mm 

1026 

700 

770 

430 

800 

TBX  1 

6  mm 

DCI 

740 

820 

— 

800 

TBX  1 

6  mm 

GCI 

640 

810 

— 

750 

TBX  2 

6  mm 

1026 

990 

930 

460 

— 

TBX  3 

6  mm 

1026 

1030 

970 

490 

1000 

a  -  Maximum  fragment  velocity 
b  -  Average  fragment  velocity 


Conclusions 
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do 


•  TBXs  are  capable  of  fragmenting  metal  casings: 


-  Fragment  mass  distributions  were  consistent  with 
literature  values;  and 

-  Fragment  velocities  were  well  predicted  using  means 
that  assume  an  instantaneous  release  of  detonation 
energy. 


Future  Work 


Determine  how  the  casing  thickness  and  material  disrupt  the 
TBX  shock  wave. 


-  Sacrificing  the  “thermobaric  effect”  to  have  fragments 


Determine  how  the  fragmentation  trends  of  the  base 
explosives  are  altered  by  the  additives. 

-  Run  trials  with  pure  NM,  and  with  silicon  as  an  additive 

Determine  why  the  casing  material  appears  to  have  little 
effect  on  fragment  velocity. 

—  High  strain  rate  failure  and  gas  dynamics  problem 
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Fragment  Energy  (J) 


Explosive 

Casing 

thickness 

Material 

Average 
mass  (g) 

Max  velocity 
(m/s) 

Energy  (J) 

Number  of 
fragments 

C4 

3.8  mm 

1026 

2.3 

1750 

3500 

180 

C4 

6  mm 

1026 

2.8 

1400 

2700 

240 

C4 

8  mm 

1026 

3.7 

1200 

2700 

260 

TBX  1 

3.8  mm 

1026 

3.0 

1000 

1500 

140 

TBX  1 

6  mm 

1026 

6.7 

800 

2100 

100 

TBX  1 

8  mm 

1026 

9.9 

650 

2100 

100 

TBX  1 

9.5  mm 

1026 

11.0 

550 

1700 

110 

TBX  1 

8  mm 

DCI 

3.8 

800 

1200 

220 

TBX  1 

8  mm 

GCI 

1.3 

750 

370 

670 

TBX  2 

3.8  mm 

1026 

1.8 

1200 

1300 

230 

TBX  2 

6  mm 

1026 

5.6 

970 

2600 

120 

TBX  2 

8  mm 

1026 

7.3 

850 

2600 

130 

TBX  3 

6  mm 

1026 

4.7 

1000 

2400 

150 

TNT 
(PMR  2A) 

9.5  mm 

GCI 

2.1 

570 

340 

740 

Fragment  Mass  Distribution 

Mott  Approach 


N(m)  =  number  of  fragments 
heavier  than  mass  m 
Mo  =  casing  mass 


MK=Bt6d 3 


1  + 


t 


\ 


d 


J 


t  =  casing  thickness 
d  =  casing  interior  diameter 
B  =  Mott  coefficient 


Fragment  Mass  Distribution 
C4  in  6  mm  thick,  1026  Steel  Casings 


Velocity  (m/s) 


Fragment  Velocity  vs.  Time 

TBX  1  in  1026  steel  casings 
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The  Gurney  Velocity: 

A  “Constant”  Affected  by 
Previously  Unrecognized  Factors 
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The  Gurney  Velocity  -  Key  to  the  Equations 

Vf  =  (2Eg) 1/2  [(M/C)  +  1/2]  ~1/2  Cylinder 

=  Vg  (Pcyl  /  Pex)  - 1'2  [(tcyl  /  Rex) 2  +  2  (tcyl  /  Rex)  +  0.5  (pex  /  pcyl)]  - 1/2 

Vf  =  (2Eg) 1/2  [(M/C)  +  1/3]  -1/2  Symmetric  Sandwich 
Vf  =  (2Eg) 1/2  [(M/C)  +  3/5]  ~V2  Spherical  Shell 


Values  for  the  Gurney  Velocity  (2Eg) 1/2  derived  from 
experiments  using  different  cylinder  materials 


Steel  (m/s) 

Copper  (m/s) 

Comp.  A-3  (RDX) 

2416 

2630 

Cyclotol  (75/25  cast) 

2320 

2790 

Comp.  B 

2310 

2700 

TNT  (cast) 

2040 

2370 

Tetryl 

2209 

2500 

Gurney  Velocity  /  Detonation  Rate  relationships: 


Vg/D  =  0.337  (P.W.  Cooper) 

Vg  /  D  =  (0.605  /  [F  —  1])  (J.  Roth  per  J.E.  Kennedy) 

where  r  =  the  adiabatic  exponent  for  the  gaseous  products 

Vg/D  =  (0.60  ([)  _  1/2  +  0.648  p0  1/2)  /  (1 .01  +1.313  p0) 

where  <|>  =  N  M  1/2  Q  1/2 ;  N  =  moles  of  gaseous  detonation  products 
M  =  average  weight  of  gases,  and  Q  =  chemical  energy  of  detonation 

(Hardesty  &  Kennedy  /  Kamlet  &  Hurwitz) 


Copper  Cylinders 


Vg/D 

Exp.  (Licht) 

Cooper 

Roth 

HK/KH 

TNT 

0.346 

0.346 

0.350 

0.351 

Comp  B 

0.345 

0.343 

0.355 

0.385 

Octol 

0.335 

0.330 

0.331 

0.328 

LX- 14 

0.326 

0.348 

PETN 

0.359 

0.355 

0.369 

0.331 

Gurney  Velocity  (m/s) 


Comparison  of  Gurney  Velocities  on  the  Basis  of 
Cylinder  Wall  Areal  Density  from 
2-inch  Diameter  Cylinder  Tests  Using  Comp.  B  Explosive 


Gurney  Velocity  /  Detonation  Rate  (Vg  /  D) 


Ratio  of  the  Gurney  Velocity  to  an  Explosive’s  Detonation 
Rate  versus  a  Radial-Projection  Areal  Density  Ratio 


0.00 


■ 

Comp.B  in  Aluminum 

♦ 

Comp.B  in  Steel 

A 

Comp.B  in  Copper 

X 

RDX  in  Steel 

• 

RDX  in  Copper 

□ 

Octol  75/25  in  Copper 

+ 

LX- 14  in  Copper 

o 

PBX9404  in  Copper 

▲ 

LX-04  in  Copper 

X 

TATB  in  Copper 

• 

LX- 17  in  Copper 

O 

PETN  in  Copper 

-Trendline 

1.000 


Vg  /  D  =  0.302  ( toy,  Pcyl  /  R*x  pex  )  (~ 5 ' 30 1 


mooo 
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Ratio  of  Cylinder  Wall  Mass  to  Explosive  Mass  in  the  Radial  Direction  ( tcyi  pcyi  /  Rex  pex ) 


Gurney  Velocity  /  Detonation  Rate  (Vg  /  D) 


Normalized  Gurney  Velocity  Data  from  Measurements  of 
Cylinder  Wall  Velocity  at  Fracture  Time 


0.10 


1.00 

Normalized  Cylinder  Wall  Areal  Density  ( tcyi  pcyi  /  Rex  pex ) 


10.00 


BRIGS  Two-Step  Detonation-Driven  Propulsion  Model 


Initial  motion: 

•  imparted  by  a  brisant  shock-dominated  process  that  depends 
upon  intimate  contact  of  explosive  with  propelled  material 

•  envisioned  as  caused  by  higher-pressure  region  of  detonation 
front  (envision  the  von  Neumann  spike  or  reaction  zone  region 
as  a  finite  thickness  of  solid  material  squeezed  at  high  pressure) 

Gas-push  (gas-dynamic)  propulsion: 

•  envisioned  similar  to  that  assumed  by  Gurney  modeling 

(gaseous  product  expansion  from  a  homogeneous  “all  burnt” 
condition  while  pushing  confining  boundaries  to  a  final 
“steady-state”  velocity  as  the  pressure  drops) 


V  initial  /  V final 


Vjnitiai  /  Vfina,  Data  Plotted  for  Cylinders  and  Plates 
of  Various  Inert  Materials  and  Explosives 


♦ 

Copper 

■ 

Iron 

▲ 

Magnesium 

X 

Tantalum 

• 

Aluminum 

X 

Fe  (est.Vi/2) 

+ 

Al  (est.Vi) 

-Cu  Cyl  Calc. 

—  -  -  Fe  Cyl  Calc.1 .9 

-  -  -Al  Cyl  Calc. 
- Ta  Plate  Calc. 

O  Cu  (BN03-5) 

□  Steel  (BN03-5) 
A  Ta  (BN03-5) 

O  Al  (BN03-5) 


Ratio  of  Cylinder  Wall  to  Explosive  Radius  ( tcyi  /  Rex )  or 
Plate  Thickness  to  Sandwich  Explosive  Thickness  ( tp  /  Tex  ) 


r  =  2.75  for  Calc. 


New  formulas  improving  insight  into 
detonation-driven  propulsion 


Initial  velocity  imparted  to  cylinder  wall 
Vi  /  D  =  0.2085  (Pcyl  /  pex)  - (tcyl  /  Rex)  (-  3'4°> 

Equation  for  cylinders  describing  ratio  of 
velocity  imparted  by  initial  coupling  to  final  velocity 

Vi  /  Vf  =  0.3446  [r  -  1]  (tcyl  pcyl/  Rex  Pex)  <5/3°>  (tcyl  /  Rex)  <-3/40> 

[(tCyl  I  Rex)  2  +  2  (tcyl 1  Rex)  +  0-5  (pex  1  Pcyl)]  1/2 


Vi  =  initial  free-surface  velocity 
Vf  =  final  “steady-state”  velocity 


Transition  Pressure  (Gpa) 

Aluminum 

20.5 

Carbon  (pressed  graphite) 

23 

Iron  (0.2  wt%  Carbon) 

14.7 

Iron  (0.5  wt%  Carbon) 

13 

Titanium 

9.4 

Velocity  /  Final  Velocity  ("Steady-State") 


Normalized  Velocity  Data  for  Eight  Explosives 

Driving  Copper  Cylinders 


Equation  for  instantaneous  velocity  (Ver)  as 
a  function  of  the  gaseous  detonation  products 
expansion  ratio  (ExR) 

Ver/Vf  =  VI  /V f  +  V2/Vf 

Vl/Vf  =  (Vi/Vf)[(e-ExR-e-ExRA3)  + 

(  {  ExR  /  ExRf  }  - 1/2  /  [  1  /  ExRf  ]  “ 1,2  )  ] 

V2/Vf  =  [1  —  {  (Vi  /  Vf)  /  (1  /ExRf)-1'2}] 

(ExR  /  ExRf)  - 1,3  [  log  (ExR)  /  log  (ExRf)  ] 

Vi  =  initial  free-surface  velocity 

Vf  =  final  “steady-state”  velocity  at  final  expansion  ratio  (ExRf) 


(Equation  was  fitted  “by  eye”  using  MathCad™  software.) 


Gurney  Velocity  /  Detonation  Rate 


0.45 


Normalized  Gurney  Velocity  Data  for  Some  Explosives 
in  Cylinders  and  Spheres  of  Different  Materials 
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0.35 


0.30 
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X 

X 


Comp  B  in  Al  2S  Cyl 
Comp  B  in  Al  Sphere 
TG  50/50  in  Al  Sphere 
PBX9404  in  Al  Sphere 
LX-04  in  A  Sphere 
RDX/TNTin  A  Sphere 
Comp  B  in  Cu  Cyl 
PBX9404  in  Cu  Cyl 
LX-04  in  Cu  Cyl 
TG  50/50  in  Cu  Cyl 
RDX/TNT  in  Cu  Sphere 
Comp  B  in  Steel  Cyl 
LX-04  in  Mild  Steel  Cyl 
LX-04  in  Stainless  Steel  Cyl 
Unk  Expl,  Stainless  St.  Sphere 
TG  50/50  in  Steel  Sphere 
RDX/TNT  in  Zn  Sphere 
RDX/TNT  in  Pb  Sphere 
Unk  Expl  in  U  Sphere 
Calc.  Comp  B  in  A  Cyl 
Calc.  Comp  B  in  A  Sphere 
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♦X 
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0.01 


0.10 


1.00 


Normalized  Boundary  Areal  Density 


The  important  message  is . 

material  properties  and  geometry 
affect  Gurney  Velocity  measurement 

The  trend  of  these  effects  was  demonstrated  in: 

•  Cylinder  tests  (copper,  aluminum,  and  steel) 

•  Fragmentation  tests  (copper,  aluminum,  etc.) 

New  data  and  analysis  are  needed: 

•  Published  data  were  not  specifically  taken 

to  reveal  the  effects 

•  Most  copper  cylinder  test  data  taken  at 

(^cyl  Pcyl  ^  ^ex  Pex)  —  ^ 

•  Gurney  (Lagrange)  assumptions  may  not  be  valid 

New  opportunities: 

•  Gurney  values  would  reflect  materials  &  geometry 

•  To  define  effect  of  cylinder  wall  phase  transitions 


Another  important  message  is . 

that  gas-dynamic  propulsion  by 
many  explosives  is  similar 


Differences  between  explosives  arise  from: 

•  Vi  (solid-state  detonation  coupling  effects) 

•  Total  convertible  energy  (i.e.  Vg  from  Vf  at  ExRf) 

•  Time  needed  to  convert  total  energy  to  kinetic 

energy  of  boundary  materials  &  detonation  products 

New,  more  accurate  data  are  needed: 

•  Vi  [grazing  wave]  using  interferometry  and 

various  materials  subject  to  phase  transitions 

•  Vf  data  at  larger  ExRf  [beyond  6.5  to  10  - 14] 

•  Not  sure  that  Gurney  (Lagrange)  assumptions  valid 
[need  more  data  from  Vi  to  ExR  =  2  &  beyond] 


Additional  Backup  Information 
Follows . 


Pressure  (JWL  Along  Isentrope)  Compared  To  Acceleration  Equation 


Comparison  of  JWL  Calculated  Pressures  to  the  Time  Derivative 
of  the  Instantaneous  Velocity  Equation  for  an  Expanding  Cylinder 


Velocity  /  Final  Velocity  ("Steady-State") 


Non-Dimensional  Velocity  as  a  Function  of  Explosive 
Volume  Expansion  for  Sixteen  "Ideal"  Explosives 


Expansion  Ratio  (Volume  /  Initial  Volume) 


Velocity  /  Final  Velocity  ("Steady-State") 


1.05 


Non-Dimensional  Velocity  as  a  Function  of  Explosive 
Volume  Expansion  for  Five  Commercial  Explosives 


0.95 
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-  -  -  0.39_ExpR10 

■  -  0.45_ExpR8 


5  7  9 

Expansion  Ratio  (Volume  /  Initial  Volume) 


3 


11 


13 


Currently  available  trend-line 

Vg  /  D  £  (0.605  /  [T-1  ])  (tcyl  pcyl  /  Rex  pex)  <-  5/30> 

(J.  Roth  per  J.E.  Kennedy)  --  expanded 

For  use  with 

Vf  =  Vg  (pcy/PeX)  [(tcy/Rex)2  +  2  (VRex)  +  °'5  (Pex^Pcyl)]  ' 


Future  work 

Vg  /  D  £  (0.60  <|>  - 1/2  +  0.648  p0  1/2)  /  (1 .01  +  1.313  p0) 

(Hardesty  &  Kennedy  /  Kamlet  &  Hurwitz) 


Vg-p  /  D  =  (0.541  /  [r-1])  (tcyl  Pcyl  /  Rex  Pex)  5/30» 

(Gas  -  push  contribution  model) 
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Overview 


❖  Objectives 

❖  XM-982  -  Excalibur  Projectile 

❖  Challenges 

❖  SRVFEM  Models 

❖  Results  &  Impacts  to  Excalibur  Program 

❖Validation  &  Verification. 

❖Utility  of  3D  FE  Models. 

❖  Conclusions  &  Recommendations. 
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Objectives 


155-mm,  extended-range,  guided  munition 


■  To  develop  and  validate  a  3D  FE  model  capable  of 
simulating  launch  environment  of  the  Excalibur 
projectile. 

■  To  predict  component  loads  and  support  the  design  of 
artillery  launched  guided  projectile  components. 

■To  support  failure  investigation  of  critical  projectile 
units  in  an  evolving  design  environment. 
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XM-982,  Excalibur  Artillery  Launched 

Guided  Projectile  ' 


Demanding  g 
Requirements 

■  In  Bore:  Axial-  15,800-gs 

■  Muzzle  Exit:  Axial  -  4,052-gso<!|, 

and  balloting  3,962-gs 


Design  Development  Methods 

■  Limited  field  tests 

1  Virtual  simulation  of  launched 
environment  using  a  detailed  3D  FEM. 
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ISSUES 


Increasing  use  of 
sensitive  electronics  in 
a  guided  projectile. 


Survivability  of  MEMS 
is  a  major  concern. 
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Challenges 
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■  Obtain  physics-based  representation  of  a  highly  complex  structural 
configurations  &  interactions 

■  Provide  real-time  design  driving  guidance  to  evolving  development  of 
projectile’s  structure 

■  Achieve  “fast-turn-around”  version  of  a  representative  3D  transient  model 
inherently  “lengthy” 

■  Obtain  accurate  &  representative  validation  test  data-  difficult  to  keep  up  with 
the  test  program 

■  Address  “test-based”  malfunctions  or  failures 
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OBR  Measured  Acceleration  (G) 
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rtillery  Launched  Guided  Projectile  — 

Typical  Launch  Response  Environment 


1  .5 


x  1  0 


Truck  2C--OBR  Muzzle  Ex  it 


0  .  5 


-0  .  5 


-1  .5 


H 


AXIAL  Response 


"H™, 


Radial 


Responses 


It-  rT7  :  ^  TT.  -n  ■—  '  ^  — 


0.002  0.004  0.006 


0.0  0  8  0.1 
Time  (sec) 


Entire  projectile  inside  gun  tube 
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0.012 


Muzzle 

Transition 


A 


0.014  0.0 


Muzzle 
Base  Exit 


ia  I 


d  ia  I  1 
d  ia  I  2 


0.01 


Demanding 

Environments 

■  Accel. 

~  16  k’g  -  Barrel  axial 

-  4  k’g  -  Balloting 

-  4  k’  g  -  Exit  axial 

■  Pressure 

~  52  ksi  -  Base 

■  Frequency  (Hz) 

33  kHz+  (exit) 


Strength  Based  design  of  subsystems 


S  Uncertain  Design  Environments 


■  High-g  transient  dynamics 
response  based  design  of 
guided  electronics  &  their 
mountings 
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3D  FEM  of  SRV  Projectile  - 

Modeling  Barriers  &  Issues 


•  Geometry 


efeatured 


-  Degree  of  defeaturing 

-  Degree  of  compatibility  of  CAD  and  CAE  platforms 

-  Art  and  Science 

Material  Models 

-  Inexact  material  properties  used  in  projectile 

-  Correlate  with  response  data 

Internal  constraints 

-  Ties  vs  contact  of  parts  within  proj 
gun 

-  Joints  and  transfer  of  loads 

Loads 

-  Base  pressure  time  hist< 

-  Gun  &  geometry 

Boundary  Conditions 

-  Trunions’ effect 
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3D  FEM  of  SRV  Projectile  - 

Modeling  Cases  &  Conditions 


❖  Over  50  modeling  cases  were  run  including: 

•  Projectile  systems  &  variations 

•  Subsystem/components  -  imu,  gnu  etc. 

•  Gun  barrel  interactions 

•  Spinning  and  balloting  effects 

•  Verification/validation/model  consistency  check 

•  Sensitivity  analyses 

•  Joint  compliance/joint  loads 

•  Gravity  gun  droop 

❖  A  selected  subset  of  cases  is  discussed  next. 
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3D  FEM  of  SRV  Projectile  - 

NO  Gun  Barrel  in  Model  (old  SRV) 


Restraints  (pinned),  u2  =  u3  =  0 


rYN  BAPPEA 


Simulated  Boundary  Conditions  (Progressive). 

Step  1 .  Both  obturator  rings  are  restrained  along  radial  directions 
Step  2.  Only  back  ring  is  restrained  along  radial  dir. 

Step  3.  Free-  BCs 


Corrected  base  pressure  with  exponential  smooting  at 
muzzle  exit 


P  =  10,860* 


Time,  sec 


Bp1-Step3-Tie-P76-V5 


High-g  impulsive 
thrust  (Typical) 


Projectile-only”  FEM  is  capable  of  predicting 
Design  loads  for  most  subsystem  design. 


#  DOFs 


#  Elements 


Total  weight  in  tbs 


C.G.  from  the  bottom  of  Boom 
Base  in  in. 


CPU  Time  in  sec,  8  CPUs 


556n  680 


123,399 


104.58 


16,87 
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3D  FEM  of  SRV  Projectile  - 

Gun  Barrel/Projectile  Interaction  Model 


Two-phase  FEA  Analyses 

❖  Phase  1.  ABAQUS  Standard  gravity  analysis 

❖  Phase  2.  ABAQUS  EXPLICIT  Dynamic  Analysis 


Model 

H  m  £R  Y-Sp  to  irri  pwt- W  djt 

#DOFs 

1,124,979 

#  Elements 

374,992 

Total  weight  in  lbs  (Projectile) 

105.78 

C.G.  from  the  bottom  of  Base  in 

in. 

13.78 

CPU  time  r.t.  No  bar 


-i- 


3-  4  X 
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Sample  Analysis  Results  3D  SRV  FEM- 


Components  shock  environment 


Viewport:  1  ODE:  E  :.'ChowG tenp-ocb/ternp-hpc >11  ew-srv-import-644.odb 


S,  Rises 

(Ave.  Crit.:  75* ) 


+1 . 19Se+04 
+9.95 4e+03 
+7 . 963e+03 
+5 . 973e+03 
+3 . 982e+03 
+1.991e+03 
+4.391e-03 


Step:  Preload  Frame:  □ 


Viewport  2  ODE:  E : /ChovjG/tenTp-odb.'ternp-hpcll ew-srv-irrport-6 44.ocfe 


Viewport  3  ODE:  E : /ClwvvG.'tenTp-odb.'tenTp-hpc'TI ew-srv-inport-6 44.ocb 


S,  Rises 

(Ave.  Crit.:  75*) 


+1. 19Se+04 
+9. 954e+03 
+7. 963e+03 
+5. 973e+03 
+3. 932e+03 
+1. 991e+03 
+4. S91e-03 


Step:  Preload  Frame:  □ 
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3-D  FEM  of  SRV  Projectile  - 


With  and  Without  Gun  Barrel  —  advantages  &  disadvantage 


M 


•  Gun  Launch  Simulation  FEM  with  Barrel 

-  Increases  the  computational  cost;  dof  s  are  almost  twice 
the  no-gun  simulation  case 

-  Gun  barrel  interaction  phenomena  is  critical  to 
electronics’  survivability  analysis 

-  High-g  transient  load  for  design  of  embedded  electronics 
components  can  be  predicted  using  this  FEM. 

•  Gun  Launch  Simulation  FEM  without  Barrel 

-  Design  loads  for  stress  adequacy  of  subsystems  can  be 
effectively  predicted. 
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FEM  Validation  -  during  launch  process  in  155-mm 

Comparison  of  Simulation  vs.  Test  response  at  OBR  Location 


"i 


Predicted  and  test  responses  are  reasonably  well  correlated  for  the  event 
BP1  Configuration  with  NO  gun  FEM 
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BARREL  TIP  DISP  (IN) 


Multiple  FEM  Verification  —  barrel  projectile  model^^,^— 

Comparison  ot  Simulated  Spin  vs.  No  Spin  Responses  using  Different  l-t  Codes  1 


" 

:  C 

J  U  - 

- NASTRAN-NOSPIN-D1244_Z-Comp. 

- NASTRAN-SPI N-D 1 244_Z-Comp. 

- ABAQUS-NO-SPIN-G-DZ 

- AB AQ US-G-DZ-spin,  in.  \ 

— I — A^AQUS-G-DZ-Revised-spin,  in.  \. 
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1  Predicted  Barrel  Tip  Dis.  (Grav.  Dir) 


-0.0100 


0.005  0.010 

TIME  (SEC) 


Predicted  Barrel  Tip  Dis.  (Perp.  To  Grav.) 


-0.005 


0.000 


0.005  0.010 

TIME  (SEC) 


0.015 


0.020 


-200Ga 


NASTRAN  FWD  BOR-X 
—  NASTRAN  FWD  BOR-X.-20hz 
ABAQUS-NO  SPIN-FB-AX 
-  SIMBAD-SPIN-FB-AX,  FB-EAX 
ABAQUS-FB-AX-spin,  g's 
ABAQUS-FB-AX-Revised-g-spin,  g's 


Qualitative  agreement 
among  predicted 
results  using  various 
FE  codes  verifies  the 
effectiveness  of  the 
FEM. 


-10000 


time,  sec 


NASTRAN  FWD  BOR-Y-20hz 
ABAQUS-NO-SPI N-FB-AY 
ABAQUS-FB-AY-spi  n,  g's 
ABAQllS-FB-AY-Revised-spin,  g's 


No  Spin 


Predicted  Radial  Accel. 


TACOM 


ARDEC 


f 


1— 


ZCAmUR- 

“14 


XM982 


Acceleration,  axial,  g's 


I  Itilitv  Qt  3D  Transient  SRV  FFM  - 


jpfcj 


w 


Locating  suitable  mounting  locations  of  sensitive  electronics  components 


"i 


0.015 


Time,  sec  *•. 


Viewport:!  Mv.*<el:  0br2-p76data-mo  del-4  Step-  .litial 


□  Comparison  of  simulated  responses: 


G’s  at  Forward  location  -  8,500 
G’s  at  AFT  location  -  3,000 


Demonstrates  the  effectiveness  of  FEM  in  evaluating  design  concepts. 
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Utility  of  3D  Transient  SRV  FEM  - 

Sub  Modeling  Analysis  -Micro-Level  failure  Investigation  of  IMU 


Viewport  1  ODB:  d:JMRCfSRV-FE453XAftJ  MU_..  JA  FT  AM  U-abq&31JJ  ob-l.od  b 


ewport:  1  Model:  0br2-p76data-m o del-d  Step:  Initial 


Z1  IQ  :  06;  OZ  Eastei:; 


Top  Driven  Nodes  from 
Global  Model 


Time  Saving 
Estimate  >  1 0X 


Predicted  global  Accel,  for  driven  nodes  are  used  for  sub  modeling 
evaluation  of  the  IMU  device. 


Local  detail  FEM  of  the  IMU  is  used  thereby  saving  computational  time. 


TACOM 


ARDEC 


fcCAlim 

L  .i  r 


XM982 


■*16 


Utility  of  3D  Transient  SRV  FEM  - 


Failure  Investigation  of  MEMS 


■I 


Acceleration  Time 
history  at  location  of 
interest 


Dynamic  Response 
(extract  dynamic  load') 


SRS  (Shock 
Response  Spectrum) 


Stress  calculations 
based  on  SRS 


Component 
Survivability  ?? 


a 
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"\f 

B  (Qs) 
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Dynamic  only  (A-B) 


PickSDOF 
response  for 
freq.  of  interest 
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predicted  ,max 


<  (7 


allowable 


Embedded  MEMS 
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Predictive  Capability  of  SRV  FEM- 

General  behavior  is  explained 
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Predicted  Displacement  &  Velocity  of  SRV 
Projectile 


- Axial  Displ.  (U3) 

—  Axial  Vel.  (V3) 


Time,  sec 


Predicted  Contact  Forces  (CF)  of  Gun 
Launched  SRV  projectile 
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CFNl(Front) 

■  CFN2  (Bottom) 
CFS3  (Front) 

■  CFS3  (Bottom) 


Time,  sec 


Conclusions  &  Recommendations 


❖  Virtual  simulation  of  a  system  behavior  using  a 
validated  and  verified  FEM  provides  a  flexible  and 
powerful  design  evaluation  tool  for  the  future 
sophisticated  systems. 

❖  Gradual  development  of  FEM  starting  with  a 
simple  model  and  then  progressively  adding 
features  to  improve  the  model  capability  is 
important  in  the  development  of  a  complex  system 
such  as  the  SRV  projectile  discussed  here. 

❖  Validation  of  FEM  with  measured  results  is  desired, 
however,  verification  with  alternative  predictive 
tools  may  be  used  as  an  alternative. 
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Clay  block 


Fabric 


Gun 
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Steel  Box 


Polycarbonate  window 


Flash  Unit 


Carrera 


MANCHESTER 

1824 


MANCHESTER 

1824 


http://www.btintemet.com/~warship/Today/wave.htm 
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Impact 
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Multi-layer  areal  density  required 
to  prevent  penetration  (kg  m  "2) 


■  Fabric  F10 


O  Fabric  F11 


Weight  of  fabric  required  to  prevent  penetration  vs.  width  of 
fabric  strip. 


Thread 

Density 


Picks/ 

Ends  / 

Fibre 

cm 

cm 

Kevlar 

129 

6.7 

6.63 

Kevlar 

49 

22.4 

22.4 

Linear  density 
of  yarn 

Weft  Warp 


(tex) 

(tex) 

Structure 

Plain 

167 

167 

weave 

Plain 

23 

23 

weave 

Multi-layer  areal  density  require 
to  prevent  penetration  (kg  m  ' 


O  Fabric  F11 

_ 


Weight  of  fabric  required  to  prevent  penetration  vs.  width  of 
fabric  strip 


Width  of  strip 

Weight 

of  fabric  required 
to  prevent  penetration 

(cm) 

(kg/sq  m) 

24 

1.12 

3.3 

0.37 

Weight  saving  (%)  | _ 67 _ 

Example  of  the  effect  of  reducing  fabric  width  for  direct 
impact  onto  a  clamped  narrow  fabric. 


> 


X 


Clay 


Fabric 


c 


Fabric  width  (cm) 


Weight  of  fabric  required  to  prevent  penetration  vs.  width  of 
fabric  strip 
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a)  cut  edges 


b)  conventional  selvedges 


Schematic  of  cross-section  through  narrow  fabrics 


Multi-layer  areal  density  required 
to  prevent  penetration  (Kg  m  *2) 


F10  -  Cut  edges, 
S  =  0.99 

FI 2  -  selvedges, 
S  =  0.63 


Weight  of  fabric  required  to  prevent  penetration  vs.  width  of 
fabric  strip 


Mechanics  of  impact.  Initially  narrow  fabric  deformations  are 
similar  to  those  found  for  a  regular  fabric  panel. 


Regular  fabric  panel 


Narrow  fabric 


Later,  in  the  case  of  narrow  fabrics,  deformation  extends 
further  in  a  direction  parallel  to  the  long  sides.  Typically, 
a  peanut  shaped  deformation  is  observed. 


This  is  reflected  in  the 
residual  deformations. 
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There  are  two  problems  in  making  practicable 
armour  from  narrow  fabrics: 

1 .  The  weakness  of  the  joins  between  fabric  strips. 

2.  The  weight  and  stiffness  of  any  structure  for 
holding  the  fabric  on  two  sides. 
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Narrow  fabrics  can  be  offset  in  different 
layers  but  the  joins  still  contribute  to  an 
overall  weakening  of  the  fabric  panel 
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Holding  the  fabric 


Various  structures  for  holding  the  fabrics  were 
investigated  but  were  eliminated,  either  because  of 
weight,  because  the  grip  was  insufficient,  or  because 
they  were  destroyed  by  impact. 
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The  solution  was  to  use  expanded  polystyrene  board. 
The  fabric  strips  were  simply  glued  along  the  short  sides 


Polystyrene  board  (Weight  =  28g)  Glued  edge 


Face  and  back  of  the  armour  panel. 
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samples 


Crater  Depth  in  clay  backing 


L 


J 


Not  only  must  armour  arrest  a  projectile,  but  it 
must  do  it  without  excessive  deformations.  For 
example,  it  is  of  no  use  stopping  a  projectile  if 
it  takes  a  metre  to  do  so. 


Crater  Depth  in  clay  backing 


Crater  Standard 

Depth  (cm)  deviation  (cm) 


Benchmark  (Fabric  F10)  1.6  0.6 


Narrow  fabric  assembly 

(Fabric  Fll)  |  1.4  | _ 0.2 

Crater  depth  for  narrow  and  benchmark  fabric  of  equal 
areal  densities. 


Performance  of  narrow  fabric  assemblies 


Multi-layer  areal  density  required  to 
prevent  penetration 


Impact 

velocity 

Benchmark  - 
Fabric  F10 

Narrow  fabric 
assembly  -  Fabric  Fll 

Weight 

reduction 

(m  s_1) 

(kg  nr2) 

(kg  nr2) 

(%) 

400 

1.6 

1.1 

31 

520 

3.4 

3.2 

6 

Weight  saving  -  Narrow  fabric  panel  compared  to  control 
fabric.  Results  based  on  weight  of  whole  assembly. 


Performance  of  narrow  fabric  assemblies 


Multi-layer  areal  density  required  to 


prevent  penetration 


Impact 

velocity 

Benchmark 
Fabric  F10 

Narrow  fabric 
assembly  -  Fabric  Fll 

Weight 

reduction 

(m  s_1) 

(kg  nr2) 

(kg  nr2) 

(°/o) 

400 

1.6 

0.6 

63 

520 

3.4 

2.6 

24 

Weight  saving  -  Narrow  fabric  panel  compared  to  control 
fabric.  Results  based  on  weight  of  fabric  alone. 
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What  we  do  not  know  about  narrow  fabrics 


The  minimum  forces  required  to  grip  the  fabrics  (Can 
the  weight  of  a  soldier's  equipment  be  utilised  to  hold 
the  narrow  fabrics  under  tension?) 

The  weave/width  combination  required  to  maximise 
ballistic  performance. 

The  selvedge  construction  required  to  maximise 
ballistic  performance. 
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What  we  do  not  know  about  narrow  fabrics  (Continued) 

The  fibre  properties  required  to  maximise  the  ballistic 
performance  (Should  they  be  different  in  warp  and  weft). 

The  frictional  properties  of  yarns  required  to  maximise  the 
ballistic  performance. 

The  relationship  between  ballistic  performance  and  fabric 
length. 

The  relationship  between  ballistic  performance  and  impact 
velocity. 


What  we  know  about  narrow  fabrics 

The  deformations  of  narrow  fabrics  undergoing  ballistic 
impact  differ  to  those  for  regular  fabrics. 

With  the  appropriate  choice  of  narrow  fabric  construction, 
width,  selvedge  and  method  of  gripping  significant 
improvements  in  ballistic  performance  can  be  achieved. 

Incorporating  narrow  fabrics  into  ballistic  panels  results  in 
improvements  in  performance-to-weight  ratios  despite  the 
weaknesses  introduced  by  fabric  joins. 

Where  some  existing  structure  is  present,  such  as  an 
aircraft  fuselage,  still  greater  performance-to-weight  ratios 
can  be  achieved. 
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Outline 


•  Introduction 

•  Numerical  Models 

•  Preliminary  Calculations 

•  Cylinder  Expansion  Test  Simulations 

•  Discussion 

•  Future  Work 
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Introduction 


•  Cylinder  Expansion  (CYLEX)  tests  conducted  to 
determine  explosive  performance 

-  New  conventional  explosives 

-  Non-standard  mixes 

•  Explosive  cased  in  a  copper  cylinder  is  detonated 

•  Resultant  expanding  wall  velocities  recorded 

•  Explosive  energy  can  be  determined 

•  Allows  for  determination  of  parameters  for  detonation 
product  gas  equation  of  state 
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Cylinder  Expansion  Test 


Images  Courtesy  of  DRDC  Suffield 


Objectives 


•  Simulation  of  cylinder  expansion  tests  for  three 
increasingly  energetic  explosives: 

-  Nitromethane 

-  TNT 

-  Composition  B3  (60  wt%  RDX,  40  wt%  TNT) 

•  Validation  of  Chinook  CFD  code 

-  Multiple-material  model 

-  Reaction  model 


CHINOOIC 


MARTEC 


Background 


•  Multiple-material  model  previously  applied  to  (non¬ 
reactive)  underwater  explosion  and  landmine 
simulations 

•  One  of  earliest  calculations  after  implementation  of 
advanced  energetic  material  modelling  capability  in 
Chinook 

•  Calculations  performed  using  Eulerian  formulation 
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Model  Theory 


•  Continuum  approach  -  mixed  cells  with  uniform 
pressure  and  velocity 

•  Allows  for  several  equations  of  state  to  be  used  in  a 
calculation 

-  Mie-Gruneisen  (HOM)  for  condensed  explosive 
and  tube  wall 

-  JWL  for  explosive  products 

-  Ideal  gas  for  air 

•  One-step  Arrhenius  reaction  model  converts 
reactants  into  products 

-  User-specified  heat  of  reaction 
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Preliminary 

Calculations 


•  0.1  m  long  ID  domain 

•  Small  initiation  region  at  CJ  state 

•  0.2  mm  cell  size 

•  CJ  pressure  and  detonation 
velocity  extracted  for  each 
explosive 


Explosive 


Nitromethane 


TNT 


Comp  B3 


Detonation 
Velocity  (m/s) 

6800  (2.8%) 

7280  (1.3%) 

8400  (3.1%) 


CJ  Pressure 
(GPa) 

13.27  (0.84%) 
19.97  (-1.8%) 
27.03  (-0.66%) 
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Cylinder  Expansion  Simulations 


•  2D  axi-symmetric  analysis  of  cylinder  tests 

•  30  cm  long  by  2.54  cm  diameter  copper  tube 

•  0.2  mm  cell  size  (structured  quads) 

•  Expanding  grids  eliminate  boundary  reflections 

•  Approximately  1.1  million  total  cells 

-  Utilize  parallel  computing 
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Cylinder  Expansion  Simulations 


•  Monitoring  array  placed  20  cm  down  tube 

•  Wall  velocities  extracted  at  6  and  19  mm 

•  Wall  position  taken  at  0.5  mass  fraction  copper 


Wall  Displacement 


MARTEC 


t=  0.00000  ms 


8 

i 
■  ■■ 


(0 
■  Hi 

■o 

(0 


Axial  Distance  (m) 


Pressure  Contours 


Comparison  of  Results 


E 


cylinder 


(mm/us)2/2 
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Discussion  /  Future  Work 

•  Reasonable  agreement  with  experimental  results 

-  CJ  pressure  /  detonation  velocity 

-  Wall  velocities 

•  Higher  mesh  resolution  ? 

•  Eulerian  approach  results  in  some  numerical  diffusion 

-  Adaptive  mesh  refinement 

-  Material  interface  tracker 

-  Lagrangian  or  moving  mesh 

•  Next  Step:  Revisit  this  study  and  investigate 
explosives  with  higher  reaction  rates 
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Material  for  Explosively  Formed 
Penetrator  (EFP)  Warhead 
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Outline 

•  Background 

•  Consolidation 
processes 

-  Hot  Isostatic  Press 

-  Extrusion 

-  Resistance  Sintered 

•  Near  Net  Process 

•  Material  Evaluation 

•  Summary 
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Explosively  Formed  Penetrator 

(EFP)  Warheads 


Single  EFP 


Warhead  hardware 


Multiple  EFPs 


Combined  Effects  EFP 


Conventional  EFP 
Fabrication  Process 


c 


i  i  i  i  i 

Ta  Bar  Cut  Into  Slugs 


c 


FORCE 


Slug 


STEP  1 


Forged  Disc 


STEP  2 


Anneal(Temperature  Treatment) 

i  i  i  1  i 


STEP  3. 


Machine  Domed  Disk 
With  Desired  Liner 
Surface  Profile 


STEP  5 J 
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EFP  Formation  Variation 


f  2nd  Cut  From  Ta  Bar  \ 

CDOCD  (D  CD  (D  CD  (D  CD  (D  CD 


/  7th  Cut  From  Ta  Bar  \ 
DDDDDDCDDDD 


9th  Cut  From  Ta  Bar 
DDDDDDDDODD 


EFP  Warhead  Powder 
Tantalum  Program 


•  Objectives:  To  develop 

•  Powder  process  specification 

•  Near  net  shape  liner  fabrication  process 

•  Advantages  of  powder  Tantalum  approach 

•  Availability 

•  Consistency 

•  Lower  costs 
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Powder  Consolidation  Processes 

•  Hot  Isostatic  Press 

•  Extrusion 

•  Resistance  Sintered 
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HIP  consolidation 
process 
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HIP  Test  Results 


Orthogonal  X-Ray  Results 


Horizontal 


Cordin  Photo  Results 


Vertical 
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Extrusion  consolidation  process 


RAW  POWDER 

j) 

COLD  ISOSTATIC  PRESS 

EXTRUSION 


FINISHED  BAR  STOCK 
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Extrusion  Test  Results 


Vertical 


Orthogonal  X-Ray  Results 


Horizontal 


Cordin  Photo  Results 


Resistance  Sintered 

consolidation  process 

+ 

RAW  POWDER 

V 

COLD  ISOSTATIC  PRESS 

V 

RESISTANCE  SINTERED 

V 

FINISHED  BAR  STOCK 
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Resistance  Sintered  Results 


Horizontal 


Vertical 


Orthogonal  X-Ray  Results 


Cordin  Photo  Results 
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Comparison  of 
Consolidation  Processes 


Near  Net  Process 
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Near-Net  Process  -  Powder 

To  “Pucks” 


Hot  Isostatic  Press  Insert 
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HIP  To  Puck  Process 

Pucks  in  Container 


Pucks  (Post-HIP)  Containers  removed 


HIP 
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Test  Data  (HIP) 


HIP  to  PUCK 


X-Ray  Results 


Cordin  Photo  Results 


HIP  to  Bar 


EFP  Design  and  Material 

Evaluation 
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EFP  Design  Progression 


Design  #1 


Design  #2 


Design  #3 
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Material  Evaluation 


To  assess  Poweder  Ta  ductility  for  advance  EFP  designs 


Advance  EPF  Design 


EFP  with  pure  Ta 


Hydrocode  Simulation 


EFP  with  TaW  alloy 
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Material  Evaluation  -  Con’t 


Powder  Ta  liner  shot  #1 


Powder  Ta  liner  shot  #2 


Powder  Ta  liner  shot  #2 


X-Ray 


X-Ray 


Cordin  Photo 


Pure  Ta  Bar 


X-Ray 
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Summary 


•  Demonstrated  lower  cost  repeatable 
EFP  Liner  fabrication  process  using 
Powder  Tantalum. 

•  Eliminated  lot-to-lot  variations  in 
material  property 

•  Generated  a  powder  tantalum  material 
process  and  liner  fabrication 
specification 
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METAL  SOLUTIONS  FOR 
TACTICAL  VEHICLES 


WEAPONS  AND  MATERIALS  RESEARCH  DIRECTORATE 


ARMOR  MECHANICS  BRANCH 


•  MANY  CURRENT  ARMORED  KITS  FOR  TACTICAL 
VEHICLES  UTILIZE  ROLLED  HOMOGENEOUS 
ARMOR  OR  HIGH  HARD  STEEL  SOLUTIONS 

•  SELECTION  HAS  BEEN  DRIVEN  BY  EXPEDIENCY 
REQUIREMENTS,  MATERIAL  AVAILABILITY  AND 
COST  CONSIDERATIONS 

•  METAL  SOLUTIONS  ALSO  PROVIDE  VERY  GOOD 
MULTIHIT  FRAGMENT  PROTECTION  AGAINST 
IMPROVISED  EXPLOSIVE  DEVICES  DUE  TO  THE 
TOUGHNESS  OF  THE  PLATES 

•  ARL  IS  RE-EXAMINING  METAL  SOLUTIONS  THAT 
COULD  PROVIDE  REDUCED  WEIGHT  AND/OR 
IMPROVED  PROTECTION 


POSSIBLE  METAL 
TECHNOLGIES 


WEAPONS  AND  MA  TERIALS  RESEARCH  DIRECTORA  TE  ARMOR  MECHANICS  BRANCH 

•  IMPROVED  STEELS 

-  ULTRA  HIGH  HARD  STEELS 

-  DUAL  HARD  STEELS 

•  HIGH  STRENGTH  ALUMINUM  ALLOYS 

•  ALTERNATE  TITANIUM  ALLOYS  TO  TI-6AL-4V 

-  BETA  ALLOYS 

-  DUAL  HARD  TITANIUM 

•  METAL  LAMINATES 

-  MILD  STEEL/RHA  OR  HIGH  HARD 

-  STEEL/ALUMINUM  LAMINATES 

-  TITANIUM/ALUMINUM  LAMINATES 

•  METAL/POLYMER  COMPOSITES 

-  METAL/POLYMERIC  LAMINATES 

-  METAL/FIBERCOMPOSITES 


DUAL  HARD  STEEL 
DEVELOPMENT 
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•  DUAL  HARD  STEEL  WAS  DEVELOPED  IN  1965  AT 
THE  WATERTOWN  ARSENAL,  NOW  ARL,  FOR 
HELICOPTER  ARMOR  APPLICATIONS 

•  MIL-S-46099A,  STEEL  ARMOR  PLATE,  ROLL- 
BONDED,  DUAL  HARDNESS  WAS  ESTABLISHED  IN 
NOVEMBER  1966 

•  PRODUCTION  FEASIBILTY  STUDIES  FOR  ROLL 
BONDED  PLATE  WERE  CONDUCTED  BY  US  STEEL 
IN  1968  AND  THEN  AT  OTHER  COMPANIES 

•  ALLEGHENY  LUDLUM  IS  THE  ONLY  CURRENT 
PRODUCER  OF  K12®  DUAL  HARD  ARMOR  STEELS 
USING  A  PROPRIETARY  PROCESS 


ALLEHGENY  L  UDL  UM  K12® 
DUAL  HARD  STEEL 
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•  K12®  DUAL  HARD  STEEL  RESULTS  FROM  ROLL 
BONDING  TWO  NI-CR-MO  PLATES  THAT  PRODUCE 
A  HIGH  HARDNESS  FRONT  SIDE  AND  A  SOFTER 
BACK  SIDE 

•  THE  FRONT  PLATE  HAS  A  HARDNESS  OF  601-712 
BRINELL  WITH  THE  REAR  PLATE  461-534 

•  K12®  IS  NORMALLY  FURNISHED  IN  THE 
ANNEALED  STATE  FOR  SHAPING  AND  CUTTING 
AND  IS  THEN  FINAL  HEAT-TREATED 

•  MIL-S-46099C  HAS  A  BALLISTIC  ACCEPTANCE 
CRITERIA  WITH  PLATES  FROM  0.170”-0.330” 
TESTED  WITH  THE  0.30-CAL  APM2  AND  PLATES 
FROM  0.290”-  0.585”  WITH  THE  0.50-CAL  APM2 


RESEARCH  GOAL 
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•  DUAL  HARD  STEEL  ISSUES: 

-  DEBONDING  OF  ROLLED  PLATES 

-  PLATE  SIZE  LIMITATIONS 

-  REAR  PLATE  DEFORMATION 
ALTERNATE  FABRICATION  TECHNOLOGIES 

•  EXAMINE  MIL-46099C  FOR  POSSIBLE  ADDITION  OF  A 
NEW  HARDNESS  SPECIFICATION 

•  PAPER  PROVIDES  THE  INITIAL  RESULTS  OF  ARL’S 
EXAMINATION  OF  DUAL  HARD  STEEL  PRODUCED  BY 
THE  ELECTROSLAG  REMELTING  PROCESS  AS 
DEVELOPED  BY  THE  INSTITUTE  OF  METAL  SCIENCE 
(ISM)  OF  THE  BULGARIAN  ACADEMY  OF  SCIENCES 

•  BASELINE  RHA  AND  HIGH  HARD  WERE  ALSO 
RETESTED  AS  HISTORIC  DATA  IS  INCOMPLETE 


B  UL  GARIAN  INSTITUTE  OF  METAL 
SCIENCE  DUAL  HARD  STEEL 
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•  ELECTROSLAG  REMELTING  (ESR)  PROCESS 
WAS  HIGHLY  DEVELOPED  IN  THE  1980’S,  BUT 
WAS  NOT  WORKED  ON  AT  THE  ISM  AFTER  1991 
FOR  LACK  OF  FUNDS 

•  ARL  EXAMINED  ISM  DUAL  HARD  PLATES  IN  2004 
AND  CONTRACTED  WITH  ISM  TO  RESTART 
THEIR  ESR  FURNACE  AND  PRODUCE  A 
MODIFIED  VERSION  OF  US  DUAL  HARD  STEEL 

•  ARL  SPECIFIED  A  FRONT  PLATE  WITH  A 
HARDNESS  OF  500-560  BRINELL  WITH  THE  REAR 
PLATE  340-370  BRINELL 

•  PAPER  EXAMINES  THE  BALLISTIC  PROPERTIES 
OF  DUAL  HARD  PLATE  COMPARED  AGAINST 
RHA  AND  HIGH  HARD  STEEL 


ISM  ESR016  ELECTROSLAG 
REMELTING  FURNANCE 


W 
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DR.  STEFAN  VODENICHAROV 
INSTITUTE  OF  METAL  SCIENCE 
BULGARIAN  ACADEMY 
OF  SCIENCES 
SOFIA,  BULGARIA 

INGOTS  WERE  ESR  CAST 
FORGED  AND  ROLLED 
HEAT-TREATED  AND  CUT 
300mm  X  300mm  PLATES 

PLATES  WERE  FABRICATED  IN 
THICKNESSES  OF  5mm  TO  80mm 

ONLY  5-1 0mm  PLATE  DATA 
PROVIDED  IN  THIS  PAPER 


ISM  ESR01 6  ELECTROSLA  G 
REMELTING  FURNANCE 


SOFT  HARD 

330-370BHN520-570BHN 


TEST  PROJECTILES 
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0.30-CALIBER  44-GRAIN  7.62mmx63  (0.30  CALIBER) 
FRAGMENT  SIMULATING  APM2  ARMOR  PIERCING 
PROJECTILE  PROJECTILE 


DUAL  HARD 
TEST  DATA 
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0.30-cal  FSP  V50  VELOCITY  AT  THICKNESS  (0°  OBLIQUITY)  (m/s) 

Nominal 

Thickness 

5mm 

6mm 

7mm 

8mm 

9mm 

1 0mm 

Actual 

Thickness 

5.39 

6.45 

7.54 

8.23 

9.96 

10.95 

Front/Rear 

Hardness 

523/365 

572/375 

561/373 

537/363 

544/340 

557/337 

Bulgarian 
Dual  Hard 

708 

763 

911 

1063 

1145 

>1402 

,30-cal  APM2  V50  VELOCITY  AT  THICKNESS  (0°  OBLIQUITY)  (m/s) 

Nominal 

Thickness 

5mm 

6mm 

7mm 

8mm 

9mm 

1 0mm 

Actual 

Thickness 

5.41 

6.45 

7.90 

8.74 

9.37 

10.36 

Front/Rear 

Hardness 

523/365 

572/375 

568/366 

547/361 

531/340 

539/345 

Bulgarian 
Dual  Hard 

512 

NA 

693 

691 

703 

705 

LIMIT  VELOCITY  (m/s) 


0.30  CALIBER  APM2  @  0° 
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THICKNESS  (mm) 


0.30 
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CALIBER  FSP  @  0° 
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THICKNESS  (mm) 


10mm  TEST  PLATE 
0.30-CALIBER  APM2 


m 

w 
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FRONT 


REAR 


10mm  TEST  PLATE 
0.30-CALIBER  FSP 
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REAR 


Oilwell  Perforators: 
Theoretical  Considerations 


Brenden  Grove 

22nd  International  Symposium  on  Ballistics 
November  17, 2005 


ScHunberier 


Overview 


•  Introduction/ Background 

•  Compressible  Jet  Model 

•  Implications: 

-  Hydrodynamic  penetration  depth 

-  Impact  pressure 

-  Target  strength 

•  Conclusions 
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Introduction  /  Background 


•  Oilwell  perforators 
->  20-60mm 
shaped  charges 

•  Perforate  casing, 
cement,  formation 
rock 

•  Perf  tunnels  must 
be  clean  to  allow 
subsequent  fluid 
flow 


Scfttanterger 
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Liner  fabrication 


Pressing  force 


Completed 

liner 


Liner 

punch 

Liner 

powder 


Liner 

diebody 


Liner 

diebody' 


Completed  liner  ejected 
from  diebody 


Liner  powder  placed 
in  ejector 


Liner  pressed  with 
high  force 


formed  by  powdered  metallurgy  (P/M) 
unsintered  (green)  'I.  ‘  i 


Metal  powder  photograph: 
http:/ /www.  mpif.  org/IntroPM/mak 
e  powder,  asp  ?linkid=5 
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Introduction  /  Background 
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Jet=millions  of  discrete  powder  particles 
How  to  model  jet  penetration? 

-  Traditional  models  underpredict  penetration 
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Introduction  /  Background 


•  Jet=millions  of  discrete  powder  particles 


•  How  to  model  jet  penetration? 


•Density  decreases 
•Compressibility  (porosity)  increases 


Compressible  Jet 

Penetration  Theory 

•  Multiple  discrete  impacts;  particles  "pile  up"? 

•  Macroscopically  -  jet  compresses 

•  Incompressible  Bernoulli  not  applicable 
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Compressible  Jet 

Penetration  Theory 


V-U 


► 
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Jet  Material  Compaction  Curve 


•  Flis  &  Crilly  (18th  ISB):  compressible  target 

•  Reverse  their  analysis  here 


p0  =  initial  (distended)  jet  density 
ps  =  solid  (pore-free)  jet  density 
^  =  1  -  P0  /  P  s  =  jet  porosity 
R  =  jet  compaction  initial  resistance 

Sdrinferger 


Compressible  Jet 

Penetration  Theory 


•  Distended  jet  at  p0  traveling  at  (V-U) 

•  Compacted  to  ps ,  decelerated  to 


mass 

momentum 


/>o /»•„/- A /"V/ 

$ 

II 

1 

C! 

Poj+Po&oj)2- 

'-P\j+ P\(wlj)2 

full  compaction 


P\j  ~  Pi  j  -  Pc 
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Apply  incompressible  Bernoulli  to  compacted  jet 
impact  ScUunkerfer 


Compressible  Jet 

Penetration  Theory 


Jet  pressure 
Target  pressure 


PD  (hydrodynamic) 


pv=R+\P„(r-ufn+^ 

2  V  Y  ' 

P2,=Y,  +  '-p„U 2j 


Compressibility 


PD 

L 


trrect 


PD  (target  strength 
important) 
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Compressible  vs.  Incompressible  Jet 

Penetration  Theory 


Incompressible 


Jet  Pressure 

Hydrodynamic 
Penetration  Depth 


PD  |^7 


►  Compressible 


Pv=R+\p„(V-Uf( 

2 


PD 

L 


X  — >  (1+^) 
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ScHunberiGr 


Compressible  vs.  Incompressible  Jet 

Penetration  Theory 


•  (1+<|>)-»A, 

•  Porous  jet  model  reduces  to  incompressible 
model  at  limits 

•  At  upper  limit  (<j)=1;  A,=2),  p  zero 

•  For  a  given  length  and  macroscopic  density, 
porosity  increases  the  following: 

-  Penetration  depth 

-  Impact  pressure 
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Jet  Density,  Length,  Porosity 

Summary 

Consider  3  cases:  (Fixed  mass,  velocity,  diameter) 

pL=constant 


Case  A:  ps,  incompressible;  w/////////////.  - y 

length=Ls 


Case  B:  p0,  compressible  to  ps; 
length=L0 


►  V 


Case  C:  p0,  incompressible; 
length=L0 


►  V 
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Jet  Density,  Length,  Porosity 

Summary 


Hydrodynamic  Penetration  Depth 


Case  A: 


Case  B: 


Case  C: 


(alternative  form) 
(alternative  form) 
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ScHunberier 
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Jet  Density,  Length,  Porosity 

Summary 

Hydrodynamic  Penetration  Depth 


pdb>pdc>pda 


Scttnkerger 
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Jet  Density,  Length,  Porosity 

Summary 

Dynamic  Jet  Pressure 


Case  A: 


Pa=\p.{V-U? 


Case  B: 


PB=-ps{V-U)\\-f) 


=—pa ( V-U )2 (1+^)  (alternative form) 


Case  C: 


pc=\p*<y-u? 


Scfehnkerfcr 


Jet  Density,  Length,  Porosity 

Summary 


Dynamic  Jet  Pressure 


P  >P  >P 

1  a'1  b^1  c 


Jet  Density,  Length,  Porosity 

Summary 

Consider  3  cases:  (Fixed  mass,  velocity,  diameter) 

pL=constant 


Case  A:  ps,  incompressible;  w/////////////.  - y 

length=Ls 


Case  B:  p0,  compressible  to  ps; 
length=L0 


►  V 


Case  C:  p0,  incompressible; 
length=L0 


►  V 
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Jet  Density,  Length,  Porosity 

Summary 


•Greatest  impact  pressure 
•Shallowest  hydrodynamic  penetration 


Case  A 


•Deepest  hydrodynamic  penetration 
•2nd  greatest  impact  pressure 


Case  B 


Case  C 
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Impact  Pressure  &  Target  Strength 


What  about  target  strength? 

Consider  some  critical  point  of  interest  (Ucrit,  Pcrit) 


Characteristic 
target  property 


P 


crit 


isobars  in  (pJf  V)  coordinates 


Incompressible;  1<^<2 


Mr-ucnl)2 


Compressible 


There  is 
no  unique 
V 

v  min 
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Impact  Pressure  &  Target  Strength 

What  do  these  isobars  look  like? 


Target 

Material 

Pt 

(g/cc) 

Yt  (kbar) 

^cnt 

Ps 

(g/cc) 

R=Yj  (kbar) 

Mild  steel 

7.86 

3  (UTS) 

2Yt 

16 

0 

Aluminium 

2.7 

3  (UTS) 

tt 

cc 

tt 

Concrete 

2.25 

0.3  (UCS) 

tt 

tt 

tt 

Sandstone 

2.25 

0.6  (UCS) 

(i 

tt 

tt 
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Impact  Pressure  &  Target  Strength 

Isobars  in  (p](V)  space 


Mild  Steel 


2.0E+04  4.0E+04  6.0E+04  8.0E+04  1.0E+05  1.2E+05 


Jet  Velocity  V  (cm/s) 


Incompressible  Bernoulli 
Continuous 
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Impact  Pressure  &  Target  Strength 

Isobars  in  (p](V)  space 


Mild  Steel 


Jet  Velocity  V  (cm/s) 


Incompressible  Bernoulli 
Particulated 
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Impact  Pressure  &  Target  Strength 

Isobars  in  (p](V)  space 


Mild  Steel 


2.0E+04  4.0E+04  6.0E+04  8.0E+04  1.0E+05  1.2E+05 

Jet  Velocity  V  (cm/s) 


Compressible  (porous) 
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ScHunberier 


Impact  Pressure 

Isobars  in 


Mild  Steel 


Jet  Velocity  V  (cm/s) 


Sandstone 


Jet  Velocity  V  (cm/s) 
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&  Target  Strength 

pj(V)  space 


Aluminum 


Concrete 


10 


-  -  -  incomp,  solid 

- incomp,  particulated 

comp.  Bernoulli 
•  comp.  F&C 

l 

1 

% 

\  *\ 

2.0E+04  4.0E+04  6.0E+04  8.0E+04  1.0E+05  1.2E+05| 
Jet  Velocity  V  (cm/s) 


Scfttanterger 


Jet  Density  rho  (glee 


Impact  Pressure  &  Target  Strength 

Isobars  in  (p](V)  space 
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Conclusions  (1) 


•  Developed  treatment  of  compressible  jet 

-  Reduces  to  well-known  expressions  for  solid  &  fully- 
particulated  jets 

-  For  a  given  mass,  velocity,  diameter: 

•  Porous  (compressible)  penetrator  penetrates  deeper 
than  incompressible  penetrator  of  same  L,  p 

•  ...also,  deeper  than  shorter  penetrator  of  higher  p 

•  ...produced  impact  pressure  which  is  intermediate 
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Conclusions  (2) 


•  Looked  at  steady-state  impact  pressure 

-  Presented  isobars  in  (pjfV)  coordinates 

•  Compressible  jet  model  interpolates  between  solid  & 
fully-particulated  incompressible  jet  model 

-  Highly-distended,  low-velocity  jets  may  effectively 
penetrate  moderate-strength  geologic  targets 


This  approach,  so  far,  neglects  transients 
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Thank  You 


Questions?? 
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Oilwell  Perforators: 
Theoretical  Considerations 


Brenden  Grove 

22nd  International  Symposium  on  Ballistics 
November  17,  2005 


SCHIUpGIrflGb 


Overview 


•  Introduction  /  Background 

•  Compressible  Jet  Model 

•  Implications: 

-  Hydrodynamic  penetration  depth 

-  Impact  pressure 

-  Target  strength 

•  Conclusions 
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Introduction  /  Background 


Oilwell 

perforators  -> 
20-60mm 
shaped  charges 

Perforate  casing, 
cement, 
formation  rock 

Perf  tunnels 
must  be  clean  to 
allow  subsequent 
fluid  flow 
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Introduction  /  Background 
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Introduction  /  Background 


Liner: 


-  formed  by  powdered  metallurgy  (P/M) 


-  unsintered  (green) 


Pressing  force 


Liner  Liner 
powder  diebody 


Metal  powder  photograph: 
http:/ /www.  mpif.  org/IntroPM/m 
akepowder.  asp  ?linkid=5 


\  /  Completed 
liner 


Liner  powder  placed 
in  ejector 


Liner  pressed  with 
high  force 


Completed  liner  ejected 
from  diebody 


5  BMG 
11/28/2005 


2CHHipGlrflGb 


Introduction  /  Background 
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Jet=millions  of  discrete  powder  particles 
How  to  model  jet  penetration? 

-  Traditional  models  underpredict  penetration 
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Introduction  /  Background 


Jet=millions  of  discrete  powder  particles 


,  ,  .  i  i  ■  -  -  •tensity decreases 

How  to  model  jet  penetration? 

J  '  •Compressibility  (porosity) 

increases 
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Incompressible  Jet 

Penetration  Theory 


Incompressible  Bernoulli 

Jet  pressure 


p2,=Y,+-P„(v-uy 


Target  pressure 


pit  ~Yt  +~PotU' 


PD  (hydrodynamic) 


pd= 

l  i  p, 


PD  (target  strength 
important) 


PD 

UPj  2(7 

L  ] 

1  p,  p,  (V-U)2 

<J=Yt-Yj 
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Compressible  Jet 

Penetration  Theory 

•  Multiple  discrete  impacts;  particles  “pile 
up”? 

•  Macroscopically  -  jet  compresses 

•  Incompressible  Bernoulli  not  applicable 
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v-u 
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Compressible  Jet 

Penetration  Theory 


Jet  Material  Compaction  Curve 


•  Flis  &  Crilly  (18th  ISB):  compressible 
target 

•  Reverse  their  analysis  here 


p0  =  initial  (distended)  jet  density 
ps  =  solid  (pore-free)  jet  density 
<|>  =  1  -  p0  /  p  s  =  jet  porosity 
R  =  jet  compaction  initial  resistance 

sscppaipcbieb 


Compressible  Jet 

Penetration  Theory 


•  Distended  jet  at  p0  traveling  at  (V-U) 

•  Compacted  to  ps ,  decelerated  to  w-,: 


mass 

momentum 


^ . 

II 

1 

P()j+  Pa(w<ij)2  - 

zPii+  P\(wijf 

full  compaction 


P\j  ~  Pi  i  -  Pc 
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Apply  incompressible  Bernoulli  to  compacted 
jet  impact  gcpinupdilM. 


Compressible  Jet 

Penetration  Theory 


Jet  pressure 
Target  pressure 

PD  (hydrodynamic) 


P2i=R+\po(v-u?Q^i 


Compressibilit 


PD 

L 


PD  (target  strength 
important) 
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Compressible  vs.  Incompressible 

Jet 

Penetration  Theory 

Incompressibl 

e 

-  Compressible 

Jet 

| ,•  =Yj+±p0(V-U)2 

^p2J=r+^Po(v-u)2(  1+0 

Pressure 

Hydrodynamic 

Penetration 

1  PD  Upj  1 

ii,  S 

1™ 

—(!+<*) 

Pt 

Depth 

/i  — > 

(i+0 
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Compressible  vs.  Incompressible 

Jet 

Penetration  Theory 

•  (1+<)))  A, 

•  Porous  jet  model  reduces  to 
incompressible  model  at  limits 

•  At  upper  limit  (<f)=1 ;  A=2),  p  zero 

•  For  a  given  length  and  macroscopic 
density,  porosity  increases  the  following: 

-  Penetration  depth 

-  Impact  pressure 
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Jet  Density,  Length,  Porosity 

Summary 

Consider  3  cases:  (Fixed  mass,  velocity, 
diameter) 

pL=constant 

Case  A:  ps,  incompressible;  mmmm  - ►  v 

length=Ls 

PSI  length=L0 

Case  C:  p0,  incompressible;  MMMWillllM  - ►  V 

length=L0 
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Jet  Density,  Length,  Porosity 

Summary 

Hydrodynamic  Penetration  Depth 


Case  A: 


Case  B: 


Case  C: 


PD  b  —  Lq  ^ 

—(i+0 

Pt 

=pda 

+  l 

(alternati/eform) 

~PDC 

■\l  ]-+</> 

(alternati/eform) 

PDC=L  & 
V  Pt 
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Jet  Density,  Length,  Porosity 

Summary 

Hydrodynamic  Penetration  Depth 


pdb>pdc>pda 


2Cflnupeb06k 
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Jet  Density,  Length,  Porosity 

Summary 


Case  B: 


Case  C: 


p»=\p.<y-uy<\-4>2) 

=—pJV-U )2(\+<fr)  (alternatveform) 

Pc=\p«<y-u? 
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Jet  Density,  Length,  Porosity 

Summary 


Dynamic  Jet  Pressure 


P  >P  >P 

1  A^1  B^1  C 


Jet  Density,  Length,  Porosity 

Summary 

Consider  3  cases:  (Fixed  mass,  velocity, 
diameter) 

pL=constant 

Case  A:  ps,  incompressible;  mmmm  - ►  v 

length=Ls 

PSI  length=L0 

Case  C:  p0,  incompressible;  MMMWillllM  - ►  V 

length=L0 
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Jet  Density,  Length,  Porosity 

Summary 


•Greatest  impact  pressure 
•Shallowest  hydrodynamic  penetration 


Case  A 


•Deepest  hydrodynamic  penetration 
•2nd  greatest  impact  pressure 


Case  B 


Case  C 
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Impact  Pressure  &  Target  Strength 


What  about  target  strength? 

Consider  some  critical  point  of  interest  ( Jcrit,  Pr  rit) 


Characteristic 
target  property 


P 


crit 


isobars  in  (pj,  V)  coordinates 


Incompressible;  1 <X<2 


=  p,uL  +2(j 

W-ucrit)2 


Compressible 
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There  is 
no 

unique 

V min 
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Impact  Pressure  &  Target  Strength 


What  do  these  isobars  look  like? 


Target 

Material 

Pt 

(g/cc) 

Yt  (kbar) 

^crit 

Ps 

(g/cc) 

R=Yj  (kbar) 

Mild  steel 

7.86 

3  (UTS) 

2Yt 

16 

0 

Aluminium 

2.7 

3  (UTS) 

it 

it 

Concrete 

2.25 

0.3  (UCS) 

(i 

it 

it 

Sandstone 

2.25 

0.6  (UCS) 

(i 

it 

n 
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Impact  Pressure  &  Target  Strength 

Isobars  in  (pj,V)  space 


Incompressible 

Bernoulli 

Continuous 
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Impact  Pressure  &  Target  Strength 

Isobars  in  (pj,V)  space 


Incompressible 

Bernoulli 

Particulated 
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Impact  Pressure  &  Target  Strength 

Isobars  in  (pj,V)  space 


Mild  Steel 


Jet  Velocity  V  (cm/s) 


Compressible  (porous) 
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Jet  Density  rho  (g/cc  I  Jet  Density  rho  (g/cc 


Impact  Pressure  &  Target  Strength 

Isobars  in  (pj,V)  space 


2.0E+04  4.0E+04  6.0E+04  8.0E+04  1.0E+05  1.2E+05 


Jet  Velocity  V  (cm/s) 

Sandstone 


Concrete 


2.0E+04  4.0E+04  6.0E+04  8.0E+04  1.0E+05  1.2E+05 

Jet  Velocity  V  (cm/s) 
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Jet  Density  rho  (g/cc' 


Impact  Pressure  &  Target  Strength 

Isobars  in  (pj,V)  space 
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Conclusions  (1) 

•  Developed  treatment  of  compressible  jet 

-  Reduces  to  well-known  expressions  for  solid  & 
fully-particulated  jets 

-  For  a  given  mass,  velocity,  diameter: 

•  Porous  (compressible)  penetrator  penetrates 
deeper  than  incompressible  penetrator  of  same  L, 

P 

•  . .  .also,  deeper  than  shorter  penetrator  of  higher  p 
...produced  impact  pressure  which  is  intermediate 
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Conclusions  (2) 

•  Looked  at  steady-state  impact  pressure 

-  Presented  isobars  in  (pj,V)  coordinates 

•  Compressible  jet  model  interpolates  between  solid 
&  fully-particulated  incompressible  jet  model 

-  Highly-distended,  low-velocity  jets  may  effectively 
penetrate  moderate-strength  geologic  targets 


This  approach,  so  far,  neglects  transients 
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The  residual  damage  in  CFRP 
composite  after  ballistic  impacts 
(experiments  &  simulations) 


Survivability  of  aircraft 


TNO  I  Knowledge  for  business 


TMO 


Contents 


•  Introduction 

•  Ballistic  test  data 

•  Residual  velocity 

•  Internal  damage 

Material  model  of  CFRP,  AS4/3501 

•  Simulation  results: 

•  Residual  velocity 

•  Internal  damage 
Conclusions  and  future  work 


Koen  Herlaar,  M.Sc.  koen.herlaar@tno.nl 


Introduction 


Survivability  of  aircraft: 

Combined  threats  :  both  blast  and  projectile  impacts 
Lighter  platforms:  more  composite  material 


c 


Combined  threats  +  new  materials 

=>  New  failure  models 
New  survivability  tools 


Koen  Herlaar,  M.Sc.  koen.herlaar@tno.nl 
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C-scan  result 


m  m  M  M 

S9 

SHE  K  ^  "  M 

FSP,  Fragment  Simula 

iting  Projectile 

D1=1660mm2(20dB) 


D2=1350mm2(20dB) 


D6=1183min2(20dB) 


D3=  1 605min2(20dB) 


D4=1263inm2(20dB) 


D5=  1 740mm2(20dB) 


Koen  Herlaar,  M.Sc.  koen.herlaar@tno.nl 


Vancouver,  November  18th,  2005  THt 


Ballistic  tests 


Residual  velocity: 

Threat :  1.1  gram  FSP 
(Fragment  Simulating  Projectile) 


5 


Internal  damage 
14.3  mm  (C-Scan): 


Koen  Herlaar,  M.Sc.  koen.herlaar@tno.nl 


3=  1 605mm2(20dB) 


D4=  1 263mm2(20(1B) 


Vimpact  597  m/s 

R damage  26-28  mm 


V impact  698  m/s 
R damage  26-27  mm 


V impact  822  m/s 
R damage  22-24  mm 


Vancouver,  November  18th,  2005  TKI 


Contents 


•  Introduction 
Ballistic  test  data 

Residual  velocity 
Internal  damage 

•  Material  model 

•  Simulation  results: 

•  Residual  velocity 

•  Internal  damage 

Conclusions  and  future  work 


Koen  Herlaar,  M.Sc.  koen.herlaar@tno.nl 
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Material  model  of  CFRP,  AS4/3501 


Advanced  Damage  Material  Model  (ADAMMO)  within  AUTODYN®  : 

•  An  orthotropic  elastic  model  : 

Parameters  from  quasi-static  material  tests  (tension  and  V-notch) 
A  linear  EOS  (assumed  linear;  inverse  flyer  plate  necessary) 

•  An  orthotropic  damage  model  with 

Orthotropic  failure  criteria  (tension  and  V-notch) 

An  orthotropic  softening  algorithm  (data  from  literature) 
Orthotropic  post  failure  response;  tensile  stresses  are  still 
allowed  in  non-failed  material  directions 


Koen  Herlaar,  M.Sc.  koen.herlaar@tno.nl 
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Residual  velocity  (my 


Simulation  results 

Residual  velocity  (threat  1.1  gram  fsp): 
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Simulation  results 


v 


Impact 


=  597  m/s 


Delamination 


AUT0DYN-3D  v5.0  from  Century  Dynamics 


Cycle  1122 

Time  1 .000E-003  ms 

Units  mm,  mg,  ms 

1.1G  FSP  IMPACTING  ON  14.3  MM  CFRP 
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=  14.3  mm 


Fiber  failure 


AUT0DYN-3D  v5.0  from  Century  Dynamics 


Cycle  1122 

Time  1 .000E-003  ms 

Units  mm,  mg,  ms 

1.1G  FSP  IMPACTING  ON  14.3  MM  CFRP 


Vancouver,  November  18th,  2005 


Simulation  results 


v 


Impact 


=  597  m/s 


Delamination 


AUT0DYN-3D  v5.0  from  Century  Dynamics 


Cycle  3205 
Time  3.000E-003  ms 
Units  mm,  mg,  ms 

1.1G  FSP  IMPACTING  ON  14.3  MM  CFRP 
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=  14.3  mm 


Fiber  failure 


AUT0DYN-3D  v5.0  from  Century  Dynamics 


Cycle  3205 
Time  3.000E-003  ms 
Units  mm,  mg,  ms 

1.1G  FSP  IMPACTING  ON  14.3  MM  CFRP 


Vancouver,  November  18th,  2005 


Simulation  results 

V 'lmpact  =  597  m/s,  thickness  =  14.3  mm 


Delamination 


Cycle  6722 
Time  5.000E-003  ms 
Units  mm,  mg,  ms 

1 .1 G  FSP  IMPACTING  ON  14  3  MM  CFRP 


Fiber  failure 


Cycle  6722 
Time  5.000E-003  ms 
Units  mm,  mg,  ms 

1 .1 G  FSP  IMPACTING  ON  14.3  MM  CFRP 


11 


Koen  Herlaar,  M.Sc.  koen.herlaar@tno.nl 


Vancouver,  November  18th,  2005  THt 


Simulation  results 


V impact  =  597  m/s,  thickness  =  14.3  mm 


Delamination 


Cycle  12862 
Time  7.000E-003  ms 
Units  mm,  mg,  ms 

1 .1 G  FSP  IMPACTING  ON  14  3  MM  CFRP 


Fiber  failure 


Cycle  12862 
Time  7.000E-003  ms 
Units  mm,  mg,  ms 

1 .1 G  FSP  IMPACTING  ON  14.3  MM  CFRP 
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Simulation  results 


V impact  =  597  m/s,  thickness  =  14.3  mm 


Delamination 


Cycle  21136 
Time  9.000E-003  ms 
Units  mm,  mg,  ms 

1 .1 G  FSP  IMPACTING  ON  14  3  MM  CFRP 


Fiber  failure 


Cycle  21136 
Time  9.000E-003  ms 
Units  mm,  mg,  ms 

1 .1 G  FSP  IMPACTING  ON  14.3  MM  CFRP 
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Simulation  results 


v 


Impact 


=  597  m/s 


Delamination 


AUT0DYN-3D  v5.0  from  Century  Dynamics 


Cycle  31776 
Time  1.100E-002  ms 
Units  mm,  mg,  ms 

1.1G  FSP  IMPACTING  ON  14  3  MM  CFRP 
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=  14.3  mm 


Fiber  failure 
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Simulation  results 


v 


Impact 


=  597  m/s 


Delamination 


AUT0DYN-3D  v5.0  from  Century  Dynamics 


Cycle  42598 
Time  1 .300E-002  ms 
Units  mm,  mg,  ms 

1.1G  FSP  IMPACTING  ON  14  3  MM  CFRP 
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;kness  =  14.3  mm 


Fiber  failure 


Vancouver,  November  18th,  2005  Til 


Simulation  results,  total  internal  damage 


V impact  =  597  m/s,  thickness  =  14.3  mm 


*~L 

impact-3-dikker-v597 
Cycle  0 

Time  0.000E+000  ms 
Units  mm,  mg,  ms 


1.1G  FSP  IMPACTING  ON  14.3  MM  CFRP 


Koen  Herlaar,  M.Sc.  koen.herlaar@tno.nl 


in 


16 


Vancouver,  November  18th,  2005 


Total  internal  damage  (comparison) 


V, 


impact 


=  597  m/s 


V: 


impact 


=  698  m/s 


V: 


impact 


=  822  m/s 
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Total  internal  damage  (comparison) 

•  Quantitative  comparison  is  very  difficult  ->  No  direct  link 
Experiments:  Simulations: 


C-scan  with  20  dB  threshold 

Damage  radius  «  22  -  27  mm 

Little  difference  in  damaged  area 
within  velocity  range 

However  small  trend: 

vimpact  increases  => 
damaged  area  decreases 


“C-scan”  with  60  %  threshold 

Damage  radius  «  20  -  25  mm 

Little  difference  in  damaged  area 
within  velocity  range 

However  small  trend: 

vimpact  increases  => 

damaged  area  increases 


Improvements  can  be  made  by  measuring  the  softening  behaviour 
in  thickness  direction  and  by  performing  inverse  flyer  plate  tests 

Vancouver,  November  18th,  2005  THt 
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Conclusions  and  future  research 


With  limited  material  tests  ->  good  material  behaviour 
Simulation  results  are  consistent  with  experimental  residual  velocity 
Internal  damage  is  difficult  to  compare  with  experiments 
The  damage  observed  is  in  the  same  order  of  magnitude,  however 
the  trend  is  inconsistent  with  the  experiments. 

Future  research: 

Further  improvement  of  material  model  (softening,  flyer  plate) 
Combine  blast  loading  with  fragment  impacts  on  construction 
Translate  damaged  area  into  strength  reduction  of  construction 
With  the  knowledge  from  the  FEM  models  -> 
improve  the  current  vulnerability  and  survivability  tools 
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Why  Soft-Recovery? 


Munitions  Director 


•  Improved  design  of  EFPs  (Explosively  Formed  Projectiles)  using  state-of- 
the-art  constitutive  descriptions  require  ab  initio  information  of  the  high- 
rate,  multi-axis  stress  deformation  path 

•  Explosive  shock  effects 

•  Predominate  crystallographic  orientation  (i.e.  texture)  and  its  dynamic 
evolution 

•  Classical  flash  radiography  and  high-speed  photography  only  captures 
geometry  information 

•  Collaboration  with  Dr.  Paul  Maudlin,  Los  Alamos  National  Lab.,  TCG-I 


Flash  Radiograph  of 
an  EFP  in  flight 


Recovered 

EFP 


Cross-section  of  Code  simulation,  EPIC 
recovered  EFP  Anisotropic  MTS 


Example  EFP  Formation 
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Developmental  Engineering  Tools 
3D,  Anisotropic  Material  Descriptions 


Munitions  Directorate 


Material  Processing 

•  Conventional 

•  Unconventional 


Physical 

Characterization 

Grain  size 
Chemistry 
Crystallography 


Yield  Surface  Projections 


in 


Initial 


Recovered  Projectile 


Calculations  by  Dr.  Paul  Maudlin,  LANL,  under  TCG-I 
“Computational  Mechanics  and  Material  Models”,  EPIC-3D 


Mechanical 

Characterization 

Low  Rate  &  High  Rate 
Temperature 
Uniaxial 
Taylor  Impact 


Explosive  Metal 
Forming 
Experiments 

Shock  loading 
High-rate 
Tri-axial  stress 


The  Approach  Taken  Here 


*  Apply  a  simplistic,  but  adequate  equation-of-motion  for  the 
projectile  through  the  media 

*  Arrange  mathematical  relationships  that  relate  to  the  physical 
experiment  and  instrumentation  capability 

*  Conduct  experiments  from  which  the  data  is  used  to  calibrate  the 
soft-recovery  media 

*  Compare  results  of  the  general  model  with  specific  experiments 

*  Use  the  model  and  calibrated  media  constants  for  tailored  design 
and  construct  of  soft-recovery  apparatus 
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Mathematical  Model 


Region  1:  Drag  Force  Model 

m v  =  -  —  CDApv2 ,  v  >  vc 

Region  2:  Poncelet  Form 

mv  =  -A(f3v2  +  i?)  ,  v  <  v( 


m  is  mass  of  the  projectile 
v  is  current  projectile  velocity, 

A  is  projectile  cross-section  area 
p  is  the  density  of  the  target  medium 
CD  is  a  dimensionless  drag  coefficient 
p  is  a  coefficient  (dim.  of  density) 

R  is  a  target  strength  factor  (dim.  of  stress) 


Allen,  etal.  J.  Applied  Physics,  1957 
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Mathematical  Model  (cont’d) 


Munitions  Director 


Assumptions  made  on  the  solution: 

•  Model  will  be  applied  for  soft-catch  sections  where  v  >  vc 

•  Drag  coefficient  is  not  dependent  on  velocity 

Direct  integration  of  the  Drag-Force  model  gives: 


and 


2  m 


where: 

•  z  is  the  displacement  into  that  section  of  the  soft  catch 

•  t  is  the  experimentally  obtained  time  at  the  displacement  z 

•  v0  is  the  entrance  velocity  to  that  section 

•  v  is  the  exit  velocity  from  that  section 
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Mathematical  Model  (cont’d) 


•  By  dividing  the  prior  two  equations,  a  single  relationship  is 
obtained,  dependent  only  on  parameters  that  can  be  directly 
obtained  during  an  experiment. 


z 

t 


f 


V, 


°— 1 

V  J 


^  (v  ' 

=  v0£n 


o 


where: 

•  z  is  the  displacement  into  that  section  of  the  soft  catch 

•  t  is  the  experimentally  obtained  time  at  the  displacement  z 

•  v0  is  the  entrance  velocity  to  that  section 

•  v  is  the  exit  velocity  from  that  section 
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Solution  Steps 


Step  1 : 


tyv  j 

Step  2: 

1  CDAp_ 

2  m 

Step  3: 

1  CDApt 

2  m 


Munitions  Directorate 


v0£n 


fv  } 
J_o_ 

V  J 


Is  applied  to  each  catch  section  with  experimental 
time-position  data,  using  Newton’s  method  for 
convergence  on  an  iterative  root  solution  to  get  the 
ratio  v/v,  where  v0  is  the  exit  velocity  from  the 
previous  section  and  v  is  our  unknown 


o 


V 


The  ratio  v/v  is  then  put  back  into  this  equation  to 
obtain  a  CD  for  each  penetrator/soft-catch  material 
combination 


1_A 

v  v0 


The  CD  for  sections  of  each  material  are  averaged, 
and  then  used  to  obtain  an  estimated  exit  time,  t,  for 
comparison  to  the  experimental  values 
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The  Soft-Catch  Process 


A  gradient  in  the  media  density  safely  decelerates  the  projectile 

Instrumented  shotline  captures  the  time,  position,  and  velocity  data  to 
feed  the  transcendental  relationship 

The  section  of  water  not  only  provides  a  gradual  increase  in  media 
density,  but  also  serves  to  quench  the  projectile 

Dual,  orthogonal  radiographs  were  used  in  pre-impact,  free-flight  to 
establish  external  geometry  and  entrance  velocity 


i2ineh  cross-section  Decelerates  to  rest  and 


12-ft  Polystyrene 
(32kg/m3) 


8-ft  Vermiculite  8-ft  Fiberboard  8-ft  Water  Sand 
(126kg/m3)  (256kg/m3)  (1000kg/m3) 


Experimental  Setup 


Munitions  Director 


1 .  Three  warhead  liner  types  were  shot 

•  3  shots  with  “Tantalum  Design  1” 

•  2  shots  with  “Tantalum  Design  2” 

•  1  shot  with  “Cu  EFP” 

2.  Fine  Grain  Octol  explosive  (65%HMX,  35%  TNT) 

3.  Design  was  for  a  simple  'fold-over’  projectile 

•  Explosive  shock  conditions 

•  Representative  strain-rate  and  strain  paths 


Rear  RP-1 
Initiator 


Ex[ 


1 


Steel  Case 


Warhead 

Liner 
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Results  from  Applying  the  Model 


TABLE  I.  DRAG  COEFFICIENTS,  CD,  FOR  ALL  EFP  TYPES 


Polystyrene 
(p  =  32.0  kg/m3) 

Vermiculite® 

(p=  126.4  kg/m3) 

Fiberboard 
(p  =  256  kg/m3) 

Ta  Design  1 
(avg.  of  3  shots) 

0.777 

1.534 

0.395 

Ta  Design  2 
(avg.  of  2  shots) 

0.84 

0.88 

0.86 

Cu  EFP 
(1  shot) 

0.77 

0.94 

0.76 
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Comparison  with  the  Theory 


TABLE  II.  MODEL  COMPARISON  WITH  “TA  DESIGN  2,  SHOT  1” 


Tantalum  Design  2,  Shot  1 

Impact  Velocity  =  1440  m/s 

Experiment 

est.  Exit 

Experiment 

est.  Exit  time 

DrfferenceLin 

Media 

Velocity  (m/s) 

Velocity  (m/s) 

Exit  time  (us) 

(us) 

'  time  (%) 

1395 

1381 

437 

432 

-1.1% 

1382 

1324 

878 

883 

0.6% 

Polystyrene 

1273 

1270 

1357 

1353 

-0.3%  \ 

1229 

1218 

1853 

1844 

-0.5% 

1212 

1168 

2356 

2355 

-0.1% 

1146 

1120 

2888 

2888 

0.0% 

1039 

941 

3475 

3482 

0.2% 

872 

791 

4174 

4190 

0.4% 

Vermiculite 

721 

665 

5019 

5032 

0.3% 

615 

559 

6010 

6033 

0.4% 

512 

469 

7201 

7225 

0.3% 

419 

395 

8656 

8644 

-0.1% 

325 

280 

10529 

10488 

\  -0.4% 

Cellotex 

235 

198 

13119 

13092 

\-0.2%/ 

•  Excellent  results  validate  the  assumption  of  simple  drag-force 
relationship,  v>  vc 

•  Note,  penetration  velocities  are  below  that  reported  in  Mayfield, 
et  al  for  vc 
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Recovered  Projectiles 


Shot 

Impact 

Velocity 

(m/s) 

Recovered 
Mass  (kg) 

(%initial) 

Location 

Recovered 

Cross- 

section 

Picture 

Ta  Design  1 
Shot  1 

1423 

0.726 

(83.4%) 

1 2-in  into 
sand 

T a  Design  1 
Shot  2 

1397 

0.720 

(82%) 

24-in  into 
sand 

+ 

Ta  Design  1 
Shot  3 

1389 

Not 

Available 

1 8-in  into 
water 

Not 

available 

Ta  Design  2 
Shot  1 

1440 

0.738 

(92.6%) 

41 -in  into 
fiberboard 

# 

;  -  \ 

T a  Design  2 
Shot  2 

1422 

0.761 

(93.4%) 

60-in  into 
fiberboard 

0 

Cu  EFP 

2030 

0.409 

(90.8%) 

fiberboard 

section 

• 
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Discussions/Conclusions 


V 


Areas  for  further  investigation: 

1 .  Determine  if  the  tacit  assumption  that  v  >  vc  used  is  valid 
—  Find  vc  in  these  media 

2.  Determine  if  CD  is  independent  of  velocity  as  assumed 
(velocities  were  in  a  relatively  narrow  band) 

—  Find  CD  for  a  wider  range  of  velocities  in  each  media  and 
compare 

3.  Use  constants  for  predictive  design  of  soft-catch  build  up 
and  capture  higher  velocity,  more  tactical  (collapsed) 
projectiles 
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Discussions/Conclusions 


V 

42# 


Munitions  Directorate 


An  alternative  approach  to  calibrating  projectile/media  constants: 

1 .  Construct  a  soft  catch  with  just  one  catch  material  at  a  time  to 
obtain  CD  for  all  velocities  above  and  below  vc 

2.  Place  velocity  screens  at  closer  intervals  for  better  resolution 
such  that  vc  might  be  obtained 

Total  penetration  depth  z  =  S  is  found  when  v  =  0,  i.e.  projectile  comes 
to  rest 

The  original  equations  of  motion  can  be  arranged  to  solve  for  total 
penetration  into  the  single  media  and  extract  vc 


S  =  ^~  In 


o 


+ 


m 


\VcJ 


2J3A 


•In 


f  n  2  A 

pvc 

^-  +  1 
R 


CDAp 


v 


J 


Summary 


V 

4iP 


*  A  simple  theory  has  been  successfully  applied  for  the 
design  of  a  soft-catch  apparatus 

*  Model,  thus  far,  yields  excellent  agreement  with 
experimental  time  data 

*  Further  experiments  needed  to  support  underlying 
assumptions 

*  These  interest  items  will  be  the  subject  of  future 
experiments 
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Motivation 

An  unexplained  non-monoticity  of  the  scaled  blast  duration  is 
evident  in  CONWEP  data. 

A  possible  explanation  might  be  related  to  non-ideality  of  the 
explosion  (post-detonation  burning). 


Hopkinson  Scaled  Distance  [ft/lb1/3] 
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Non-Ideal  Blast  Waves 


I  D  F 


Ideal  Explosion 

Non  Ideal  Explosion 

All  chemical  energy  is 
extracted  in  the  detonation 
front 

Some  energy  remain  due  to 
negative  oxygen  balance 

Gas  diffusion  and  surface 
turbulence  are  neglected. 

Gas  diffusion  and  surface 
turbulence  taken  into  account 

Detonation  products  do  not 
bum 

Detonation  products  bum 
when  oxygen  is  supplied 

22nd  International  Ballistics  Symposium,  2005 


The  Goal  of  our  Study 


To  model  effectively  the  non  ideal 
blast  effects  (i.e.  gas  diffusion  and 
burning)  in  order  to  understand 
their  influence  on  the  blast 

profile. 


time 
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Spherical  Explosion  Shock  Wave  Dynamics 


C  -  Det.  products  -  air  surface. 

51  -  Primary  shock  wave 

R  -  Rarefaction  shock  wave. 

52  -  Secondary  shock  wave. 


©  D.W.  Boyer,  J.  Fluid  Mech.,  9(3)  1960 
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Gas  Diffusion 

Classically,  the  diffusion  of  two  gases  controlled 
by  three  gradients: 

-  Concentration 

-  Pressure 

-  Temperature 


Concentration  estimated  diffusion  velocity: 
500-1500  m/s 
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Model  Hypothesis 


All  gases  are  ideal  gases  with  various  adiabatic  constants. 

Adiabatic  constant  of  mixed  gas  is  a  concentration 
weighted  average. 

The  burning  process  affects  only  the  internal  energy  and 
the  adiabatic  constant  of  the  gas. 


Relative  concentrations  inside  detonation  products  cloud: 


•  r;  -  air 

•  £  -  Pre-burned  gas 


A 


rj  +  %  +  /3  =  1 


•  P  -  burning  products 


c£!ms£] 
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Energy  release  dynamics 


I  D  F 


Initial  conditions  (det.  products  cloud):  rj  —  0;  ^  —  1;  [5  —  0 


Diffusion: 


r/(r,t)  =  7j0  exp 


r 


V  ut  J 


Rate  of  change  in  internal  energy:  Eb  =  min  (y  77  ,%)-E 


Burning  Products  concentrations: 


m= 


Eb{t) 

p  final 

^ b 


Pre-bumed  gas  concentration: 


l=i -P-n 


Eu  =  min 


f 


W 


7  rj(r9t)9  1  -rj(r9t)- 


V 


V 


final 

^ b  j  j 


.  p0 
^ b 
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Numerical  Simulation 


I  D  F 


Simulation  conducted  using  AUTODYN  ver.  6.0. 

New  EOS  (based  on  ideal  gas)  was  defined  in  order 
to  implement  the  model. 

Two  stages  simulation: 

-  Detonation 

-  Blast  propagation  &  burning 

The  simulation  set-up  was  a  spherically  symmetrical 
explosion  of  5kg  TNT  charge. 
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Numerical  Simulation 

ID  axially  symmetrical  explosion 


t  [ms] 
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Simulation  results 


I  D  F 


Without  burning  With  Burning 
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Simulation  Results 


I  D  F 


(i)  (ii) 

r—0.8  m  r=1.6m 


Without  Burning 
With  Burning 
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Experimental  Set  Up 


I  D  F 


•  Two  energy  equivalent  free 
field  explosion  tests  conducted: 

-  5  kg  TNT 

-  4.2  kg  C-4 


c£lfe£] 
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Detonation  Products  of  TNT  and  C-4 


I  D  F 


TNT 

C-4 

[%mole] 

[%mole] 

H  2 

21.4 

12.6 

CO 

54.5 

22.0 

co2 

5.5 

10.5 

h2o 

1.5 

18.0 

n2 

12.5 

28.0 

Others 

5.5 

9.9 

co+h2 

76% 

35% 
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Experimental  Results 


I  D  F 


Pulse  durations 


Pulse  pressure  peak 
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Hopkinson  Scaled  Distance  [ft/lb1/3]  Hopkinson  Scaled  Distance  [ft/lb1/3] 
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Model  Calibration 


Calibration  parameters: 

u,  TJq, 

The  model  was  calibrated  employing: 

-  CONWEP  data 

-  Experimental  results 
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Results 
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Summary 


An  effective  model  is  proposed,  capable  of 
reproducing  the  burning  effect  on  a  non-ideal  blast 
waves. 

Good  agreement  obtained  between  numerical  and 
experimental  results. 

Post  detonation  burning  affects  the  blast  profile  in 
the  near  field  (scaled  distance  <  4  ft/lb1/3). 

When  non  ideal  blast  effects  are  important  the  TNT 
equivalence  convention  must  be  reconsidered. 
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ABOUT  VARYING  SHAPED  CHARGE  LINER  THICKNESS 


Pierre  Y.  Chanteret,  Andre  Lichtenberger 

ISL,  French-German  Research  Institute  of  Saint-Louis, 
BP  34  F68301  Saint-Louis,  France 


The  influence  of  linear  wall  thickness  on  shaped  charged  jet  behavior  is  investigated. 
X-ray  pictures  are  presented  of  jets  from  60°  copper  liners  with  thickness  ranging  from 
0.4%  to  7%  of  the  charge  diameter.  First,  it  is  shown  that  under  a  minimal  thickness, 
about  0.25  mm,  no  compact  copper  jet  is  formed.  Then,  the  influence  of  liner  wall 
thickness  on  jet  fragmentation  is  studied  and  a  variation  of  the  characteristic  velocity 
difference  between  consecutive  fragments  is  reported.  Introducing  this  variation  in  a  1 D- 
code  allows  for  reproducing  the  jet  fragmentation  behavior  for  the  whole  range  of  liner 
thickness  considered. 


#  =  0.4  riMTT 


4  3  0.6  mm 


4  a  1.0  mm 


o  a  1 .9  mm 


«  a  3*0  mm 
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0.20  mm 


7.7  hm/a 


7.2  krn/t 


6.0  km/s 


5.4  km/D 


4.2  knV* 
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Experiments  in  plate  cutting  by  shaped 
high  explosive  charges 

By  A.I.O.  Zaid,  J.B.  Hawkyard  and  W.  Johnson 
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Charge  configuration,  using  high  quality  explosive  PBX  -  LX07, 
based  on  90%  HMX  and  10%  Viton. 
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Same  charge  configuration  with  different  explosives  PEX01,  based  on  93  wt%  HMX 
and  siliconic  binder  and  formed  by  injection  molding  into  a  cylindrical  shape. 

After  polymerization  (curing)  the  hollow  cone  is  machined  to  its  final  dimension. 
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Same  charge,  same  explosive  PEX01 ,  but  this  time 
the  charge  is  formed  by  direct  injection  of  the 
explosive  on  the  liner. 
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WATER 


OIL 
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WATER 


Simple  demonstration  of  Rayleigh  Taylor 
Instability 
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Demonstration  of  the  relevance  of  the  RTI 
to  shaped  charge  jet  formation. 


Initial  grid  setup  with  a  sinusoidal  liner. 
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The  growth  of  the  initial 
grid  during  collapse  with 
a  liner  made  of  copper 
(nominal  strength). 


No  change  in  shape 
during  collapse  if  the 
simulation  is  performed 
with  increased  strength 
by  a  factor  of  4. 


The  same  initial  charge 
configuration  except  for  the  liner. 
This  is  the  simulation  results  with 
aluminum  liner  (nominal  strength). 
Actually,  the  interface  between  the 
liner  and  the  explosive  tends  to  be 
smoothed  to  a  straight  line. 
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Repeating  the  experiment  with  the  machined  cone  in 
the  PEX01  explosives,  but  with  Molybdenum  liner. 
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Repeating  the  experiment  with  the  machined  cone  in 
the  PEX01  explosives  with  2%  copper  liner. 
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The  crater  opening  on  aluminum  target  by  a  configuration  with  3%  CW  thickness  copper  liner. 


Same  as  above,  but  with  a  variable  thickness  copper  liner  from  1%  CW  near  the  apex  to  3%  CW  near  the  base. 
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summary 
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The  RTI  might  have  a  destructive  role  in  jet  formation 
unless  the  charge  configuration  is  designed  beyond 
certain  threshold  limits: 

Minimum  liner  thickness  which  depends  on  liner 
density  and  strength. 

Explosive  quality,  especially  near  the  contact 
surface  with  the  liner. 

Linear  charges  experiments  have  proven  to  be 
much  more  sensitive  to  initiate  instabilities.  Some 
reasonable  assumptions  have  been  offered  to 
explain  the  difference  between  conical  and  linear 
charges.  Further  investigation  is  needed  in  order 
to  gain  a  conclusive  understanding. 
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The  charged  configuration  with  PBX-LX07  scaled 
down  by  factor  6.  Liner  thickness  around  0.08  mm. 
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Repeating  the  experiment  with  the  machined  cone  in 
the  PEX01  explosives,  but  with  a  thin  layer  of  lexan  (2 
mm  thick)  inserted  between  the  liner  and  the 
explosive. 
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Background 


•  Rocket  based  lift  to  orbit  costs  between  $4000  to  $14,000  per  pound 
depending  on  launcher  and  payload 

•  For  the  Mars  mission,  NASA  anticipates  a  need  for  large  quantities  of 
supplies  and  material  delivered  to  earth  orbit 

•  NASA  is  interested  in  the  possibility  of  earth  launch  to  space  of 
acceleration  tolerant  payloads 

o  Low  Earth  Orbit  (LEO) 
o  500  kg  payload 
o  1000-2000  kg  total  launch  mass 

•  ARL  asked  to  lead  examination  of  selected  launcher  technologies: 

o  Slingatron 

o  Blast  Wave  Accelerator 
o  EM  Coil  Gun 
o  EM  Rail  Gun 
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Launch  Calculations 


Consider  orbit  of  International  Space  Station  (ISS): 

o  359  km  at  51.4°  inclination  to  the  equator 
o  Orbital  velocity  7.7  km/s 

Flight  Body: 

o  10°  half  angle  blunted  conical  nose 
o  Base  diameter,  D  =  0.46  m 
o  Cylindrical  section,  L  =  2.5  D 
o  Nose  radius,  rn  =  0.0575  m 
o  Total  flight  mass,  m  =  1000  kg 


Launch  Calculations 
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Assume: 

o  Projectile  just  reaches  apogee  at  ISS  orbit  Vm 

o  Newtonian  flow  theory  applies  =>  CD  =  0.119 

◦  Use  Rice,  et  al  [5]  to  estimate  velocity  at  exit  from  sensible 
atmosphere: 

Ve/Vm  =  exp[-C0po  A/2mpsin6] 
where  p  is  the  atmospheric  lapse  rate  =  1.1  x  104/m 
o  50%  structural  mass,  85%  rocket  mass  fraction,  I  =  305  s 


Compute  muzzle  velocity,  insertion  velocity,  and  payload 
mass  versus  elevation  angle  for  1000  kg  launch  mass. 
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13.48 

7.88 

5.74 

4.52 

3.81 

3.35 
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2.93 

2.88 

AV  (km/s) 

1.03 

2.62 

3.89 

4.91 

5.66 

6.28 

6.80 

7.28 

7.70 

mpi  (kg) 

329 

156 

72.5 

26.3 

0.94 

Propellant  system  mass  equals  or  exceeds 
available  non-structural  mass 
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High  Altitude  Research  Project 


HARP: 

First  serious  attempt  at  gun 
launch  to  space. 

Double  length  16”  Naval 
gun. 

Achieved  altitude  of  180  km 
(Vm  =  2.1  km/s)  in  1966  firing 
at  Yuma  Proving  Ground. 

Plans  underway  for  rocket 
propulsion  for  orbital 
insertion. 

Project  cancelled  in  1967. 


Martlet  2:  launch  mass  up  to 
215  kg  capable  of  25  kGee. 
Cost:  $3000  each 
Launch  Interval:  1  per  hour 
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Sling  atron 


Graphics  courtesy  of  J.  Jones,  NASA 


Concept:  Derek  Tidman,  USA,  1995 
Table  top  demo  =>  100  m/s 
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Slingatron  Functioning 
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ARL  MATLAB  Simulation 


u  =  0.007 


0  20  40  60  80 


time  (s) 

ARL  simulation  agreed  with  Tidman’s  analysis.  Slingatron  can  in  theory 
accelerate  a  body  to  high  velocity,  but  there  are  significant  unknowns  - 
one  of  which  is  high  speed  friction. 
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Slingatron  Launcher 


•  For  1,000  kg  launch  mass  at  8,000  m/s,  a  Slingatron  solution 

Circuit  Diameter  =  300  m 
Tube  Diameter  =  0.4  m 
Tube  Wall  Thickness  up  to  0.25  m 
Gyrational  Frequency  =  9  cycles/sec 
Weight  =  20,000  tons 

•  Issues: 

Structural  integrity 
Friction  and  surface  wear 
In-Bore  stability 

Fabrication  of  curved  guide  tube 

Drive  and  support  system  design  integrity 

Synchronization  of  multiple  drives 

Power  to  bring  to  speed  and  maintain  during  firing 

Reacting  the  acceleration  and  centrifugal  recoil  load 
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Blast  Wave  Accelerator 


Graphics  courtesy  of  J.  Jones,  NASA 
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Blast  Wave  Accelerator  Concept 


Sequentially  detonated  explosive  rings 


ARL  calculations  in  agreement. 


Max  Pressure  on  Projectile,  MPa 


RDECOM 

’iL 


Surface 

Deformation 


ARL  CTH  Results 


Y  (m) 


I 

O 


I 

o 


o 
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o 

CO 


Internal 

Cracking 


Simulations  showed  structural  damage  to  mild  steel  projectile. 
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Blast  Wave  Launcher 


•  For  1,000  kg  launch  mass  at  8,000  m/s: 

Charge  mass  per  Stage  =  10  kg 
Number  of  Stages  =  2870 
Gun  Length  =  861  m 

•  Issues: 

Projectile  integrity 

Launcher  integrity  (repetitive  firings) 
Timing  of  sequential  detonation 
Explosive  detonation  uniformity 
In-bore  stability 

Detonation  of  30,000  kg  of  explosive 
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Electromagnetic  Guns 


Sandia  Electromagnetic  Coilgun 


Greenfarm  Electromagnetic  Railgun 


Coilgun  has  launched  230  g  to 
1000  m/s 


Railguns  have  launched  ~2  kg  to 
3000  m/s  (1  g  to  7000+  m/s) 
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EM  Coilgun 


A  magnetic  travelling  wave  is  created  by  synchronized  sequential  switching  of 
individual  coils. 

Strong  centering  force  “levitates”  projectile  for  minimal  wall  contact 


Graphics  courtesy  of  Sandia 


K#ESiC_4g 

'^mRDECOM 

1L 


Sandia  Coilgun  Concept 


Satellite 


Graphics  courtesy  of  Sandia 
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Coilgun  Launcher 


•  For  837  kg  launch  mass  at  7,000  m/s: 

Launcher  Length  =  400  m 
Tube  Diameter  =  1.0  m 

•  Issues: 

Structural  integrity  of  projectile  and  launcher 

Pulsed  power  supply 

Switching 

In-bore  stability  of  levitating  projectile 
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EM  Railgun 


Rail  Repulsion  Force 


Magnetic  Field  (B) 


a 


Armature 


Current,  J 


For  velocity  >  3  km/s  a  plasma 
armature  is  used. 


Firing  similar  to  coilgun  -  near 
equator,  at  high  altitude,  and 
from  an  evacuated  launch  tube. 


Lorentz  Force  (JxB) 


Graphics  courtesy  of  IAT 
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Railgun  Launcher 


•  For  1,250  kg  launch  mass  at  7,500  m/s: 

Launcher  Length  =  1600  m 
Tube  Diameter  =  1.1  m 

•  Issues: 

Structural  integrity  projectile  and  launcher 
Rail  life 

Pulsed  power  supply 
Switching 

Plasma  armature  performance 
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Summary 


Comparison  of  launchers: 


vm  (m/s) 

Mproj  (kg) 

Em  (GJ) 

1-tube 

Amax  (kG) 

Slingatron 

8000 

1000 

32 

D  =  300  m 

43 

BWA 

8000 

1000 

32 

861 

55 

Coil  Gun 

7000 

837 

21 

400 

7 

Railgun 

7500 

1250 

35 

1600 

2 

Conclusions: 

•  Achieving  an  8  km/s  muzzle  velocity  did  not  violate  any  laws 
of  physics 

•  All  had  serious  engineering  and  materials  issues 

•  Significant  research  is  required 

•  Facilitization  costs  would  be  high 

•  All  are  high  risk 


Caseless  Ammunition  and  Advances  in  the 

Characterization  of 
High  Ignition  Temperature  Propellant 

(HITP) 

Erin  K.  Hardmeyer 

Armaments  Engineering  and  Technology  Center 
US  Army  Armament  Research  Development  Engineering  Center 

Ben  Ashcroft 

Alliant  Technical  Systems 
Thiokol  Propulsion 

NDIA  22nd  International  Symposium  on  Ballistics 

November  14-18,  2005 
Vancouver  BC,  Canada 
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Outline 


•  Objective 

•  Background 

•  Overall  Program  Requirements 

•  Advantages 

•  Technical  Approach 

•  Status 

•  Summary 

•  Future  Plans 
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Caseless  Program  Objectives 

•  Characterize  the  HITP 

•  Demonstrate  a  function  and  production 
capability 

•  Identify  &  address  technical  challenges/potential 
risk  areas 

•  Deliver  prototype  Caseless  ammunition  for 
ballistic  demonstration  in  support  of  the 
Lightweight  Small  Arms  Technologies  (LSAT) 
Defense  Technology  Objective 

•  5.56mnn  cartridge  configuration 

•  Transfer  technology  to  industry 
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Background 


•  Original  work  performed  under  the 
Advanced  Combat  Rifle  (ACR)  Program 

•  Technology  Development  funded  by  US 
(ARDEC)  and  Germany  to  Heckler  & 
Koch(H&K)/Dynamit  Nobel(DNAG) 

•  Successful  Demonstration  of  a  Caseless 
Ammunition  Rifle  System 

•  Technology  Licensed  &  Transferred  to  the 
US  at  ARDEC 
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G-11  Open  Source  Data 


COMPONENTS  OF  A 
4.73  X  33MM  CASELESS  ROUND 


AMMO 


PRIMER 


PLASTIC 

CAP 


PROJECTILE 


PROPELLENT 
BODY 


Caseless  ammunition 


Length 

Cross-section 
Total  weight 
Projectife  weight 


33  mm/' 
8x8  mrr, 
5.20  g/0, 
3.25  g/0.1 


in. 


Ignition 

Mean  gas  pressure 
Muzzle  velocity  VD 


mechanical 
3850  bar 

approx.  93D  m/sec, 
3051  ft. /sec. 
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Overall  Program  Requirements 


Threshold  requirements 

•  25%  Decrease  in  ammunition  weight 

•  Same  lethality  as  the  5.56mm  M855  cartridge 

•  Environmentally  friendly  alternatives  for  processing 

•  Low  life-cycle  costs 


Desired  requirements 

•  Additional  decrease  in  ammunition  weight 

•  Increased  lethality  over  the  5.56mm  M855  cartridge 

•  Same  cost/round  as  current  5.56mm  M855 
ammunition 
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Why  Caseless  Ammunition? 


•  Lightweight 

-  Force  Multiplier 

-  Decreased  Logistics  Burden 

•  High  Ignition  Temperature 
Propellant  (HITP)  Provides 
Improved  Propellant 
Characteristics  &  Energetic 
Behavior 


•  Technology  with  other  potential 
applications 
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Weight  Reduction  Potential 


•  The  brass  cases  account  for  50  %  of  total  ammunition  weight! 

•  Current  Ammo  Weight  (600  rds):  20.8  lbs 

•  Target  Ammo  Weight  (600  rds):  10.2  lbs  «  50%  Weight  Reduction 


Why  Use  an  HUP? 


•  Advantages  of  Brass  Casing 

-  Barrel  &  chamber  temperature: 

•  Brass  is  a  heat  sink  that  is  discharged  with  each  round  fired 

-  Brass  has  a  long  history  in  small  caliber  ammunition 

-  Brass  provides  structural  strength  to  ammunition 

-  Brass  provides  gas  seal  in  weapon 


•  Advantages  of  HITP 

-  Weight  savings  with  elimination  of  the  brass  case 

•  Potential  for  high  firing  rates  with  the  elimination  of  the  case 
ejection  cycle 


-  High  thermal  stability  over  typical  ball  powder 

-  Structural  integrity  over  a  wide  temperature  range 
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Initial  Caseless  Development 


ARDEC 
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Caseless  Program  Status 


•  Overall  ARDEC  In-house  Progress 

•  HITP  formulation  verified  by  ARDEC/ATK 

•  Viable  material  sources  identified 

•  Producibility  demonstrated 

-  Several  hundred  rounds  have  been  produced  from 
lab-scale  propellant  mixes 

•  Conducted  three  ballistic  firings 

•  Continue  manufacturing  processes  and  prototyping 
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Caseless  Program  Status 


•  Conducted  three  ballistic  firings  in  Mann  Barrel 
at  Armament  Technology  Facility  at  ARDEC 

•  Test  Firing  I 

-  Five  Shots 

-  Chamber  Pressure  (8190-11314  psi) 

-  Muzzle  Velocity  (589-1204  fps) 


*  Test  Firing  II  (Addition  of  Booster  Charge) 

-  12  Shots 

-  Chamber  Pressure  (14372-62225  psi) 

-  Muzzle  Velocity  (1489-2795  fps) 

•  Test  Firing  III  JP"' 

-  14  Shots 

-  Chamber  Pressure  (15992-24579  psi) 

-  Muzzle  Velocity  (1000-1 686  fps) 
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Caseless  Program  Status 


•  Material 
Characterization 

•  Cooperative  effort 
between  ARDEC  &  ATK 

•  Original  &  New  HUP 

-  Chemical  Analysis 

-  Density 

-  Thermal 
Properties/Ignition 
Temperature 

-  Heat  of  Explosion 
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Caseless  Program  Status 


•  HITP  Characterization 

-  Differential  Scanning  Calorimetry  -  Thermal  Properties  of  HITP  in  comparison  to 
Nitrocellulose  based  Propellant 


Formulation  &  Process 
Development 


ATK  Thiokol 


ARES  he 


ATK 


Balteiie 


St.  Marks  Powder 

A  GENERAL  DYNAMICS  COMPANY 


Caseless  Program  Status 

•  Overall  Formulation  &  Process  Development 

Progress 

•  HITP  Characterization 

-  Formulation  verified 

•  Material  sources  identified  -  CONUS  and  OCONUS 

•  Raw  materials  procured  or  synthesized 

•  Compatibilities  verified  between  ingredients 

•  Hazard  Analysis 

•  HITP  Mixing  &  Processing 

-  Identifying  manufacturing  processes  &  limitations 

•  Closed  Bomb  Testing  and  Modeling 


Inc 


ATK 


St.  Marks  Powder 

A  GENERAL  DYNAMICS  COMPANY 


Caseless  Program  Status 


•  Hazard  Analysis 

•  Safety  Data  Determined 

•  Propellant  ingredients 

•  HUP 


ABL  Impact 
(cm) 

ABL  Friction 
(lbs  @  fps) 

TC  ESD 

(J) 

>80 

>800  @  8 

>8 
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Caseless  Program  Status 


•  HITP  Characterization-  Mixing 


•  Four  HITP  hand  mixes  produced 

•  Safety  testing/Lab  analysis 


-  Original  HITP  re-processed  for  lessons-learned 

-  Multiple  sub-scale  HITP  mixes  in  progress-Design  of 
Experiment  (  DOE) 


Caseless  Program  Status 


•  HITP  Processing 

-  HITP  Rheology  and  molding  studies 

•  Material  characteristics  and  processability  testing  with 
Capillary  Rheometer 


AJU  a 


Capillary  Rheometer 
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Burn  R  ate  (in/sec) 


Caseless  Program  Status 


•  HITP  Testing 

-  Low  &  High  Pressure  Bomb  analysis 
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Caseless  Program  Status 


•  HITP  Modeling 

-  Initial  closed  bomb  testing  gave 
repeatable  but  unexpected 
results 

•  The  CB  data  from  the  HITP  tests 
show  the  vivacity  plots  to  be 
inconsistent  with  regressive 
burnback 


P/P  max 


Typical  Regressive  Vivacity  Plot 


•  Grains  did  not  appear  to  burn  with 
the  expected  type  form  function 

•  Additional  testing  is  planned  to 
understand  this  issue 


ARES  he 


ATK 


Baiteiie 


Obtained  Progressive  Vivacity  Plot 
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Caseless  Program  Status 


•  HITP  Modeling  -  cont 

-  JHU/APL  Status 

•  ACR  1 B  sequence  recreated 

•  Modeled  as  single  perf  grain 

•  HITP  burn  data 

-  Reasonable  results  given  limited  data 

•  Pressure  &  velocity  close 

•  Sensitivity  analysis  of  several  parameters  generally 
consistent 

-  Model  will  be  updated  as  additional  data  becomes 
available 
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HITP  Prcpellart 
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Caseless  Program  Status 


•  Future  Plans 

-  Complete  HITP  Development  DOE 

•  Identify  critical  formulation  and  process  parameters 

-  Produce  first  prototype  Spring  2006 

•  Begin  Mann  barrel  testing  to  screen  DOE  mixes 

-  1st  Delivery  of  test  articles  for  validation  testing  by  Summer  of 
2006 

•  Primary  goal  is  to  deliver  cartridges  for  proof-of-concept 

-  Secure  material  sources  for  process  scale-up  work 

•  Some  materials  require  a  6-month  lead  time 

-  Support  follow-on  process  improvements 
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Caseless  Program  Summary 


•  Technical 


-  Initial  mixing  and  processing  very  promising 

-  Fabrication  solutions  appear  viable 

-  Interior  ballistic  modeling  is  underway  that  is  critical  to 
cartridge  design  -  can  be  accomplished  with 
additional  testing 

•  Schedule 


-  Mix  process  DOE  has  been  prepared  and  set  for 
completion  in  Spring  2006 

-  Smaller  DOE  will  likely  follow  to  finalize 
process/formulation  parameters  -  Summer  2006 

-  Cartridge  fabrication  for  testing  and  demonstration  - 
Summer  2006 
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Supercavitation 


^  Stresses  exerted  on  ambient  fluid 
drops  below  a  critical  value 

^  A  gaseous  cavity  envelopes  the  body 
and  moves  with  it 

^  Viscous  interaction  between 
projectile  surface  and  water  is 
alleviated  (i.e.  drag  reduction) 

^  Large  underwater  speeds  possible 
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Supercavitation  Inception 


^  Pure  Hydrodynamic  Cavitation  - 

Naturally  occurring  nuclei  (small  gas 
bubbles)  explosively  grow  due  to 
fluctuating  pressure  field  in  the 
separated  flow  region. 

Extremely  High  Speeds  Required  For 
Pure  Hydrodynamic  Supercavitation 


&  Artificial  (Ventilated)  Cavitation  - 

Pressurized  gas  is  injected  behind  the 
cavitator.  Supercavitation  at  lower 
speeds  is  possible. 

Requires  a  Pressurized  Source  of  Gas  to 
be  Carried  On-Board  or  Plumbing  to 
Re-Direct  Exhaust  Gases 
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Inception  via  Gas  Entrainment 


How  Can  We  Achieve 
the  Best  of  Both 
Worlds? 


^  Pressurized  gas  is  injected  during 

launch 


^  The  pressurized  gas  fills  the 
separated  flow  region,  resulting  in  a 
fully  developed  cavity 


^  Velocity  required  to  initiate 
supercavitation  in  drastically  reduced 

Allows  a  state  of  supercavitation 
to  be  prematurely  attained  until 
pure  hydrodynamic 
supercavitation  can  ensue 


(not  accounting  for  injected  gas) 
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Gas  Gun  Method 


^  The  projectile  leaves  the  launcher  under  the  force  of  expanding  gases  (i.e.  gun 
launch) 

^  The  driving  gases  are  the  source  of  cavitation  nuclei 

^  Degree  of  supercavitation  initiation  is  aspect  ratio  dependent  -  reattachment 
point  is  important 
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Vortex- Body  Interaction 


The  formation  of  a 
vortex  ring. 


Cavitating 
Vortex  Ring 


Generation  of  a  Interaction  of  the 

cavitating  vortex  vortex  ring  and 

r'n9-  projectile. 
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Gas  Gun  Method 


Blunt  Cylindrical 
Projectile 

Aspect  Ratio  5 

U0  =  27  m/s 

a  =  0.29 

(not  accounting  for 
injected  gas) 


A  large  degree  of 
fluid  structure 
interaction  takes 
place 

Asymmetric 
vortices  can  cause 
instability  in 
projectile  flight 
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Gas  Injection  Method 


a?  Pressurized  gas  is  injected  in  the  trajectory  of 
the  projectile  prior  to  launch 

^  The  supercavity  is  fully  developed  before  the 
projectile  interacts  with  the  surrounding  fluid 


U0  =  12  m/s 

a  =  1.48  (not  accounting 
for  injected  gas) 
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Gas  Injection  Effect 


U0  =  12.8  m/s 

a  =  1.48 

(not  accounting  for 

injected  gas) 

U0  =  13.8  m/s 
ct  =  1.12 


Gas  Injection  Method 
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Gas  Injection  Method 


Projectile  Trajectories 


Average  Projectile  Ending 
Location 


Higher  Gas  Injection  Pressures  Results  in  Stronger 
Cavity  Development  and  Better  Projectile  Stability 
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Conclusions 


<z?The  gas  injection  method  allows  the  projectile  to 
reach  a  state  of  “instant  supercavitation,”  thereby 
mitigating  any  fluid  structure  that  could  lead  to  its 
instability  during  launch 

&The  gas  injection  method  allows  supercavitation 
at  greatly  reduced  speeds 

^The  gas  injection  method  increases  projectile 
lateral  stability  and  decreases  viscous  drag 
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Smart  Material  Canard  Actuator  Design 


Actuator  Cliaractenstics 

•  Hollow  NACA  0018  aeroshell  surrounds 
actuator  mechanism. 

•  Actuator  Mechanism  consists  of  two  Macro 
Fiber  Composite  (MFC)  patches  bonded  to  e- 
glass  plate. 

•  Differential  voltage  applied  to  MFC 
patches  produces  deflection  at  trailing 
edge  of  plate  and  rotation  of  aeroshell 


Benefits  of  actuator  concept 

•  Minimizes  volume  intrusion  of  actuator  into 
munition  payload 

•  Potential  weight  savings 
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Smart  Material  Canard  Actuator  Design 
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Macro  Fiber  Composite  Patches 


•  Active  layer  consists  of  Macro  Fiber 
Composites  (FMC)  produced  by  Smart 
Material  Corp. 

•  FMCs  utilize  uniaxially  aligned  fibers 
surrounded  by  a  polymeric  matrix. 

•  Interdigitated  electrode  pattern  is  used  to 
deliver  an  electric  field  along  the  length  of 
the  fiber. 


Interdigitated 
electrode  pattern 
on  polyimide  film 
{top  and  bottom) 


Structural  epoxy 


Sheet  of  aligned 
rectangular 
piezoceramic 
fibers 


Extruded 

(cylindrical) 

piezoceramic 

fibers 


Fibers  encased  in 
otective  epoxy  matrix 


Interdigitated 
electrodes  on 
packaging  film 
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Design  Approach 


•  Multi-disciplinary  design  approach  -  coupled  structures  and 
aerodynamics. 

•  In  general,  there  is  a  trade-off  between  fin  deflection  angle  and 
available  torque  to  overcome  aerodynamic  hinge  moment. 

•  ARL/University  of  Delaware  design  -  maximize  fin  deflection  angle 
Mach  0.3-0. 5. 

•  It  is  possible  to  obtain  more  deflection  if  aerodynamic  torque  is 
ignored. 
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Applied  Torque  -  (N-mm) 


Design  Approach  (Cont.) 


Deflection  Angle  -  (Degrees) 
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Applied  Torque  -  (N-mm) 


Trade-Off  Between  Maximum  Deflection 
Angle  and  Torsional  Stiffness 


"1 
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Applied  Torque  -  (N-mm) 


4Qo3T& 


Deflection  Angles  in  Presence  of 
Aerodynamic  Torque 
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Deflection  Angle  -  (Degrees) 
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Wind  Tunnel  Testing  at  U.  of  Maryland 


•  Sponsored  by  ARL,  U.  Delaware 

•  Fin  mounted  to  balance  on  tunnel  floor  -  half  missile  body  fairing  used. 

•  Purpose  of  Wind  Tunnel  Test: 

•  Determine  response  of  fin  in  presence  of  aerodynamic  load. 

•  Measure  and  quantify  aerodynamic  loads. 

•  Determine  whether  flutter  is  an  issue. 


9 


Wind  Tunnel  Testing  at  Texas  A&M 


•  Collaboration  between  AMRDEC  Aero.  Branch,  ARL,  U.  Delaware 

•  Testing  at  Texas  A&M  University 

•  AMRDEC  provided  splitter  plate,  force/moment  balance,  testing  time. 

•  Purpose  of  Wind  Tunnel  Test: 

•  Determine  response  of  fin  in  presence  of  aerodynamic  load 

•  Measure  and  quantify  aerodynamic  loads 

•  Determine  whether  flutter  is  an  issue. 
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Optical  Encoder 


•  Measures  fin  deflection  angle 

•  Sensor  mounted  on  fin  base 

•  Sensor  reads  deflection  from  graduated  glass  disc  glued  to  bottom 
of  fin 

•  Time-dependent  sensor  output  acquired  on  laptop 
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Fin  Power  Sources 


•  Three  power  sources  provide  voltage  to  both  smart  material 
patches 

•  Power  controlled  through  laptop 

•  Continuous  range  of  voltage  available  -500V  to  +1500V 
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Wind-Off  Response 


•  Several  wind-off  tests  were  run  to  obtain  the  response 
of  the  fin  without  aerodynamic  load 

•  Results  have  some  bearing  on  experimental  approach 
and  interpretation  of  results 

•  Static  Deflection  Tests 


Peak-to-Peak 

Deflection 

Torsional  Stiffness 

250°  F  Fin 

5.3  deg 

179  N-mm/deg 

Loctite  Fin 

5.5  deg 

279  N-mm/deg 

•  Hysteresis 

•  Time-Dependent  Response 
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Encoder  Angle  Reading  (deg) 


Wind-Off  Hysteresis 


Delta  Voltage  Setting  (PS1  -  PS2) 
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Deflection  Angle  -  Deg 


Time-Dependent  Response 
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Rigid  Fin  Tests 


"1 


•  The  pure  aerodynamic  response  of  the  fin  was 
examined  using  a  rigid  fin  (no  actuator): 

•  Aerodynamic  forces  and  moments  measured  to  ±25 
deg  angle  of  attack  (AoA). 

•  “Hinge”  line  located  at  20%  of  chord  length  from  leading 
edge. 

•  Small  difference  between  rigid  fin  and  fins  with  acutator 

•  Fin  was  mounted  %”  inch  higher  off  splitter  plate 
than  for  fins  with  actuator. 

•  No  optical  encoder  mounted. 

•  Considered  two  velocities  1 80  ft/s  and  280  ft/s 

•  Measured  six  force/moment  components 


11/18/2005 


Aerodynamics  Branch 


16 


Hinge  Moment  Coeff 


Hinge  Moment 
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Smart  Material  Actuator  Fin  Tests 


•  The  smart  material  actuator  fin  was  tested  in  number  of 
modes: 

•  Unpowered  (No  applied  voltage) 

•  ±  Maximum  Voltage 

•  Varying  Voltage 

•  Tests  performed  at  180ft/s,  280ft/s 

•  AoA-14°  to  14° 

•  Two  versions  of  actuator  tested 

•  Differ  only  by  adhesive  used  to  bond  patch  to  host 

•  250°F  epoxy 

•  Loctite  epoxy 

•  Force/moment  balance  measurements  taken. 

•  Fin  deflection  measured  with  optical  encoder. 
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Hinge  Moment  Coeff 


Hinge  Moments 
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Lift  Coeff 


Lift  Force  Generated  by 
Deflection 
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Fin  Deflection  Angle  vs.  Angle  of 

Attack 
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Fin  Deflection  Angle  (Degrees) 


Fin  Deflection  Angle  vs.  Angle  of 

Attack 
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Conclusions 


"1 


•  Within  one  year,  a  canard  actuator  based  on  smart 
materials  was  designed,  constructed  and  wind-tunnel 
tested. 

•  Design  shows  promise  for  application  to  munition 
systems. 

•  Coupled  structural/aerodynamic  analysis  required 
for  the  design  of  these  flexible  structures. 

•  MFC  actuators:  rugged  and  robust  -  a  requirement 
for  gun-launched  munitions. 

•  Flutter  not  an  issue  for  velocity  range  tested. 

•  Improved  performance  possible  with  existing  design 
using  better  manufacturing  techniques. 

•  Alternative  designs,  improved  MFC  patches  being 
investigated. 


11/18/2005 


Aerodynamics  Branch 


23 


The  Study  on  Lethality  Simulation 
Method  for  Fragmentation  Warhead 


Yang  Yunbin  Qu  Ming  Qian  Lixin 


4^- 


10 go  Institute  of  Structural  Mechanics 


November,  2005 


Table  of  contents 


■4r 


EQ  INTRODUCTION 


CQ  LETHALITY  SIMULATION  METHOD 


CQ  EXAMPLE  ANALYSIS 


CQ  CONCLUSION 


CQ  FUTURE  WORK 


CQ  THANKS 


*  ISB2005 


2th  Page 


INTRODUCTION 


*r 


^  Warhead  lethality 
assessment  is  an 
important  factor  of 
warhead  study.  It 
always  includes 
warhead  lethality 
description,  missile- 
target  intersection 
and  vulnerability. 


A  work  flow  of  warhead  study 
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INTRODUCTION 


&  With  the  development  of  numerical  simulation, 
many  new  methods  of  lethality  assessment  come 
forth 

Lethality  test  simulation  could  examine  whether 
lethality  target  data  is  satisfied  with  design 
requirement  before  making  the  product,  so 
development  cycle  is  shorten  and  outlay  is  reduced. 

SSPTest  method,  analytical  method  and  numerical 
simulation  method  are  widely  adopted  . 
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Analytical  method 


Analytical  methods  are  often  used.  Such  as  SplitX  and 
LATFW,  they  have  been  applied  largely  in  design  and 
assessment  of  warhead . 

r  SplitX  of  CONDAT  Company  is  an  expert  system  for 
the  design  of  fragmentation  warhead. 

>  Assess  warhead  by 
analytical  method 


ISB2005 


(a)  Fragment  Velocity  Distribution 
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Analytical  method 


>  LATFW  (Lethality  Assessment  of  Typical 
Fragmentation  Warhead)  is  achieved  by  our  institute 
in  China.  \  V.  11  I  n 


(b)  damage  effect  of  one  cabin  (c)  Shotline  simulation 
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Analytical  method 


&  Virtue 


>  It  can  be  used  to  rapidly  assess  a  new  or  existing 
warhead  design. 

>Warhead  design  is  accomplished  based  on  the 
method  in  a  short  space  of  time. 


Shortcoming 

>  The  deviation  can  be  introduced  if  the  methods  do 
not  be  verified  or  modified  by  test  data  . 

>The  methods  are  limited  by  specifically  condition. 
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Numerical  simulation  method 


The  prediction  results  by  the  numerical  simulation 
method  are  more  accurate  than  those  by  forenamed 
methods.  It  can  make  structural  analysis,  detonation 
calculation  and  so  on. 

The  method  is  fit  for  the  multifarious  warheads,  such  as 
the  axisymmetric  warhead  and  unsymmetrical  warhead. 

Thinking  over  calculation  scale  and  complication.  The 
method  is  adopted  in  one  stage  of  lethality  assessment 
process. 
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Test  method 


Early,  the  lethality  assessment  method  is  warhead  test. 
However,  a  large  number  of  lethality  tests  may  prolong 
development  cycle  and  increase  outlay. 

In  present,  the  integrated  method  including  test 
method,  analytical  method  and  numerical  simulation 
method  is  commonly  adopted. 
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Our  Study 


s3rHow  do  we  make  lethality  simulation  for  fragmentation 
warhead  when  efficiency  and  precision  must  be 
considered? 

ssrBy  integrating  Numerical  simulation  method  and 
analytical  method,  lethality  simulation  method  for 
fragmentation  warhead  is  established. 

•^By  numerical  simulation  method,  the  fragment  initial 
fields  are  gained.  Based  on  analytical  method,  fragment 
movement  model  and  fragment  impact  model  on  target  are 
established. 
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Our  Study 


The  whole  process  description  of  fragment  field 
formed,  fragment  movement,  and  fragment  impacting  on 
target  is  achieved  by  using  the  method. 

^  In  my  topic,  I  only  present  my  outlook  of  our  study.  I 
want  to  validate  whether  the  method  is  feasibility.  We 
analyze  the  powerful  parameters  that  affect  warhead 
lethality  based  on  example  analysis. 
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LETALITY  SIMULATION  METHOD 


Lethality  simulation  method  of  fragmentation 
warhead  consists  of  numerical  simulation  model, 
fragment  movement  model,  lethality  parameter 
model,  damage  analysis  model  and  fragment  field 
simulation  model. 

>  Fragment  initial  fields  are  gained  based  on 
numerical  results. 

>  Fragment  movement  model  calculates  fragment 
trajectory.  Damage  analysis  model  analyzes  fragment 
damage  performance  to  target. 
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LETALITY  SIMULATION  METHOD 


>  Lethality  parameter  model  analyzes  fragment 
hitting  density  distribution  and  the  dispersal  angle 
distribution  of  fragments. 

>  Fragment  field  formed,  fragment  movement, 
and  fragment  impacting  on  target  are  achieved  by 
fragment  field  simulation  model. 
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Numerical  simulation  model 


>  Numerical  simulation  model  includes  numerical 
simulation  and  Interface  middleware.  Numerical 
simulation  is  accomplished  by  LS-DYNA.  Because  LS- 
DYNA  only  outputs  node  velocity  and  location,  interface 
middleware  is  developed  to  achieve  fragment  velocity 
and  location.  Using  interface  middleware,  fragment  initial 
field  is  gained  and  saved. 


Fragment 

No 

Location 

Velocity 

X 

m 

Y 

m 

Z 

in 

Vx 

m/s 

Vy 

m/s 

Vz 

m/s 

1 

0.0011 

0.6993 

0.01557 

29.68 

2627.00 

36.60 

i 

0.0089 

0.7145 

0.07378 

34.43 

2723.90 

27.81 
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Numerical  simulation  model 


3  Interface  middleware  is  builded  based  on 
modified  former  code.  Numerical  results  are  dealed 
through  calling  interface  middleware,  so  data 
exchange  between  numerical  simulation  model  with 
fragment  movement  model. 
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Fragment  movement  model 


3U  In  fragment  movement  stage,  fragment  characteristic 
parameters  are  calculated  in  fragment  field.  Because  of 
small  fragment  mass  and  short  fragment  trajectory  as  well 
as  high  fragment  velocity,  fragment  gravitational  effect  is 
ignored  and  we  assume  fragment  trajectory  is  linear. 


[5  Assuming  the  air  drag  coefficient  is  invariable, 
fragment  residual  velocity  is  calculated  by 


VT  =V0exp  (-(Cvp»ffV)A'g)r) 

2  Q 


c„  =  + 


S  -S 

n  nsphere 


(^Dnature  sphere  ) 


'D  D  sphere  ^  v  Dnature  D  sphere , 

J  nnature  J  nsphere 
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Fragment  movement  model 


H  In  this  analytical  model,  we  don’t  think  over  the  effect 
of  fragment  tumbling. 

3  In  the  paper  (Fragment  Shot-line  Model  for  Air- 
Defense  Warhead,  PEP  25,  No.2  2000),  the  analytical 
model  including  the  effect  of  fragment  tumbling  and  air 
drag  on  fragment  trajectory  is  presented  by  our  team. 
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Damage  analysis  model 


3H  THOR  method  is  adopted  to  assessment  the  maximum 
penetration  thickness,  the  residual  mass  and  the  residual 
velocity. 

>THOR  fits  for  steel  fragment  impacting  on  the  metallic 
and  nonmetallic  target,  especially  it  fits  for  fragment  of 
high  detonator  warhead 

Vr  =  Vs -0.3048  xlOCn(61O23.75/iA)c,2(l5432.1mjyu(sec0yH(3.28O84F,yt5 

< 

mr  =  ms  -6.48xl0C51(61023 . 75/l4 !;(15432. lmj^^sece1  )r;4(3. 28084 Fj"25 


/;max  =  1.638706  x  10‘~ll/<:i2 .4“1(l5432.lOTi  )“C‘3/Cl1  (sec0)_Cl4/c“  (3 ,28084F5 )(l‘ 


'U 


)A 


12 
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Damage  analysis  model 


Saving  damage  assessment  parameter 


Work  flow  of  damage  assessment 
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Lethality  parameter  model 


•*“  Definitions  of  the  variables  I 

W  L 


W — Width  of  statistical  area  L — -Length  of  statistical  area. 
0X  — -Fragment  projection  angle  ft,  — -W arhead  radius  angle. 

D — Target  distance. 
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Lethality  parameter  model 


3  Fragment  density  is  affected  by  statistical  area.  In 
order  that  statistical  results  are  accurate  the  length  of 
statistical  area  is  ascertained  according  to  warhead  radius 
angle.  According  to  the  dispersal  angle  of  fragments  the 
width  of  statistical  area  is  calculated,  the  work  flow: 


(2)  L  =2  D  ■  tan(^-) 


(1)  6.= 2  arctan(  — ) 

ii-// 


(4)  Based  on  warhead  radius  angle  the  statistical  number  is 
ascertained,  for  the  whole  target.  Fragment  density  of  each 
statistical  area  is  gained  by  the  step(l) — (3) 
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Simulation  model 


System  simulation  model  divides  into  three  stages. 


System  simulation  model 
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Simulation  model 


H  System  simulation  model  calculates  trajectory, 
velocity  ,  residual  velocity  and  residual  mass. 

3  In  order  to  investigate  fragment  field  in  different  time 
simulation  data  is  saved,  and  fragment  field  simulation  is 
displayed  in  other  CAD  software. 
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EXAMPLE  ANALYSIS 


On  the  basis  of  the  simulation  method,  lethality 
simulation  code  is  developed.  One  example  of  an 
aimable  warhead  is  examined. 

t  System  visualization 
t  Analysis  of  fragment  density 
#  Analysis  of  fragment  projection  angle 
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System  visualization 


Fragment  field  simulation 


(a)  t  =  0.25  ms  (b)  x  =  0.46  ms 


ISB2005 
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System  visualization 


F ragment  field  simulation 


(a)  0.25  ms  (b)  1.98  ms  (c)  2.88  ms 

>The  aimable  area  is  not  obvious  at  0.25  ms.  The  f 
velocity  in  aimable  area  is  higher  than  those  in  non-aimable 
area  at  2.88  ms.  Aimable  fragments  firstly  hit  target. 


26th  Page 


ISB2005 


System  visualization 


“^Fragment  field  simulation 


>This  graph  shows  visualization  graph  of  fragment 
distribution  on  target.  Target  height  is  10m  and  target 
distance  is  20m.  The  circle  target  is  unwrapped  in 
order  to  represent  fragment  distribution  clearly. 
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Analysis  of  fragment  density 


Two  cases 


02  =  30  °  02  =9(7 

>At  30°  fragment  density  fluctuates  heavily  in  each 
statistical  area,  The  aimable  area  is  not  obvious. 

>At  90°  the  aimable  area  is  obvious. 

28th  Page 


ISB2005 


Analysis  of  fragment  density 


;s  Results 

>  According  to  the  analytical  results,  the  variational 
trend  of  fragment  density  is  similar  when  02  is  same 
even  if  the  target  distance  is  different. 

>  Suppose  fragment  trajectory  is  linear,  when  02  is 
same  fragment  density  is  in  inverse  proportion  to  the 
square  of  target  distance  that  is  match  with  analytical 
results. 
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Analysis  of  fragment  projection  angle 


“^Considering  the  relationship  between  fragment 
projection  angle  and  warhead  radius  angle,  fragment 
projection  angle  corresponding  is  the  average  value 
of  projection  angle  in  different  statistical  area. 

“^Fragment  projection  angle  is  not  sensitive  to 
02  and  the  statistical  error  is  less  than  1 0 .  Fragment 
projection  angle  ranges  from  42 0  to  43 0 . 
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Analysis  of  fragment  projection  angle 


Results 

> Fragment  projection  angle  in  different  statistical 
area  fluctuates  when  02  is  same. 

>At  the  same  time  fragment  project  angle  fluctuates 
with  the  change  of  warhead  radius  angle  but  the 
extent  is  small. 

>Fragment  project  angle  distribution  of  aimable 
warhead  is  unsymmetrical  and  fluctuates  heavily. 
The  factor  must  be  considered.  It’s  better  to  point  out 
which  area  projection  angel  is  in  aimable  area,  non- 
aimable  area  or  all  the  area. 
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CONCLUSION 


Lethality  simulation  method  includes  numerical 
simulation  and  theoretical  analysis,  so  the  method  fits  for 
different  kinds  of  the  warhead  by  generalizing  it. 

Based  on  the  method,  lethality  simulation  system  is 
constructed,  and  numerical  simulation  of  fragmentation 
warhead  lethality  are  achieved  by  the  system,  and 
analytical  results  are  directly  applicable  for  warhead 
lethality  assessment. 

^  Based  on  the  test,  the  method  is  an  effective  method 
of  validating  numerical  simulation. 
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FUTURE  WORK 


•*"  How  to  assess  Virtual  Prototyping  (Numerical  Products)  of 
warhead  thinking  over  precision  and  efficiency? 

How  to  integrate  lethality  assessment  and  warhead  design, 
and  establish  the  compositive  design  environment  of  warhead 
study. 

Our  thinking 

>  Lethality  simulation  method  and  code  are  studied 
by  integrating  analytical  method,  numerical  simulation 
method  in  warhead  lethality  description,  missile-target 
intersection  and  vulnerability. 
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FUTURE  WORK 


Warhead  lethality  description 

>  Warhead  lethality  description  for  different  kinds  of 
the  warhead  is  studied  by  the  integrated  method. 

•^Missile-target  intersection 

>  Static  intersection  ->•  Dynamic  intersection. 
•^Vulnerability 

> Vulnerability  analysis  method  is  traced  and  analyzed. 

>  Structural  damage  and  functional  damage. 
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FUTURE  WORK 


•^Lethality  Simulation  Code 


>This  code  is  easy  to  integrate  with  the  compositive 
design  environment. 
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